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NOMENCLATURE 


Symbols 
Main  variables 


A 

Ac 

A, 

A2 

A3 

An 

A5 

B 

b 

C 

Cc 

CD 

CD,e 

Cr 

c, 

c2 


area 

area  of  vena  contracta 

upstream  flow  area  or  area  at  position  1 

downstream  flow  area  or  area  at  position  2 

area  at  position  3 

area  at  position  4 

area  at  position  5 

diameter  of  cylinder  wall 

width  of  piston  ring 

radial  deflection  of  piston  ring 

coefficient  of  contraction 

discharge  coefficient  for  sudden  contraction 

discharge  coefficient  for  sudden  enlargement 

crown  height  of  piston  ring 

constant  or  sonic  velocity  in  liquid  phase 

constant 


XIV 


Main  variables  (cont.) 

c3 

constant 

COP 

coefficient  of  performance 

COPB 

COP  of  base  cycle 

COPe 

COP  of  Economizer  cycle 

COPed-i 

COP  based  on  expansion  through  the  ED-WOW  (by  considering  only 
lower  enthalpy  effect) 

COPed-2 

COP  based  on  expansion  through  the  ED-WOW  (by  considering  both 
of  lower  enthalpy  and  work  output  effects) 

COPexv 

COP  based  on  isenthalpic  expansion 

COP, 

COP  of  Intercooler  cycle 

D 

diameter 

De 

equivalent  diameter 

Dexh 

orifice  diameter  of  exhaust  valve 

Dfit 

inner  diameter  of  fittings 

Dint 

orifice  diameter  of  intake  valve 

e 

surface  roughness 

E 

Young’s  modulus  of  piston  ring  or  energy 

Ecv 

energy  in  a  control  volume 

F 

friction  force  or  function 

F mass 

force  due  to  the  mass 

F scale 

force  measured  by  spring  scale 

F' 

derivative  of  function  F 

f 

friction  coefficient 

Main  variables  (cont.) 


G 

Gc 

G\ 

G2 

h 

h-B\ 

hB2 

h-B3 

tlB4 

hcom.exi 

hcom.in 

h  comp,  in 

hcomp.out 

hcomp,out,s 

hED.i 

hED-WOW,in 
h  ED-WOW,  out 
hED-WOW,out,s 
hevap.i 

h 

rievap.o 

hexh 


mass  velocity 


mass  velocity  defined  as 


£2 

Cc 


upstream  mass  velocity 

downstream  mass  velocity 

film  thickness  of  lubricant  or  enthalpy 

enthalpy  at  the  inlet  to  the  compressor  in  the  Base  cycle 

enthalpy  at  the  inlet  to  the  gas  cooler  in  the  Base  cycle 

enthalpy  at  the  outlet  of  the  gas  cooler  in  the  Base  cycle 

enthalpy  at  the  inlet  to  the  evaporator  in  the  Base  cycle 

enthalpy  at  the  outlet  of  the  compressor 

enthalpy  at  the  inlet  to  the  compressor 

enthalpy  at  the  inlet  to  the  compressor 

enthalpy  at  the  outlet  of  the  compressor 

enthalpy  at  the  outlet  of  the  compressor  for  an  isentropic  process 
enthalpy  at  the  inlet  of  the  ED-WOW 
enthalpy  at  the  inlet  to  the  ED-WOW 
enthalpy  at  the  outlet  of  the  ED-WOW 

enthalpy  at  the  outlet  of  the  ED-WOW  for  an  isentropic  process 
enthalpy  at  the  inlet  of  the  evaporator 
enthalpy  at  the  outlet  of  the  evaporator 
enthalpy  at  the  outlet  of  the  expansion  device 


XVI 


Main  variables  (cont.) 


hE\ 

hEV 

Hei" 

flE2 

h£3 

hm 

/l£3” 

7l£4 

hE5 


enthalpy  at  the  inlet  to  the  lst-stage  compressor  in  the  Economizer 
cycle 

enthalpy  at  the  outlet  of  the  ls,-stage  compressor  in  the  Economizer 
cycle 

enthalpy  at  the  inlet  to  the  2nd-stage  compressor  in  the  Economizer 
cycle 

enthalpy  at  the  inlet  to  the  gas  cooler  in  the  Economizer  cycle 

enthalpy  at  the  outlet  of  the  gas  cooler  in  the  Economizer  cycle 

enthalpy  at  the  inlet  to  the  accumulator  in  the  Economizer  cycle 

enthalpy  of  liquid  phase  CO2  at  the  accumulator  in  the  Economizer 
cycle 

enthalpy  at  the  inlet  to  the  evaporator  in  the  Economizer  cycle 

enthalpy  of  vapor  phase  CO2  at  the  accumulator  in  the  Economizer 
cycle 


hin 

h/\ 

hiv 

ha 

hn 

hu 


enthalpy  at  the  inlet  to  the  intake  valve 

enthalpy  at  the  inlet  to  the  lsl-stage  compressor  in  the  Intercooler  cycle 

enthalpy  at  the  outlet  of  the  lst-stage  compressor  in  the  Intercooler 
cycle 

enthalpy  at  the  inlet  to  the  gas  cooler  in  the  Intercooler  cycle 
enthalpy  at  the  outlet  of  the  gas  cooler  in  the  Intercooler  cycle 
enthalpy  at  the  inlet  to  the  evaporator  in  the  Intercooler  cycle 


him 


head  loss 


hm  minimum  film  thickness  of  lubricant 

hml  minimum  film  thickness  of  lubricant  calculated  by  using  internal 

halving  method  for  comparison 


XVII 


Main  variables  (cont.) 


minimum  film  thickness  of  lubricant 


piston  ring  face  profile 


enthalpy  at  the  outlet  of  the  ED-WOW  for  the  isentropic  expansion 
process 

film  thickness  of  lubricant  far  away  from  slippery  bearing 
minimum  film  thickness  of  lubricant  at  (n-l)th  iteration  step 
minimum  film  thickness  of  lubricant  at  nth  iteration  step 
heat  removal  capacity  of  the  evaporator 

analytical  integration  term 
analytical  integration  term 
analytical  integration  term 
constant 

minor  loss  coefficient  for  sudden  contraction 
minor  loss  coefficient  for  sudden  enlargement 
loss  coefficient  for  exit 
length  of  pipe  or  length  of  connecting  rod 
mass  inside  chamber 


amount  of  mass  flowing  into  the  chamber 
amount  of  mass  flowing  out  from  the  chamber 
mass  inside  chamber  before  the  piston  movement 
mass  inside  chamber  after  the  piston  movement 


mass  flow  rate 


Main  variables  (cont.) 


W-COi-leak 


H-in 

Ylexh 

P 

Pbot 

Pbk 

PB\ 

PB2 
PB3 
Perk 
P  c,in 
P  c.out 
P  ED,i 
P  ED,  in 

Ped,o 
P  ED, out 
Pel 
pexh 


mass  flow  rate  of  CO2 
CO2  leakage  rate  through  the  piston  rings 
mass  flow  rate  into  the  chamber 
mass  flow  rate  from  the  chamber 

number  of  intake  valves 

number  of  exhaust  valves 

pressure  of  lubricant 

pressure  at  trailing  edge  of  piston  ring 

pressure  at  the  back  side  of  piston  ring 

pressure  at  the  inlet  to  the  compressor  in  the  Base  cycle 

pressure  at  the  inlet  to  the  gas  cooler  in  the  Base  cycle 

pressure  at  the  outlet  of  the  gas  cooler  in  the  Base  cycle 

pressure  inside  enclosure 

pressure  at  the  inlet  of  the  compressor 

pressure  at  the  outlet  of  the  compressor 

pressure  at  the  inlet  of  the  ED-WOW 

pressure  at  the  inlet  of  the  ED-WOW 

pressure  at  the  outlet  of  the  ED-WOW 

pressure  at  the  outlet  of  the  ED-WOW 

piston  ring  elastic  pressure 

pressure  at  the  outlet  of  the  expansion  device 


XIX 


Main  variables  (cont.) 

Pevap.o  pressure  at  the  out  of  the  evaporator 

Pevap.out  pressure  at  the  out  of  the  evaporator 

Pe\  pressure  at  the  inlet  to  the  lst-stage  compressor  in  the  Economizer 

cycle 

Pev  pressure  at  the  outlet  of  the  lst-stage  compressor  in  the  Economizer 

cycle 

Pei"  pressure  at  the  inlet  to  the  2nd-stage  compressor  in  the  Economizer 

cycle 

Pei  pressure  at  the  inlet  to  the  gas  cooler  in  the  Economizer  cycle 

Pe3  pressure  at  the  outlet  of  the  gas  cooler  in  the  Economizer  cycle 

Pev  pressure  at  the  inlet  to  the  accumulator  in  the  Economizer  cycle 

Pet  liquid  pressure  in  the  accumulator  in  the  Economizer  cycle 

Pea  pressure  at  the  inlet  to  the  evaporator  in  the  Economizer  cycle 

Pes  vapor  pressure  in  the  accumulator  in  the  Economizer  cycle 

Pgc,in  pressure  at  the  inlet  of  the  gas  cooler 

Pgc.out  pressure  at  the  outlet  of  the  gas  cooler 

Phigh  high-side  pressure 

Pin  pressure  at  the  inlet  to  the  intake  valve  or  inlet  pressure  of  the  ED- 

wow 

pn  pressure  at  the  inlet  to  the  lst-stage  compressor  in  the  Intercooler  cycle 

piv  pressure  at  the  outlet  of  the  lst-stage  compressor  in  the  Intercooler 

cycle 

P12  pressure  at  the  inlet  to  the  gas  cooler  in  the  Intercooler  cycle 

Pit  pressure  at  the  inlet  to  the  2nd-stage  compressor  in  the  Intercooler  cycle 


Main  variables  (cont.) 


P/3 
Pi 4 
Pm 
Pout 

Ppp 

Pp,in 

P p,out 

P  ss 

Ptop 

Ppump.i 

Ppump.o 


P\ 

P2 

P3 

P4 

P5 

Pi 

Pi 

pT 

pT 

Pi 


pressure  at  the  outlet  of  the  gas  cooler  in  the  Intercooler  cycle 

pressure  at  the  inlet  to  the  evaporator  in  the  Intercooler  cycle 

external  mean  hydrodynamic  pressure 

outlet  pressure  of  the  ED- WOW 

power  of  the  hydraulic  pump 

inlet  pressure  of  the  hydraulic  pump 

outlet  pressure  of  the  hydraulic  pump 

power  calculated  using  the  torque  of  the  ED-WOW 

pressure  at  leading  edge  of  piston  ring 

static  pressure  at  the  inlet  of  the  pump 

static  pressure  at  the  outlet  of  the  pump 

upstream  pressure  or  pressure  at  position  1 

downstream  pressure  or  pressure  at  position  2 

pressure  at  position  3 

pressure  at  position  4 

pressure  at  position  5 

pressure  at  position  1  at  the  nth  iteration 

pressure  at  position  2  at  the  nth  iteration 
pressure  at  position  2  at  the  (n-l)th  iteration 
pressure  at  position  2  at  the  (n-l)th  iteration 
pressure  at  position  5  at  the  nth  iteration 
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Main  variables  (cont.) 

p  uniform  internal  gas  pressure 


Q 

qmx 

Qed 

Qisenth 

qiHX.max 

Qe,B 

Qe,E 

Qe,I 

R 

r 

S 

Sb\ 

SED.i 

SEV 

SEV’ 

SE2 

SE3 

SE3' 


heat  loss  rate  or  volumetric  water  flow  rate 

heat  transfer  from  hot  CO2  to  cold  CO2  through  the  internal  heat 
exchanger  per  unit  mass 

evaporator  capacity  rate  with  expansion  through  the  ED-WOW 

evaporator  capacity  rate  with  isenthalpic  process 

maximum  heat  transfer  from  hot  CO2  to  cold  CO2  through  the  internal 
heat  exchanger  per  unit  mass 

heat  capacity  of  evaporator  in  the  Base  cycle 

heat  capacity  of  evaporator  in  the  Economizer  cycle 

heat  capacity  of  evaporator  in  the  Intercooler  cycle 

crank  radius 

radius  of  piston  ring,  or  cylinder  bore  radius,  or  radius  of  rod 
engine  speed 

entropy  at  the  inlet  to  the  compressor  in  the  Base  cycle 
entropy  at  the  inlet  of  the  ED-WOW 

entropy  at  the  outlet  of  the  lst-stage  compressor  in  the  Economizer 
cycle 

entropy  at  the  inlet  to  the  2nd-stage  compressor  in  the  Economizer 
cycle 

entropy  at  the  inlet  to  the  gas  cooler  in  the  Economizer  cycle 
entropy  at  the  outlet  of  the  gas  cooler  in  the  Economizer  cycle 
entropy  at  the  inlet  to  the  accumulator  in  the  Economizer  cycle 


xxn 


Main  variables  (cont.) 

Se3„  entropy  of  liquid  phase  C02  in  the  accumulator  in  the  Economizer 

cycle 

sE4  entropy  at  the  inlet  to  the  evaporator  in  the  Economizer  cycle 

Se5  entropy  of  vapor  phase  C02  at  the  accumulator  in  the  Economizer 


si  i 
sir 

Sl2 

Sl2' 

STD 

T 

Tair 

Tbi 

Tb2 

Tb3 

Tc,in 

TC,OUt 
Tedj 
Ted, in 

T evap 
Tevap.o 
7 evap, out 


cycle 

entropy  at  the  inlet  to  the  lst-stage  compressor  in  the  Intercooler  cycle 

entropy  at  the  outlet  of  the  lst-stage  compressor  in  the  Intercooler 
cycle 

entropy  at  the  inlet  to  the  gas  cooler  in  the  Intercooler  cycle 

entropy  at  the  inlet  to  the  2nd-stage  compressor  in  the  Intercooler  cycle 

standard  deviation 

tangential  ring  tension 

air  temperature 

temperature  at  the  inlet  to  the  compressor  in  the  Base  cycle 

temperature  at  the  inlet  to  the  gas  cooler  in  the  Base  cycle 

temperature  at  the  outlet  of  the  gas  cooler  in  the  Base  cycle 

temperature  at  the  inlet  of  the  compressor 

temperature  at  the  outlet  of  the  compressor 

temperature  at  the  inlet  of  the  ED-WOW 

temperature  at  the  inlet  of  the  ED-WOW 

evaporation  temperature 

temperature  at  the  outlet  of  the  evaporator 

temperature  at  the  outlet  of  the  evaporator 


Main  variables  (cont.) 


temperature  at  the  outlet  of  the  expansion  valve 

temperature  at  the  inlet  to  the  lst-stage  compressor  in  the  Economizer 
cycle 

temperature  at  the  outlet  of  the  lst-stage  compressor  in  the  Economizer 
cycle 

temperature  at  the  inlet  to  the  2nd-stage  compressor  in  the  Economizer 
cycle 

temperature  at  the  inlet  to  the  gas  cooler  in  the  Economizer  cycle 
temperature  at  the  outlet  of  the  gas  cooler  in  the  Economizer  cycle 
temperature  at  the  inlet  to  the  accumulator  in  the  Economizer  cycle 
liquid  temperature  in  the  accumulator  in  the  Economizer  cycle 
temperature  at  the  inlet  to  the  evaporator  in  the  Economizer  cycle 
vapor  temperature  in  the  accumulator  in  the  Economizer  cycle 
CO2  temperature  at  the  gas  cooler  outlet 
temperature  at  the  inlet  of  the  gas  cooler 
temperature  at  the  outlet  of  the  gas  cooler 

temperature  at  the  inlet  to  the  intake  valve  or  inlet  temperature  of  the 
ED-WOW 


T indoor 


indoor  temperature 

temperature  at  the  inlet  to  the  lst-stage  compressor  in  the  Intercooler 
cycle 

temperature  at  the  outlet  of  the  lst-stage  compressor  in  the  Intercooler 
cycle 


temperature  at  the  inlet  to  the  gas  cooler  in  the  Intercooler  cycle 
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Main  variables  (cont.) 


Tit 


Tn 

TnA 

Tn 
Tn+ 1 

Tout 

Toutdoor 


tr 

TSat 

r, 

t2 

t3 

r4 

u 

u, 


Ufl 

UfZ 

Ugl 

Ug2 


temperature  at  the  inlet  to  the  2nd-stage  compressor  in  the  Intercooler 
cycle 

temperature  at  the  outlet  of  the  gas  cooler  in  the  Intercooler  cycle 

temperature  at  (n-l)th  iteration  step 

temperature  at  n*  iteration  step 

temperature  at  (n+l)th  iteration  step 

outlet  temperature  of  the  ED- WOW 

outdoor  temperature 

torque 

piston  ring  thickness 
C02  saturation  temperature 
temperature  at  position  1 
temperature  at  position  2 
temperature  at  position  3 
temperature  at  position  4 
piston  speed 
impact  velocity 
x-component  of  velocity 
upstream  velocity  of  liquid  phase 
downstream  velocity  of  liquid  phase 
upstream  velocity  of  vapor  phase 
downstream  velocity  of  vapor  phase 
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Main  variables  (cont.) 


internal  energy  inside  chamber  before  the  piston  movement 

down  stream  velocity  or  internal  energy  inside  chamber  after  the 
piston  movement 

velocity  of  piston  ring  in  y-coordinate 
chamber  volume 

velocity  or  y-component  of  velocity 
volumetric  flow  rate  of  water 
chamber  volume  at  bottom  dead  center 
specific  volume  of  liquid  phase 

specific  volume  difference  between  vapor  and  liquid  phases 
specific  volume  of  vapor  phase 

volume  of  fitting  attached  to  the  intake  and  exhaust  valves 

chamber  volume  at  top  dead  center 

chamber  volume  before  piston  movement 

chamber  volume  after  piston  movement 

velocity  at  position  2 

specific  volume  of  liquid  downstream 

specific  volume  difference  between  vapor  and  liquid  phases 
downstream 

velocity  at  position  3 

velocity  at  position  4 


flow  velocity 


flow  velocity  at  position  1 


XXVI 


Main  variables  (cont.) 

v  flow  velocity  at  position  2 

W  load  imposed  on  piston  ring,  or  work  done  by  system,  or  mass  flow 

rate 

Vk  work  per  unit  time  done  by  CO2 

Wa  total  uncertainty  associated  with  the  dependent  variable  A 

Wcgm  compression  power 


'  comp 

WComp 

WComp,s 

WcOP 

WC,B 

WclE 

Wc2,E 

Wc 

WC2,, 

Wd 

WED 

WED-WOW 


compression  work 
specific  compression  work 

specific  compression  work  for  an  isentropic  process 
uncertainty  for  COP 

compression  work  in  the  Base  cycle  per  unit  mass 

compression  work  of  the  lst-stage  compressor  in  the  Economizer  cycle 

compression  work  of  the  2nd-stage  compressor  in  the  Economizer 
cycle 

compression  work  of  the  lst-stage  compressor  in  the  Intercooler  cycle 
per  unit  mass 

compression  work  of  the  2nd-stage  compressor  in  the  Intercooler  cycle 
per  unit  mass 

compression  work  with  different  pressure  ratio 
power  extracted  through  the  ED-WOW 
specific  work  extracted  through  the  ED-WOW 

specific  work  extracted  through  the  ED-WOW  with  isentropic  process 
mass  flow  rate  of  liquid  phase 


wed-wow, s 
Wf 
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Main  variables  (cont.) 

Wf  specific  compression  work  with  fixed  pressure  ratio 

Wg  mass  flow  rate  of  vapor  phase 

whED  uncertainty  for  the  enthalpy  at  the  inlet  of  the  ED-WOW 

wh  _  uncertainty  for  the  enthalpy  at  the  outlet  of  the  evaporator 

w n  uncertainty  for  the  CO2  mass  flow  rate 

W, iet  net  work  output  done  by  CO2 

w#  uncertainty  for  the  work  produced  through  the  ED-WOW 

wpn> .  uncertainty  for  the  pressure  at  the  inlet  of  the  ED-WOW 

wPed  uncertainty  for  the  pressure  at  the  outlet  of  the  ED-WOW 

w  uncertainty  for  pressure  at  the  inlet  of  the  pump 

Ppump.i  r  1 

wp  uncertainty  for  pressure  at  the  outlet  of  the  pump 

Wg  uncertainty  for  the  heat  capacity  of  the  evaporator 

Wj  £D  ideal  work  output  in  an  isentropic  expansion  process 

ws,h£Da  uncertainty  for  the  enthalpy  at  the  outlet  of  the  ED-WOW  for  the 

isentropic  expansion  process 

wifij  uncertainty  for  the  entropy  at  the  inlet  of  the  ED-WOW 

wTed  uncertainty  for  the  temperature  at  the  inlet  of  the  ED-WOW 

uncertainty  for  the  volumetric  water  flow  rate 
Vi?#  ed  uncertainty  for  the  ideal  work  produced  through  the  ED-WOW 
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Main  variables  (cont.) 

wz  uncertainty  for  variable  z, 

Wr  uncertainty  for  the  torque  of  the  compressor 

w  uncertainty  for  the  rotational  speed  of  the  compressor 

*  quality  of  C02  or  x  direction  in  rectangular  coordinate 

y  y  direction  in  rectangular  coordinate 

z.  independent  variable  which  affects  the  dependent  variable  A 


Greek 

a 

void  fraction 

(X\ 

upstream  void  fraction 

0C2 

downstream  void  fraction 

Ap 

pressure  difference  across  the  hydraulic  pump 

Apbend 

pressure  loss  due  to  bending 

Apf 

frictional  pressure  loss  across  a  sudden  contraction 

Appump 

pressure  difference  across  the  hydraulic  pump 

Ap  total 

total  pressure  loss 

APl^2 

pressure  loss  due  to  sudden  contraction 

Ap2^ 

pressure  loss  due  to  sudden  enlargement  from  position  2  to  position 

Ap-i^ 

pressure  loss  due  to  sudden  enlargement  from  position  3  to  position 

Greek  (cont.) 


AP4-». 5 
At 

At  intake 

A  U 
AV 
Ax 

£ 

£ed 

£ihx 

Vcomp 

Ved-wow 

Vpump 

71 

e 

p 

p\ 

pi 

Pi 

A 

p 


pressure  loss  at  exit 
time  step 

intake  valve  opening  time 
change  of  internal  energy 
change  of  chamber  volume 
mesh  size 

convergence  tolerance 

isentropic  efficiency  of  the  ED-WOW 

effectiveness  of  the  internal  heat  exchanger 

isentropic  efficiency  of  the  compressor 

isentropic  efficiency  of  the  work  producing  expansion  device 

efficiency  of  the  pump 

constant 

radial  direction  in  cylindrical  coordinate 
density 

density  of  liquid  phase 

density  at  position  2  or  density  after  piston  movement 
density  at  position  3 
density  at  position  4 
viscosity 


Greek  (cont.) 


Vc 

P 

P  tot 

a 


Oe 


T 

CO 

COed 


coefficient  of  friction 
power  loss 
total  power  loss 

area  ratio  or  composite  roughness 
material  limiting  stress 

shear  stress  acting  on  piston  ring  or  torque  of  the  compressor 

revolutionary  speed 
revolutionary  speed  of  ED- WOW 


Executive  Summary 


Develop  an  advanced  lightweight  highly  efficient  Environmental  Control  Unit  (ECU)  by 
utilizing  CO2  cycle  technology  for  AEF  applications  that  reduces  airlift  and  logistics 
requirements.  The  technical  challenges  in  the  development  of  this  technology  include 
development  of  compact  evaporator  and  gas  cooler  based  on  micro-channel  technology, 
development  of  an  energy  recovery  expansion  device,  and  integration  of  the  expansion 
device  and  compressor  into  a  hermetically  sealed  package. 

This  report  include  the  design,  construction,  and  testing  of  a  piston-cylinder  type 
expansion  device  with  output  work  (ED-WOW)  to  be  used  in  a  trans-critical  carbon 
dioxide  (CO2)  air  conditioning  system  that  is  intended  to  replace  R-22-based  or  R-134a- 
based  Field  Deployable  Environmental  Control  Units  (FDECUs).  The  ED-WOW  is  based 
on  a  highly  modified  small  four-cycle,  two-piston  engine  with  a  displacement  of  2  x 
13.26  cm3  that  is  commercially  available.  The  ED-WOW  replaced  the  expansion  valve  in 
an  experimental  trans-critical  CO2  cycle  and  increased  the  system  performance  by  up  to 
10%.  In  addition,  this  document  reports  on  the  development  of  simulation  models  for  the 
ED-WOW  as  well  as  for  the  trans-critical  CO2  cycle.  The  ED-WOW  model  was  used  to 
verify  the  proper  operation  of  the  device  in  the  prototype  application  and  to  identify  the 
major  losses  of  the  device.  The  model  can  be  used  to  re-design  the  next  generation  ED- 
WOW.  The  cycle  model  was  used  to  predict  the  performance  of  several  modifications  to 
the  trans-critical  CO2  cycle  and  identify  which  modification  is  most  suitable  for 
implementation  in  a  carbon  dioxide-based  FDECU. 


ABSTRACT 


Baek,  Joo  Seok.  Ph.D.,  Purdue  University,  August  2002.  Development  of  a  Work 
Producing  Expansion  Device  for  a  Transcritical  Carbon  Dioxide  Cycle.  Major 
Professors:  Dr.  Eckhard  A.  Groll  and  Dr.  Patrick  B.  Lawless,  School  of  Mechanical 
Engineering. 

Because  of  its  high  volumetric  heat  capacity,  CO2  offers  the  potential  for  reduced 
weight  and  volume  in  packaged  systems,  which  is  a  major  focus  for  automobile  and 
military  applications.  The  proceeding  literature  review,  however,  has  shown  that  the 
overall  efficiency  of  the  transcritical  CO2  cycle  decreases  as  the  weight  and  volume  of  the 
system  is  decreased.  System  improvements  must  be  achieved  to  meet  the  goals  of  the 
reduction  of  weight  and  volume,  while  still  retaining  competitive  efficiency. 

A  thermodynamic  analysis  was  performed  for  several  modified  transcritical  CO2 
cycles,  with  a  goal  of  determining  the  proper  configuration  to  achieve  the  optimum 
system  performance.  The  analysis  shows  that  the  system  efficiency  of  the  transcritical 
CO2  cycle  can  be  dramatically  increased  by  modifications  of  the  system  configuration, 
and  that  the  optimum  performance  is  a  strong  function  of  the  operating  conditions. 

The  system  efficiency  of  the  transcritical  CO2  cycle  is  increased  considerably  by 
nearly  100  %  when  the  low-side  pressure  increases  from  2.5  MPa  to  4.5  MPa.  There  is  a 
specific  high-side  pressure  at  which  the  maximum  COP  occurs  in  the  CO2  cycle.  There 
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are  also  specific  pressure  ratios  across  the  compressors  for  the  minimum  compression 
work  and  the  maximum  COP  for  the  two-compressor  system. 

The  employment  of  a  work  extraction  device  in  the  role  of  an  expander  is 
expected  to  give  the  largest  improvement  of  the  system  efficiency.  This  device  will  not 
only  increase  the  heat  removal  capacity  of  the  evaporator  but  also  can  reduce  the 
compression  work.  As  part  of  this  research,  an  expansion  device  with  output  work  (ED- 
WOW)  was  designed,  constructed  and  installed  in  a  prototype  transcritical  C02 
refrigeration  system,  replacing  the  simple  throttling  valve.  During  cycle  tests,  the  use  of 
the  device  provided  between  7%  and  10%  improvement  in  COP.  These  figures  are 
considered  to  be  highly  conservative,  as  the  work  extraction  measurement  was  biased 
toward  a  lower  value  due  to  uncertainty  in  the  load  device  performance. 
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1  INTRODUCTION 


1.1  Background 

Many  nations  of  the  world  have  agreed  on  regulations  prohibiting  the  use  of 
chlorofluorocarbons  (CFCs)  and  hydroclorofluorocarbons  (HCFCs)  as  refrigerants  due  to 
their  potential  to  destroy  the  ozone  layer  of  the  earth’s  atmosphere.  Hydrofluorocarbons 
(HFCs)  are  considered  as  the  replacements  for  CFC  and  HCFC  refrigerants  since  they 
have  zero  ozone  depletion  potential.  However,  HFCs  still  have  a  significant  global 
warming  potential.  Therefore,  substances  with  no  ozone  depletion  potential  and 
negligible  global  warming  potential,  e.g.,  natural  fluids,  are  being  considered  as 
alternative  refrigerants.  Carbon  dioxide,  ammonia,  hydrocarbons,  water,  and  air  are 
gaining  more  and  more  attention.  Some  of  these  refrigerants,  i.e.,  hydrocarbons  and 
ammonia,  show  unfavorable  safety  behaviors  due  to  their  toxicity  and/or  flammability. 
Carbon  dioxide  (CO2)  is  nontoxic  and  non-flammable  and  has  zero  ozone  depletion 
potential  and  negligible  global  warming  potential  as  a  refrigerant.  It  thus  is  close  to  being 
the  ideal  working  medium  for  applications  like  air  conditioning  and  heat  pump  systems. 

Carbon  dioxide  is  also  receiving  strong  consideration  for  automobile  and  military 
applications  due  to  its  high  volumetric  heat  capacity  and  the  elimination  of  recovery  and 
recycling  equipment  and  procedures.  Because  of  its  high  volumetric  heat  capacity,  C02 
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offers  the  potential  for  reduced  weight  and  volume  in  packaged  systems,  which  is  a  major 
focus  for  transportable  field  units  in  the  military.  For  example,  the  U.S.  Army  is  seeking 
an  air  conditioning  system  with  significant  weight  and  size  reduction  and  higher 
efficiency,  resulting  in  reducing  tactical  power  generator  sizes  [Patil,  1998]. 

In  order  to  reduce  the  weight  and  volume  of  the  system,  the  weights  and  volumes 
of  the  evaporator  and  gas  cooler  need  to  decrease.  However,  the  overall  system 
efficiency  will  also  decrease  if  the  heat  exchanger  sizes  are  reduced  due  to  the  increase  in 
approach  temperatures  between  the  refrigerant  and  the  heat  sink  and  source  fluids. 
Robinson  and  Groll  compared  the  ‘Coefficient  Of  Performances’  (COPs)  of  a  packaged 
air  conditioner  using  HCFC-22  as  the  refrigerant  with  those  of  a  carbon  dioxide  based 
packaged  air  conditioner  by  using  a  computer  model  [Robinson  and  Groll,  2000].  Their 
analysis  showed  that  for  a  packaged  air  conditioner  application  using  the  same  evaporator 
size  and  capacity,  the  transcritical  CO2  air  conditioner  will  operate  with  a  COP  which  is 
similar  to  the  simulated  (not  actual)  performance  of  an  U.S.  Army  packaged  air 
conditioner  which  uses  HCFC-22.  As  previously  mentioned,  however,  if  the  size  of  the 
evaporator  is  reduced,  the  heat  removal  capacity  of  the  evaporator  decreases,  which 
would  result  in  a  reduction  in  the  thermodynamic  system  performance. 

Robinson  conducted  a  parametric  study  which  compared  the  actual  performance 
of  an  HFC-134a  based  automobile  air  conditioner  and  the  performance  of  a  carbon 
dioxide  based  automobile  air  conditioner  [Robinson,  2000].  His  analysis  showed  that  for 
an  automobile  application,  the  prototype  transcritical  C02  cycle  system  introduced  by 
Pettersen  [Pettersen,  1995]  will  operate  with  a  COP  which  is  66-75%  (average  70%)  of 
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the  experimentally  determined  COP  of  a  production  HFC-134a  based  automobile  air 
conditioner  while  providing  the  same  evaporator  capacity. 

As  discussed  above,  the  performance  of  a  carbon  dioxide  based  air  conditioning 
system  may  be  less  than  that  of  an  HFC  or  HCFC  based  air  conditioning  system, 
especially  if  the  size  of  the  system  is  reduced.  System  improvements,  therefore,  need  to 
be  realized  to  meet  the  goal  of  the  reduction  of  weight  and  volume  while  still  maintaining 
the  same  or  higher  system  efficiency.  There  are  several  approaches  to  increase  the 
system  efficiency.  These  include  an  increase  of  heat  exchanger  frontal  area,  the  addition 
of  a  liquid-to-suction  line  heat  exchanger,  the  use  of  two-stage  compression  with  an 
intercooler  or  an  economizer,  and  the  use  of  an  expansion  device  that  provides  for  work 
extraction.  The  concepts  behind  these  improvements  are  explained  in  detail  in  the 
following  sections. 


1.2  Basic  Transcritical  CO?  Cycle 

Figure  1.1  shows  a  schematic  of  the  basic  transcritical  CO2  cycle  and  Figure  1.2 
illustrates  the  cycle  on  a  T-s  diagram.  As  shown  in  Figure  1.1,  the  basic  transcritical  CO2 
cycle  consists  of  a  compressor,  a  gas  cooler,  an  expansion  valve,  and  an  evaporator.  The 
cycle  is  composed  of  four  basic  processes;  compression  (1-2),  heat  rejection  (2-3), 
expansion  (3-4)  and  heat  absorption  (4-1)  as  shown  in  Figure  1.2.  In  the  expansion 
process,  the  paths  3-4s,  3-4h  and  3-4w  represent  isentropic  expansion,  isenthalpic 
expansion,  and  expansion  through  an  expansion  device  with  output  work  (ED-WOW), 
respectively.  In  the  compression  process,  the  paths  l-2s  and  1-2  stand  for  the  isentropic 
and  actual  processes,  respectively. 
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The  process  that  differs  most  from  current  conventional  refrigeration  systems  is 
the  supercritical  heat  rejection  process.  The  critical  temperature  and  pressure  of  C02  are 
30.82°C  and  7.353MPa  [ASHRAE,  1997].  During  hot  weather  periods,  the  temperature 
outside  may  be  higher  than  the  critical  temperature  of  C02.  Therefore,  since  the  C02 
temperature  in  the  gas  cooler  must  be  higher  than  the  outside  temperature  to  reject  heat  to 
atmosphere,  the  C02  temperature  needs  to  be  much  higher  than  the  critical  temperature 
(i.e.  in  a  supercritical  state).  Hence,  the  heat  will  be  rejected  from  C02  during  a 
supercritical  process  while  the  heat  is  rejected  from  the  conventional  fluorocarbon-based 
refrigerant  during  a  constant  temperature  and  pressure  two-phase  process.  During  the 
supercritical  heat  rejection  process  the  thermodynamic  properties  of  C02  will  change 
from  a  ‘vapor-like’  to  a  ‘liquid-like’  fluid. 

1 .3  Cycle  Improvements 

In  order  to  increase  the  efficiency  (or  performance)  of  the  basic  cycle  shown  in 
Figure  1.1,  some  processes  need  to  be  modified  or  additional  processes  need  to  be  added 
to  the  basic  cycle.  A  literature  review  was  conducted  on  studies  of  the  transcritical  C02 
cycle  in  cases  where  the  expansion  valve  was  replaced  by  an  expansion  device  producing 
work,  and  in  cases  where  the  transcritical  C02  cycle  was  modified  to  achieve  better 
system  performance. 
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1.3.1  Transcritical  CO2  Cycles  with  Expansion  Device  with  Work  Extraction 

Little  research  has  focused  on  the  topic  of  using  work-producing  expansion 
devices  instead  of  the  conventional  expansion  valve.  Those  studies  available  are 
reviewed  and  summarized  in  this  section. 

Robinson  and  Groll  have  suggested  the  possibility  of  using  a  two-phase 
expansion  turbine  in  place  of  the  expansion  valve  [Robinson  and  Groll,  1998],  The 
authors  developed  a  thermodynamic  computer  model  to  analyze  the  transcritical  CO2 
cycle  with  and  without  an  expansion  work  recovery  turbine.  They  discovered  from  their 
analysis  that  replacing  the  expansion  valve  with  an  expansion  work  recovery  turbine 
reduced  the  process’s  contribution  to  total  cycle  irreversibility  by  35%  and  the  COP  was 
on  average  of  25%  higher  than  the  COP  of  the  transcritical  CO2  cycle  with  an  expansion 
valve  and  maximum  internal  heat  exchanger. 

Li  et  al.  introduced  a  vortex  tube  expansion  device  and  expansion  device 
producing  work  [Li  et  al.,  2000].  Carbon  dioxide  leaving  a  gas  cooler  enters  a  vortex 
tube.  The  vortex  tube  separates  the  carbon  dioxide  into  saturated  liquid  and  superheated 
vapor  at  evaporation  pressure  due  to  the  Ranque-Hilsch  effect.  The  superheated  vapor 
enters  an  auxiliary  heat  exchanger,  dissipates  heat  to  the  environment  and  then  mixes 
with  the  saturated  liquid.  The  quality  of  the  CO2  becomes  less  than  that  of  isenthalpic 
expansion  so  that  the  capacity  of  the  evaporator  increases.  The  authors  also  suggested  a 
piston-cylinder  type  expansion  device  instead  of  the  expansion  valve  that  recovers  some 
energy  while  CO2  expands  through  the  device.  The  energy  can  be  used  to  recover 
mechanical  shaft  work  that  either  directly  drives  a  compressor  or  that  could  be  utilized  to 
generate  electric  power. 
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Figure  1.3  shows  the  schematic  of  the  transcritical  C02  cycle  with  a  work- 
producing  expansion  device  substituted  for  the  expansion  valve.  For  simplicity,  the  work 
producing  device  will  be  referred  to  as  the  ED-WOW  (Expansion  Device  With  Output 
Work).  If  C02  is  expanded  through  the  ED-WOW  (process  3-4w),  the  enthalpy  at  the 
inlet  to  the  evaporator  will  be  lower  than  the  case  of  isenthalpic  expansion  through  an 
ordinary  expansion  valve  (process  3-4h),  as  shown  in  Figure  1.2,  so  that  the  capacity  of 
the  evaporator  will  increase.  When  energy  is  extracted  from  the  C02  through  the  ED- 
WOW  during  the  expansion  process  and  used  to  assist  the  compression  process,  the 
compression  work  will  be  reduced  and  the  COP  of  the  system  will  increase  significantly. 

Li  et  al.  developed  a  computer  model  and  performed  thermodynamic  calculations 
of  the  transcritical  C02  cycle  using  three  expansion  options  [Li  et  al.,  2000].  The  authors 
discussed  that  the  COPs  of  the  cycles  with  a  vortex  tube  and  an  ED-WOW  are  increased 
by  37%  (with  the  expansion  work  output  device  at  an  isentropic  efficiency  of  0.9  and 
with  the  vortex  tube  expansion  process  at  an  efficiency  of  0.5)  compared  to  that  of  the 
cycle  with  an  expansion  valve.  It  was  concluded  that  the  vortex  tube  and  piston-cylinder 
type  expansion  device  have  good  potential  for  cycle  improvement. 

Maurer  and  Zinn  experimentally  investigated  axial,  swash-plate  piston  expanders 
and  internal  gear  expanders  in  a  specially  constructed  liquid-loop  test  stand  [Maurer  and 
Zinn,  1999].  Both  expanders  were  first  tested  as  off-the-shelf  components.  Based  on  the 
experiences  gathered  during  the  first  round  of  testing  the  design  of  both  expanders  was 
modified.  The  suction  port  of  the  axial,  swash-plate  piston  expander  was  redesigned  and 
better  piston  rings  were  introduced.  With  respect  to  internal  gear  expander,  the  tolerance 
of  the  sealing  gap  between  the  gears  and  the  housing  was  reduced.  After  these 
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modifications  were  implemented,  the  isentropic  and  volumetric  efficiencies  that  were 
achieved  with  both  expanders  were  up  to  50%. 

Heyl  et  al.  proposed  an  expander-compressor  unit  to  improve  the  thermodynamic 
efficiency  of  the  transcritical  CO2  cycle  [Heyl  et  al.,  1998].  Figure  1.4  shows  the 
expander-compressor  unit.  As  shown  in  the  figure,  this  unit  is  constructed  as  a  free 
piston  machine  with  two  double-acting  pistons,  which  are  connected  by  a  piston  rod. 
During  the  piston  movement  from  left  to  right  in  cylinder  1,  the  refrigerant  is  expanded 
from  the  gas  cooler  pressure  to  the  evaporator  pressure.  In  cylinder  2,  the  refrigerant  is 
first  compressed  then  discharged.  In  cylinder  3,  the  remaining  refrigerant  is  expanded 
and  additional  refrigerant  is  taken  in  for  the  next  compression  process.  In  cylinder  4  the 
expanded  refrigerant  is  swept  out.  The  authors  mentioned,  however,  that  the  piston 
would  stop  due  to  the  balance  of  the  four  forces  in  the  cylinders  of  Figure  1.4.  For  the 
continuous  operation  of  the  piston,  the  authors  proposed  to  terminate  the  expansion  by 
opening  the  exhaust  valve  in  a  position  of  the  stroke  where  the  resulting  balance  of  forces 
stops  the  movement  of  the  piston.  The  expander  operates  according  to  the  ‘full  pressure 
principle’.  The  authors  mentioned  that  the  resulting  balance  of  forces  of  the  expander- 
compressor  unit  operating  with  the  ‘full  pressure  principle’  was  positive  over  the  entire 
piston  stroke  and  they  ensured  the  functioning  of  the  unit.  In  this  case,  however,  the 
expansion  remains  incomplete.  There  is  unavoidable  loss  in  the  operating  method 
proposed  instead  of  a  complete  expansion.  The  authors  discussed  that  the  calculated 
efficiency  of  the  unit  was  78%  in  relation  to  a  complete  expansion  due  to  the  loss. 
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1.3.2  Transcritical  C02  Cycle  with  Other  Options 

Instead  of  replacing  the  expansion  valve  with  a  work  producing  expansion  device, 
many  other  authors  suggested  other  modifications  of  the  transcritical  C02  cycle  to 
improve  the  system  performance. 

Lorentzen  and  Pattersen  introduced  the  idea  of  an  internal  heat  exchanger  for  a 
transcritical  C02  cycle  [Lorentzen  and  Pattersen,  1993].  Figures  1.5  and  1.6  show  a 
schematic  of  the  transcritical  C02  cycle  with  an  internal  heat  exchanger  and  the  cycle  on 
a  T-s  diagram,  respectively.  The  C02  exiting  the  gas  cooler  rejects  heat  to  the  C02 
exiting  the  evaporator  through  the  internal  heat  exchanger  (also  known  as  a  ‘liquid-to- 
suction  line’  heat  exchanger).  In  the  figures  the  processes  3’-3  and  1-1’  correspond  to  the 
internal  heat  exchange.  There  are  two  major  functions  of  the  internal  heat  exchanger 
[Robinson,  2000]:  1)  It  ensures  that  the  fluid  entering  the  compressor  is  devoid  of  any 
entrained  liquid  C02;  2)  It  reduces  the  temperature  of  the  fluid  prior  to  expansion  that 
results  in  lower  vapor  quality  of  the  C02  entering  the  evaporator.  This  in  turn  increases 
the  capacity  of  the  evaporator. 

G.T.  Voorhees  proposed  a  transcritical  C02  cycle  with  an  economizer  in  1905 
that  is  also  known  as  the  ‘Multiple-Effect  cycle’  [Voorhees,  1979].  Figure  1.7  illustrates 
this  cycle.  The  expansion  process  occurs  through  two  expansion  valves  with  an 
accumulator  that  separates  vapor  phase  C02  from  liquid  phase  C02.  The  vapor  phase 
C02  mixes  with  the  superheated  C02  exiting  the  lst-stage  compressor  and  the  mixture 
enters  the  2nd-stage  compressor.  The  liquid  phase  C02  expands  again  and  enters  the 
evaporator  as  shown  in  Figures  1.7  and  1.8.  Figure  1.8  shows  the  system  on  T-s  diagram. 
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As  shown  in  Figure  1.8,  the  enthalpy  of  CO2  entering  the  evaporator  is  reduced  and  it 
results  in  increased  capacity  of  the  evaporator. 

The  compressor  work  input  can  also  be  reduced  by  means  of  multistage 
compression  with  intercooling.  Figure  1.9  illustrates  a  two-stage  compressor  with  an 
intercooler  [Moran  and  Shapiro,  1992].  The  accompanying  p-v  and  T-s  diagrams  show 
the  states  for  internally  reversible  processes.  Process  1-c  denotes  an  isentropic 
compression  from  state  1  to  state  c  where  the  pressure  is  />,.  In  process  c-d  the  gas  is 
cooled  at  constant  pressure  from  temperature  Tc  to  Tj.  Process  d-2  is  an  isentropic 
compression  from  state  d  to  state  2.  The  work  input  per  unit  mass  flow  is  represented  on 
the  p-v  diagram  by  area  l-c-d-2-a-b-l.  Without  intercooling  the  gas  would  be 
compressed  isentropically  in  a  single  stage  from  state  1  to  state  2’  and  the  work  would  be 
represented  by  area  l-2’-a-b-l.  The  crosshatched  area  on  the  p-v  diagram  represents  the 
reduction  in  work  that  would  be  achieved  with  intercooling.  The  size  of  the  crosshatched 
area  on  the  p-v  diagram  of  Figure  1.9  representing  the  reduction  of  compression  work 
depends  on  both  the  temperature  Tj  at  the  exit  of  the  intercooler  and  the  intercooler 
pressure  pt.  By  properly  selecting  7>  and  p^  the  total  work  input  to  the  compressor  can  be 
minimized.  This  concept  can  be  applied  to  the  transcritical  CO2  cycle. 

Figure  1.10  shows  a  schematic  of  the  transcritical  CO2  cycle  with  two-stage 
compression  with  intercooling.  Figure  1.11  illustrates  the  T-s  diagram  of  the  system. 
The  compression  work  will  be  reduced  with  this  cycle  as  discussed  above  but  the 
capacity  of  the  evaporator  does  not  change.  Therefore,  the  coefficient  of  performance 
(COP)  of  this  cycle  will  increase  compared  to  the  basic  transcritical  CO2  cycle. 


10 


Heyl  et  al.  analyzed  several  modifications  of  the  basic  transcritical  C02  cycle 
[Heyl  et  al.,  1998].  The  author  considered  single-stage  and  two-stage  compressions.  The 
single-stage  compression  cycle  included  the  basic  transcritical  C02  cycle  (Figure  1.1),  the 
basic  cycle  with  an  internal  heat  exchanger  (Figure  1.5),  the  transcritical  C02  cycle  with 
an  expander  (Figure  1.3),  and  the  transcritical  C02  cycle  with  an  expander  and  an  internal 
heat  exchanger.  The  two-stage  compression  cycle  included  the  transcritical  C02  cycle 
with  two-stage  compressions  and  intercooling  (Figure  1.10),  the  transcritical  C02  cycle 
with  an  economizer  (Figure  1.7),  and  the  transcritical  C02  cycle  with  an  expander,  two 
compressors,  and  intercooling.  Their  analysis  showed  that  the  COP  increased  with  the 
modification  of  the  basic  cycle  and  the  COP  became  the  maximum  when  the  work  output 
extracted  during  the  C02  expansion  was  used  to  run  the  2nd-stage  compressor.  The  COP 
improvement  was  about  65%  compared  the  basic  transcritical  C02  cycle. 

1.4  Research  Objectives 

Because  of  its  high  volumetric  heat  capacity,  C02  offers  the  potential  for  reduced 
weight  and  volume  in  packaged  systems,  which  is  a  major  focus  for  automobile  and 
military  applications.  The  proceeding  literature  review,  however,  has  shown  that  the 
overall  efficiency  of  the  transcritical  C02  cycle  decreases  as  the  weight  and  volume  of  the 
system  decrease.  System  improvements  must  be  achieved  to  meet  the  goals  of  the 
reduction  of  weight  and  volume,  while  still  retaining  competitive  efficiency. 

As  previously  discussed,  the  employment  of  a  work  extraction  device  in  the  role 
of  an  expander  is  expected  to  give  the  largest  improvement  of  the  system  efficiency.  This 
device  will  not  only  increase  the  heat  removal  capacity  of  the  evaporator  but  also  can 
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reduce  the  compression  work.  Most  of  the  research  above  suggested  the  conceptual 
design  of  the  expansion  device.  However,  only  Heyl  tested  on  expander  in  a  transcritical 
CO2  cycle  which  had  several  limitations  as  outlined  above  [Heyl,  2000].  A  major 
objective  of  this  research  was  to  design,  construct,  test,  and  analyze  a  work  extraction 
expansion  device.  This  device  was  installed  in  a  prototype  transcritical  CO2  refrigeration 
system,  replacing  the  simple  throttling  valve.  Experiments  were  performed  to 
characterize  the  effect  of  the  expansion  device  on  the  cycle  performance.  This  data  was 
compared  to  analytical  predictions  of  the  cycle  performance. 

Another  objective  of  this  research  is  to  perform  thermodynamic  analyses  for  the 
transcritical  CO2  cycles  using  an  economizer  (Economizer  cycle)  and  an  intercooler 
(Intercooler  cycle),  which  are  described  above,  with  the  goal  to  determine  the 
configuration  that  achieves  optimum  system  performance.  Especially,  intermediate 
pressures  were  investigated  for  two  compressor  systems  which  gave  the  maximum 
system  efficiency.  Such  an  analysis  of  the  optimum  intermediate  pressure  of  the 
transcritical  CO2  cycle  with  intercooling  has  not  been  published  to  date. 


1.5  Organization  of  This  Document 

In  Chapter  1,  the  background  of  this  research  is  mentioned  and  the  literature 
review  has  been  calculated  on  the  topic  of  the  improvement  of  the  thermodynamic  system 
performance  of  the  transcritical  CO2  cycles.  The  objectives  of  this  research  are  discussed 
in  Chapter  1  as  well.  In  Chapter  2,  the  technical  approach  to  conduct  this  research  is 
described.  First,  the  methods  of  the  cycle  analyses  are  given  for  transcritical  CO2  cycles 
with  an  economizer,  an  intercooler,  and  an  expansion  device  producing  work.  Then,  the 
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determination  of  a  work  producing  expansion  device  is  discussed.  A  detailed  description 
of  a  theoretical  model  for  the  expansion  device  is  also  given  in  Chapter  2.  In  addition, 
Chapter  2  provides  a  summary  of  the  literature  review  on  the  topics  of  the  minor  flow 
losses,  piston  ring  scuffing,  and  lubrication.  Furthermore,  the  testing  procedure  for  the 
experiments  with  the  expansion  device  in  a  transcritical  CO2  cycle  is  explained  in 
Chapter  2.  The  design  of  the  expansion  device  is  discussed  in  detail  in  Chapter  3.  The 
results  of  cycle  analyses  are  discussed  in  Chapter  4.  The  results  of  the  experiments  are 
also  discussed  in  detail  in  Chapter  4.  The  expansion  device  is  analyzed  using  the 
theoretical  model  for  its  design.  The  experimental  and  simulation  results  are  compared  in 
Chapter  4  as  well.  A  summary  and  conclusions  of  the  research  accomplishments  and 
recommendations  for  design  improvements  of  the  expansion  device  are  given  in  Chapter 
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Evapaorator 


Figure  1.1:  Schematic  of  basic  transcritical  carbon  dioxide  cycle. 


s 

Figure  1.2:  Temperature-entropy  diagram  of  transcritical  carbon  dioxide  cycle. 


Gas  Cooler 


Figure  1.3:  Schematic  of  transcritical  C02  cycle  with  expansion  device  with  output  work 

(ED-WOW)  replacing  expansion  valve. 
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Figure  1.4:  Free  piston  expander-compressor-unit  [Heyl  et  al.,  1988]. 
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Figure  1.5:  Schematic  of  transcritical  C02  cycle  with  liquid-to-suction  line  (internal) 

heat  exchanger. 


Figure  1.6:  Temperature-entropy  diagram  of  transcritical  C02  cycle  with  internal  heat 

exchange. 
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Evapaorator 

Figure  1.7:  Schematic  of  transcritical  C02  cycle  with  an  economizer. 


Figure  1.8:  Temperature-entropy  diagram  of  transcritical  C02  cycle  with  an  economizer 
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Compressor  Compressor 


Figure  1.9:  Two-stage  compression  with  intercooling  [Moran  and  Shapiro,  1992], 


Evapaorator 

Figure  1.10:  Schematic  of  transcritical  CO2  cycle  with  two-stage  compression  with 

intercooling. 
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Figure  1.11:  Temperature-entropy  diagram  of  transcritical  CO2  cycle  with  two-stage 

compression  with  intercooling. 
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2  TECHNICAL  APPROACH 

2.1  Cycle  Analysis 

One  of  the  objectives  of  this  research  was  to  perform  the  thermodynamic  analyses 
to  investigate  the  thermodynamic  behavior  of  the  transcritical  CO2  cycle  and  to  find  the 
optimum  operating  conditions  for  the  maximum  system  performance.  A  thermodynamic 
model  was  developed  to  perform  theses  analyses. 

The  analyses  have  been  performed  in  two  parts.  In  the  first  part,  the  cycles 
without  the  expansion  device  producing  work  have  been  considered.  A  transcritical  CO2 
cycle  with  an  economizer  (Economizer  cycle)  or  an  intercooler  (Intercooler  cycle)  has 
been  analyzed.  Since  pressure  is  one  of  the  most  important  parameters  in  the  air 
conditioning  system,  the  pressure  effects  are  examined  in  detail.  Especially,  the  effects 
of  the  pressure  ratios  across  the  compressors,  the  low-side  pressure,  and  the  high-side 
pressure  on  the  system  performance  of  the  transcritical  CO2  cycle  are  investigated. 

In  the  second  part,  the  effect  of  an  expansion  device  with  output  of  work  (ED- 
WOW)  on  the  system  performance  of  the  transcritical  CO2  cycles  has  been  analyzed. 
There  are  two  advantages  in  using  the  ED-WOW  instead  of  a  throttling  expansion  valve 
on  the  system  performance.  The  first  advantage  is  that  lower  enthalpy  can  be  obtained  at 
the  inlet  of  the  evaporator.  The  other  one  is  that  work  can  be  extracted  through  the 


20 


expansion  device.  If  this  work  is  used  to  run  the  compressor,  the  input  work  for  the 
compression  of  carbon  dioxide  can  be  reduced.  In  the  analysis,  these  effects  have  been 
considered.  First,  the  use  of  the  work  output  extracted  from  C02  during  the  expansion 
process  has  not  been  considered.  That  is,  only  the  enthalpy  decrease  through  the 
expansion  device  has  been  considered.  Secondly,  the  work  output  through  the  ED-WOW 
has  also  been  considered  with  the  decrease  of  the  enthalpy  and  it  has  been  assumed  that 
the  work  is  used  to  run  the  compressor. 


2.1.1  Assumptions 

The  main  assumptions  for  the  analysis  of  the  thermodynamic  performance  are  as 
follows: 

1.  Flow  is  steady  state. 

2.  Compression  processes  is  isentropic  but  the  expansion  process  is  isenthalpic. 

3.  There  are  no  pressure  losses  in  the  heat  exchangers  and  piping. 

4.  The  atmosphere  is  the  heat  sink.  The  C02  temperature  equilibrates  with  the 
ambient  air  temperature  due  to  the  heat  transfer  through  the  gas  cooler  and 
intercooler. 

2.1.2  Transcritical  C02  Cycles  with  an  Economizer  and  Intercooler 
When  the  transcritical  C02  cycle  features  a  two-stage  compression  system,  the 
COP  can  be  optimized  by  choosing  different  pressure  ratios  for  each  stage.  The  system 
performances  of  the  Economizer  cycle  and  Intercooler  cycle  were  examined  and 
compared  with  that  of  the  Base  cycle.  The  investigated  parameters  are  the  effect  of  the 
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pressure  ratios  across  the  compressors,  the  effect  of  low-side  pressure,  and  the  effect  of 
high-side  pressure  on  system  performance. 

Two  basic  CO2  cycles  are  analyzed  here;  namely  the  transcritical  CO2  cycle  with 
an  economizer  (Economizer  cycle,  Figure  1.7)  and  the  transcritical  CO2  cycle  with  an 
intercooler  (Intercooler  cycle,  Figure  1.10).  Both  cycles  include  two  compressors.  As 
discussed  in  Chapter  1,  both  cycles  are  modified  from  the  basic  transcritical  CO2  cycle 
(Base  cycle)  of  Figure  1.1.  The  Base  cycle  consists  of  a  compressor,  a  gas  cooler,  an 
evaporator,  and  an  expansion  valve. 

The  Economizer  cycle  is  composed  of  two  compressors,  a  gas  cooler,  an 
evaporator,  two  expansion  valves  and  an  accumulator.  After  the  first  expansion,  the 
saturated  vapor  of  C02  at  the  intermediate  pressure  in  the  accumulator  flows  to  the  inlet 
to  the  2nd-stage  compressor  and  mixes  with  the  CO2  coming  from  the  outlet  of  the  1st- 
stage  compressor.  The  saturated  liquid  C02  expands  again  through  another  expansion 
valve.  Following  this  valve,  the  two-phase  flow  of  carbon  dioxide  enters  the  evaporator. 
The  enthalpy  of  the  carbon  dioxide  at  the  inlet  to  the  evaporator  decreases  as  shown  in 
Figure  1.8  and  it  will  increase  the  capacity  of  evaporator.  This  results  in  an  increase  of 
COP  of  the  C02  cycle. 

Figure  1.10  represents  the  Intercooler  cycle.  As  shown  in  this  figure,  the 
Intercooler  cycle  comprises  two  compressors,  a  gas  cooler,  an  expansion  valve,  an 
evaporator  and  an  intercooler.  The  CO2  coming  from  the  outlet  of  the  lst-stage 
compressor  is  cooled  down  by  removing  heat  through  the  intercooler.  Then  it  flows  to 
the  inlet  to  the  2nd-stage  compressor.  As  discussed  in  Chapter  1,  it  results  in  saving  the 
compression  work  and  increase  of  COP. 
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2. 1 .2. 1  Calculation  Methods 

The  calculation  methods  for  the  thermodynamic  performance  of  the  CO2  cycles 
are  given  below.  The  method  depends  on  the  type  of  transcritical  C02  cycle  being 
analyzed. 


2. 1.2. 1.1  Calculation  Method  for  Base  Cycle 

In  this  analysis,  the  subscript  ‘B’  is  used  on  all  variables  to  represent  the  Base 

cycle  (Figure  1.1). 

The  C02  temperature  and  pressure  at  the  inlet  to  the  compressor  (T[n  and  pB\)  and 
the  air  temperature  ( Tair )  are  given  as  input  variables.  The  temperature  Tb\  is  higher  than 
or  equal  to  the  saturated  temperature  ( Tsat )  at  pB\.  The  entropy  (sB  1)  and  enthalpy  (hB  1)  at 
TBl  and  pB  1  are  found  using  the  C02  thermodynamic  property  data  by  Span  and  Wagner 
[Span  and  Wagner,  1996].  In  the  Base  cycle,  the  high-side  pressure  (pB2 )  is  the  same  as 
the  pressure  (pn  or  Pi2>  see  the  following  sections)  at  the  outlet  of  the  2  -stage 
compressor  of  the  Economizer  cycle  or  Intercooler  cycle.  For  the  Base  cycle 

Pbi  =  Pb\  x  (pressure  ratio)  .  (2-1) 

The  temperature  after  compression  (Tbt)  is  found  for  the  same  sB\  and  pB 2  using 
the  thermodynamic  property  data  [Span  and  Wagner,  1996].  Then  the  enthalpy  (hB2)  at 
the  outlet  of  the  compressor  is  found  at  TB2  and  pB2.  The  compression  work  per  unit  mass 
flow  is  calculated  by  subtracting  the  enthalpy  of  C02  at  the  inlet  to  the  compressor  from 
the  enthalpy  at  the  outlet  of  the  compressor  as  follows: 


Wc  B  ~  h-Bl  ^B\  • 


(2-2) 
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After  heat  rejection  to  the  ambient  atmosphere  through  the  gas  cooler,  the  CO2 
temperature  (7^3)  becomes  Ta\r.  Then  the  enthalpy  (Jibs)  at  the  outlet  of  the  gas  cooler  is 
found  at  Tb3  and  pB3  =  pB2 .  The  enthalpy  after  expansion,  hE 4,  at  pE\  becomes  the  same 

as  Hb3  based  on  the  assumption  of  isenthalpic  expansion  process.  Then,  the  heat  capacity 
of  the  evaporator  is  calculated  as  follows: 

Qe,B  =  hB\  ~  kB4  ■  O3) 

The  COP  of  the  Base  cycle  then  becomes 

COPB=-^-.  (2-4) 

W 

vvc,B 


2. 1 .2. 1 .2  Calculation  Method  for  Economizer  Cycle 

The  subscript  l£”  is  employed  to  represent  the  Economizer  cycle  in  this  section 
(Figure  1.7). 

The  solution  proceeds  by  first  specifying  the  CO2  temperature  and  pressure  at  the 
inlet  to  the  lst-stage  compressor  (TE]  and  pe\)  and  the  ambient  air  temperature  (Tair).  The 
entropy  (sa)  and  enthalpy  (hE\)  are  found  at  7a  and  pE  1  using  the  CO2  thermodynamic 
property  data  [Span  and  Wagner,  1996].  Identical  or  different  (based  on  analysis) 
pressure  ratios  across  the  lst-stage  and  2nd-stage  compressors  are  assumed,  and  then  the 
pressures  at  the  outlets  of  the  compressors,  pEy  and  pE 2,  are  calculated.  The  temperature 
at  the  outlet  of  the  lst-stage  compressor,  TEy,  is  found  based  on  the  assumption  of 
isentropic  compression  (^£r  =  s£1)  and  the  known  pEy.  Then  the  enthalpy  at  the  outlet  of 
the  lst-stage  compressor  ( hEy )  is  found  at  TEy  and pEy. 
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The  enthalpy  at  the  outlet  of  the  gas  cooler  (hE 3)  is  found  at  Tajr  and  pE3  —  pE2  ■ 
Then  the  saturated  vapor  and  liquid  enthalpies  (hE5  and  hE y)  are  found  at 
Pei'  ~  Pe r  =  Pes  •  The  quality  of  C02  in  the  accumulator  becomes 

x  -  ^£3  ~  ^£3'  (2-5) 

^£5  _  ^£3' 

The  enthalpy  at  the  outlet  of  the  2nd  expansion  valve  at  pE  1,  hE 4,  becomes  the  same 
as  ha”  based  on  the  assumption  of  isenthalpic  expansion  process.  The  heat  capacity  of 
the  evaporator  is  calculated  as  follows: 

Qe,E  =(hEi-hEi)x(l-x).  (2-6) 

The  compression  work  of  the  lst-satge  compressor  is  calculated  as  follows: 

Wcl'E={hEV-hE])x(l-x).  (2-7) 

The  enthalpy  at  the  inlet  to  the  2nd-stage  compressor  (hEr)  after  the  mixing  of  the 
saturated  vapor  C02  flowed  from  the  accumulator  and  the  superheated  C02  at  the  outlet 
of  the  lsl-stage  compressor  is  calculated  at  the  intermediate  pressure  {pEy  —  pEl- )  as 
follows: 

hEr  =hE5xx  +  hEV  x (1  - x) .  (2-8) 

Then,  the  temperature  of  the  C02  mixture  (TEr)  is  found  with  hEy<  and  pEy  using  the 
thermodynamic  property  data  [Span  and  Wagner,  1996].  The  entropy  of  the  C02  mixture 
(sEy)  is  found  with  TEy  and  pEy.  The  temperature  of  C02  at  the  outlet  of  the  2nd-stage 
compressor  (Tei)  is  found  with  the  assumption  of  isentropic  compression  process 
(5£2  =  sey)  an<^  high-  side  pressure,  pez .  Then  the  enthalpy  at  the  outlet  of  the 
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2nd-stage  compressor  (He 2)  is  found  at  Tej  and  pEi-  The  compression  work  of  the  2nd- 
stage  compressor  is  calculated  as  follows: 

Wc2,E=hE2-hEl..  (2-9) 

The  COP  of  the  Economizer  cycle  is  then  calculated  as 

COPE  = - ^ - .  (2-10) 

wchE  +  wc2'E 

2. 1.2. 1.3  Calculation  Method  for  Intercooler  Cycle 

The  subscript  7’  is  used  on  all  variables  to  represent  the  Intercooler  cycle  in  this 
section  (Figure  1.10). 

The  solution  proceeds  by  first  specifying  the  CO2  temperature  and  pressure  at  the 
inlet  to  the  lst-stage  compressor  (7/j  and  pn)  and  the  ambient  air  temperature  (Tair).  The 
entropy  (sn)  and  enthalpy  (hn)  are  found  with  the  given  Tn  and  pn  using  the  CO2 
thermodynamic  property  data  [Span  and  Wagner,  1996],  Identical  or  different  (based  on 
analysis)  pressure  ratios  across  the  lst-stage  and  2nd-stage  compressors  are  assumed,  and 
then  the  pressures  at  the  outlets  of  the  compressors,  p/y  and  pn,  are  calculated.  The 
temperature  at  the  outlet  of  the  lst-stage  compressor,  Tjy,  is  found  based  on  the 
assumption  of  isentropic  compression  (slV  =  sn)  and  the  known  pn>.  Then  the  enthalpy 
at  the  outlet  of  the  lst-stage  compressor  (h/y)  is  found  at  Tn-  and  pjy.  The  compression 
work  of  the  lst-stage  compressor  is  calculated  by  subtracting  the  enthalpy  at  the  inlet  to 
the  lst-stage  compressor  from  the  enthalpy  at  the  outlet  of  the  ls,-stage  compressor  as 
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After  heat  transfer  to  the  ambient  atmosphere  through  the  intercooler,  the  C02 
temperature  {Tit)  becomes  Tair  at  the  outlet  of  the  intercooler.  The  entropy  O/2O  and 
enthalpy  (/1/2O  at  the  inlet  to  the  2nd-stage  compressor  are  found  with  Tit  and  p!2-  =  p,v 
using  the  C02  thermodynamic  property  data  [Span  and  Wagner,  1996].  Then,  the 
temperature  at  the  outlet  of  the  2nd-stage  compressor,  7>2,  is  found  with  the  assumption  of 
the  isentropic  compression  process  {sn  =  s/r)  and  the  known  high-side  pressure,  pn. 
The  enthalpy  at  the  outlet  of  the  2nd-stage  compressor,  hI2,  is  found  at  Ta  and  pn  using  the 
CO2  thermodynamic  property  data  [Span  and  Wagner,  1996].  The  compression  work  of 
the  2nd-stage  compressor  is  calculated  by  subtracting  the  enthalpy  at  the  inlet  to  the  2 
stage  compressor  from  the  enthalpy  at  the  outlet  of  the  2nd-stage  compressor  as  follows: 

Wc2J=hI2-hIT.  (2-12) 

After  heat  rejection  to  the  ambient  atmosphere  through  the  gas  cooler,  the  C02 
temperature  (T/3)  becomes  Tair  at  the  outlet  of  the  gas  cooler.  Then  the  enthalpy  at  the 
inlet  to  the  expansion  valve,  hn,  is  found  with  7/3  and  pl3  =  pn  using  the 
thermodynamic  property  data  [Span  and  Wagner,  1996].  After  expansion,  the  enthalpy, 
hlA,  at  pi 4  (=  pn)  becomes  the  same  as  hR  by  assuming  isenthalpic  expansion  process. 
Therefore,  the  heat  removal  capacity  of  the  evaporator  can  be  calculated  as  follows: 

Qe,i=hn~hu •  (2'13) 

The  COP  of  the  Intercooler  cycle  is  then  calculated  as  follows: 
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2.1.3  Transcritical  CO2  Cycles  with  Work  Extraction 
The  transcritical  CO2  cycle  with  an  ED-WOW  has  been  analyzed.  The  main 
purpose  of  this  analysis  is  to  examine  the  effect  of  the  ED-WOW  on  the  system 
performance.  The  system  efficiencies  are  calculated  for  the  cycles  without  the  expansion 
device  and  the  cycles  with  the  ED-WOW.  For  the  cycles  with  the  ED-WOW,  two  cases 
have  been  considered.  First,  the  decrease  of  the  CO2  enthalpy  through  the  ED-WOW  has 
been  considered  for  the  system  performance.  Then,  the  work  output  through  the  device 
as  well  as  the  decrease  of  the  enthalpy  has  been  considered. 


2, 1.3.1  Calculation  Method 

The  calculation  method  for  this  analysis  is  the  same  as  that  for  the  cycles  without 
an  expansion  device  producing  work  in  previous  sections.  The  only  differences  are  that 
the  isentropic  efficiencies  are  considered  for  the  compressors  and  the  ED-WOW  and  the 
effectiveness  of  the  internal  heat  exchanger  is  assumed.  The  assumed  isentropic 
efficiency  of  the  compressor  ( t]comp )  is  0.7.  The  isentropic  efficiency  of  the  ED-WOW 
( Ved-wow )  is  considered  to  be  0.5  and  the  effectiveness  of  the  internal  heat  exchanger 
( £ihx )  is  assumed  to  be  0.7.  These  component  efficiencies  are  defined  as  follows: 
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Where,  the  descriptions  about  the  variables  are  given  in  the  nomenclature  section. 

The  cycles  that  were  analyzed  are  the  Base  cycle  (Case  1),  Base  cycle  with  the 
ED-WOW  (Case  2),  internal  heat  exchanger  cycle  (IHX  cycle.  Case  3),  IHX  cycle  with 
the  ED-WOW  (Case  4),  Intercooler  cycle  (Case  5),  Intercooler  cycle  with  the  ED-WOW 
(Case  6),  and  Economizer  cycle  (Case  7).  Schematics  of  these  cycles  are  shown  in 
Figures  2.1  through  2.7. 


2.2  Expansion  Device 

2.2.1  Overview 

A  literature  and  technical  review  has  been  conducted  to  determine  the  suitability 
of  various  expansion  devices  which  produce  work.  Both  continuous  flow  and  positive 
displacement  devices  were  investigated.  The  operating  conditions  for  a  typical  expansion 
device  in  a  prototype  transcritical  CO2  cycle  are  given  in  Table  2.1.  These  design 
operating  conditions  were  determined  based  on  a  study  of  the  transcritical  CO2  cycle  by 
Li  and  Groll  [Li  and  Groll  2001]. 

Because  of  their  high  inherent  efficiency,  turbomachinery  type  continuous  flow 
expanders  are  considered  a  prime  candidate  CO2  cycle  application.  Several  issues  must 
be  addressed  before  their  suitability  can  be  properly  judged.  First,  erosion  by  liquid 
impingement  may  be  an  issue  due  to  the  two-phase  flow  expected  in  the  expansion 
device.  Second,  the  rotational  speed  required  for  efficient  performance  should  be 
determined  to  evaluate  the  required  bearing  technology. 
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Thiruvengadam  presents  a  theory  for  prediction  of  blade  erosion  based  on  the 
characteristics  of  the  flow  and  the  turbine  material  [Thiruvengadam,  1966].  This  theory 
assumes  erosion  damage  by  repetitive  impact  by  spherical  droplets,  and  expresses  the 
limiting  criteria  as: 


Where: 

ae  -  material  limiting  stress 
px  =  density  of  liquid  phase 
Cj  =  sonic  velocity  in  liquid  phase 
U,  =  impact  velocity 

Alvarez  made  a  preliminary  study  of  the  potential  for  turbo-type  expanders  for  a 
transcritical  C02  cycle  application  and  applied  this  theory  to  an  expander  with  a  nominal 
inlet  C02  pressure  and  temperature  of  10  MPa  and  40°C  [Alvarez,  1999].  With  steel  as 
the  blade  material,  velocities  of  up  to  800  m/s  were  determined  to  be  tolerated  without 
danger  of  erosion.  These  values  are  far  in  excess  of  any  local  velocity  to  be  expected  in  a 
machine  of  this  class  and  type,  and  thus  erosion  due  to  liquid  impingement  was  ruled  out 
as  a  concern  for  this  device. 

Alvarez  also  applied  specific  speed  and  specific  diameter  correlations  [Balje, 
1981]  to  estimate  the  rotational  speed  required  for  a  turbine  operating  between  the  10 
MPa  inlet  pressure  and  0.4  MPa  exit  pressure  at  a  flow  rate  of  0.08  kg/s.  For  this 
application,  a  rotational  speed  on  the  order  of  500,000  rpm  was  required  to  achieve  a  high 
efficiency  design.  Estimations  for  a  partial  admission  turbine  determined  that  speeds  as 
low  as  70,000  rpm  could  be  achieved  with  a  25  mm  diameter  radial  turbine.  However, 
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this  required  an  extremely  low  4%  admission  resulting  in  a  maximum  efficiency  of  50%. 
As  correlations  in  that  range  are  extrapolations  of  data  taken  from  higher  admission 
device,  the  actual  efficiency  may  well  be  lower  even  if  the  machine  is  well  designed  and 
manufactured  to  close  tolerances.  In  addition,  the  Balje  correlations  do  not  include  the 
impact  of  two-phase  flow.  Liquid  in  the  turbine  flow  passages  can  be  expected  to 
migrate  due  to  centrifugal  effects  and  will  thus  have  significant  impact  on  true 
performance. 

Although  continuous  flow  machines  have  potential  for  the  expander  in  large  scale 
CO2  cycles,  for  the  capacity  of  concern  here  the  small  dimensions  and  high  rotational 
speeds  make  them  a  high  risk  design  proposal.  Thus,  a  positive  displacement  device  was 
chosen  as  the  type  of  expander  to  employ  in  the  prototype  system.  Because  of  its 
inherent  simplicity  and  the  vast  knowledge  base  available  on  inherent  losses,  a  piston- 
cylinder  assembly  device  was  chosen  as  the  appropriate  expansion  device  for  the 
prototype  application. 

The  piston  device  consists  of  a  cylinder,  piston  and  connecting  rod  as  shown  in 
Figure  2.8.  The  piston  is  connected  to  a  crank  mechanism  that  drives  the  shaft.  The  shaft 
in  turn  is  connected  to  a  hydraulic  pump  used  as  a  mechanical  loading  device.  The 
expansion  process  is  controlled  using  fast-acting  solenoid  valves  as  intake  and  exhaust 

valves. 

After  determination  of  the  appropriate  type  of  device,  a  representative  example 
was  designed  and  fabricated  to  form  the  test  subject  for  the  first  part  of  this  research 
effort.  Initial  testing  was  performed  on  this  device  in  isolation,  the  shaft  connected  to  a 
torque  meter  and  loading  device.  Measurements  of  mass  flow  rate,  expansion  ratio,  and 
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work  extraction  was  conducted  for  the  purpose  of  characterizing  component  efficiency. 
This  effort  is  presented  in  detail  in  Chapter  3. 


2.2.2  Thermodynamic  Processes  of  the  Piston  Expander 
There  are  three  processes  that  occur  within  the  time  period  of  one  revolution: 
intake,  expansion  and  exhaust.  During  the  intake  process,  the  intake  valve  opens  and  the 
supercritical  high  pressure  CO2  flows  into  the  cylinder  as  the  piston  moves  down  from 
Top  Dead  Center  (TDC)  to  the  location  of  expansion  initiation  (LOE)  as  shown  in  Figure 
2.8.  When  the  piston  reaches  LOE,  the  intake  valve  closes  and  the  expansion  process 
begins  with  the  piston  moving  down  further,  as  the  CO2  expands  from  high  pressure  and 
temperature  to  low  pressure  and  temperature.  When  the  piston  reaches  Bottom  Dead 
Center  (BDC),  the  exhaust  valve  opens.  During  the  second  stroke  the  expanded  C02  is 
swept  out  to  the  inlet  of  the  evaporator.  As  the  piston  reaches  TDC,  the  exhaust  valve 
closes.  The  next  revolution  begins  by  opening  the  intake  valve  again.  Figure  2.9  shows 
the  p-v  diagram  of  the  piston  process  when  the  process  is  ideal. 


2.2.3  Theoretical  Modeling  of  the  ED-WOW  Device 
In  order  to  analyze  the  performance  of  the  prototype  expansion  device,  a 
theoretical  model  was  developed.  This  model  was  used  to  both  verify  proper  operation  of 
the  device  in  the  prototype  application,  and  to  identify  loss  mechanisms  in  the  device  so 
as  to  guide  design  improvements.  This  model  should  feature  sufficient  physical  detail  to 
model  the  unique  characteristics  of  the  carbon  dioxide  application.  No  research  has  been 
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published  on  modeling  the  expansion  device  for  a  transcritical  C02  cycle.  Hence,  this  is 
the  first  model  to  predict  the  thermal  and  mechanical  behavior  of  the  expansion  device. 

2.2.3. 1  Basic  Loss  Mechanisms 

The  minor  losses  other  than  boundary  friction  have  been  reviewed  for  the  mass 
flow  into  and  from  the  cylinder  chamber.  The  topics  for  the  piston  scuffing  and 
lubrication  have  also  been  reviewed  to  predict  the  piston  ring  behaviors. 

2.2.3. 1.1  Minor  Flow  Losses 

In  order  to  simulate  the  C02  flow  into  and  from  the  cylinder  chamber,  minor 
losses  due  to  sudden  contractions,  sudden  enlargements,  and  bends  of  flow  path  must  be 
characterized.  Figure  2.10  shows  the  configuration  of  the  solenoid  valve,  fittings  and  the 
piston-cylinder  section.  The  C02  pressure  loss  from  the  outlet  of  the  solenoid  valve  to 
the  cylinder  chamber  is  assumed  to  consist  of  losses  due  to  sudden  contraction  from  point 

1  to  point  2  in  the  solenoid  valve,  the  bend  of  flow  path  in  the  solenoid  valve  (from  point 

2  to  point  3),  sudden  enlargement  from  solenoid  valve  to  the  fitting  (from  point  3  to  point 
4),  and  sudden  enlargement  from  the  fitting  to  the  cylinder  chamber  (from  point  4  to 
point  5).  The  pressure  loss  due  to  the  wall  friction  on  few  straight  sections  of  the  channel 
is  considered  much  smaller  than  those  associated  with  the  above  mentioned  sections. 


2.2.3 . 1 . 1 . 1  Single-Phase  Minor  Losses 
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Most  of  the  analytical  and  experimental  expressions  for  head  losses  are  based  on 
single-phase  phenomena.  In  this  section  the  minor  losses  for  the  single-phase  are 
reviewed.  Then  two-phase  head  losses  are  reviewed  in  the  next  section. 

The  pressure  losses  due  to  sudden  contraction  and  enlargement  of  the  flow  area 
are  expressed,  respectively,  as  follows  [Fox  and  McDonald,  1992]. 

(2-19) 

K  =  K,  y  (2-20) 

where,  Kc,  Ke,  and  v  stand  for  the  minor  loss  coefficients  for  the  sudden  contraction  and 
enlargement,  and  the  flow  velocity  at  the  narrow  path  between  upstream  and  downstream 
of  the  flow,  respectively. 

The  minor  loss  coefficients  for  sudden  expansions  and  contractions  in  circular 
ducts  are  given  in  Figure  2.11  [Fox  and  McDonald,  1992],  Note  that  both  loss 

v2 

coefficients  are  based  on  the  larger  — .  Thus  losses  for  a  sudden  expansion  are  based  on 

v2  v2 

—  ,  upstream  condition,  and  those  for  a  contraction  are  based  on  —  ,  downstream 
2  2 

condition. 

There  are  also  head  losses  due  to  the  flow  entrance  and  exit.  The  entrance  losses 
occur  when  a  fluid  enters  a  conduit  from  a  large  tank  or  reservoir  as  shown  in  Figure  2.12 
[Evett  and  Liu,  1987].  The  amount  of  head  loss  is  significantly  dependent  on  the  shape 
of  the  entrance.  If  the  entrance  is  well-rounded,  the  entrance  loss  becomes  very  small. 
Entrance  loss  coefficients  can  be  found  in  Figure  2.13.  Exit  losses  occur  when  a  liquid 
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exits  a  conduit  and  enters  a  large  tank  or  reservoir,  as  shown  in  Figure  2.14.  The  exit  loss 
coefficient  is  1.0  for  all  cases  (i.e.,  regardless  of  the  shape  of  the  exit)  [Evett  and  Liu, 
1987]. 

Gradual  expansion  losses  occur  when  there  is  a  gradual  increase  in  conduit  size, 
and  gradual  contraction  losses  happen  when  there  is  a  gradual  decrease  in  conduit  size. 
Loss  coefficients  for  gradual  expansions  can  be  determined  from  Figure  2.15  [Evett  and 
Liu,  1987].  Table  2.2  illustrates  loss  coefficients  for  gradual  contractions. 

The  head  loss  of  a  bend  is  larger  than  for  fully  developed  flow  through  a  straight 
section  of  equal  length.  The  additional  loss  is  primarily  the  result  of  secondary  flow,  and 
is  represented  most  conveniently  by  an  equivalent  length  of  straight  pipe.  The  equivalent 
length  depends  on  the  relative  radius  of  curvature  of  the  bend.  Because  they  are  simple 
and  inexpensive  to  erect  in  the  field,  miter  bends  often  are  used,  especially  in  large  pipe 
systems  [Fox  and  McDonald,  1992].  Figure  2.16  shows  design  data  for  miter  bends. 


2.2.3. 1.1.2  Two-Phase  Minor  Losses 

As  the  C02  expands  in  the  piston-cylinder  device,  its  phase  changes  from  a 
supercritical  state  to  single-phase  state  (compressed  liquid  or  superheated  vapor)  and  then 
to  a  two-phase  state.  This  means  that  the  C02  becomes  two-phase  at  the  end  of  the 
expansion,  leaving  the  piston  device  and  entering  the  evaporator  as  a  two-phase  fluid.  In 
order  to  simulate  the  discharging  of  C02  from  the  expansion  device,  literature  on  head 
losses  in  two-phase  flow  was  reviewed.  Unfortunately,  one  of  the  least-studied  aspects  of 
two-phase  vapor-liquid  flow  is  that  of  the  change  in  static  pressure  and  the  energy  loss 
with  such  features  as  enlargements,  contractions,  orifices,  bends,  and  valves. 
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Consider  a  two-phase  flow  passing  through  a  sudden  enlargement  as  shown  in 
Figure  2.17.  The  general  method  used  in  single-phase  flow  is  applicable  to  a  one¬ 
dimensional  separated  two-phase  flow  [Collier,  1972].  Collier  discussed  that  the  pressure 
just  downstream  of  the  enlargement  (plane  1)  is  found  to  be  equal  to  that  at  the  end  of  the 
smaller  section  of  pipe  and  a  simplified  momentum  balance  for  the  combined  flow  yields: 

Pi^i  ~  Pi^i  ~  Wf  —  Ug\ )  W/  (m/2  —  uf  \  )•  (2-21) 

Substituting  for  Wg,  Wf,  ug  and  «/  in  terms  of  G,  x,  a,  v/,  and  vg  using  the  continuity 
relationships  and  defining  a  as  ratio  (Aj/A2),  Equation  (2-21)  becomes 


Pi  Pi  ~  f 


(2-22) 


where  W  and  u  are  the  mass  flow  rate  and  velocity  and  subscripts  g  and /denote  the  gas 
phase  and  liquid  phase,  respectively.  G  is  mass  velocity  and  is  defined  as  the  rate  of 
mass  flow  divided  by  the  flow  area: 

_  W  u 

G  —  —  =  pu  —  — .  (2-23) 

A  v 


x,  a ;  and  v  represent  the  mass  quality,  void  fraction,  and  specific  volume,  respectively. 

If  the  void  fraction  remains  unchanged  across  the  feature  (a\-ao=d).  Equation  (2- 
22)  becomes 


Pi  ~  Pi  =  G?cr(l-  cr)vf 


(2-24) 


Comparing  this  result  with  that  for  the  total  flow  assumed  to  be  single-phase  liquid 
reveals  that  the  term  in  the  square  brackets  acts  as  a  two-phase  multiplier.  For 
homogeneous  flow  Equation  (2-24)  reduces  to 
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Measurements  of  void  fraction  made  in  the  vicinity  of  sudden  enlargements  by  some 
other  authors  [Richardson,  1958,  for  horizontal  flow  and  Petrick  and  Swanson,  1959,  for 
vertical  flow]  indicate  that  the  assumption  ar,=ar2=«is  only  approximately  true. 

Figure  2.18  shows  a  sudden  contraction  in  the  cross-sectional  area  of  the  pipe 
two-phase  flow.  The  general  approach  used  for  single-phase  flow  can  again  be  applied  to 
a  two-phase  separated  flow  passing  through  a  sudden  contraction.  As  the  fluid  passes 
from  plane  1  to  plane  c,  at  the  vena  contracta  it  is  accelerated  and  pressure  energy  is 
converted  to  kinetic  energy  with  little  or  no  friction  dissipation  [Collier,  1972].  Beyond 
plane  c,  however,  the  conditions  are  similar  to  those  for  a  sudden  enlargement  and 
frictional  dissipation  occurs.  Collier  mentioned  that  equations  derived  above  could 
therefore  be  applied  to  this  ‘enlargement’  from  plane  c  to  plane  2  to  obtain  an  expression 
for  the  frictional  pressure  loss.  The  coefficient  of  contraction,  Cc,  is  used  to  denote  the 

A  g  'j 

ratj0  — ^and  replace  a,  likewise  the  mass  velocity  Gc  =  —■  replaces  G\  in  Equation  (2- 

4  V  CJ 

22).  Rearranging  and  simplifying  the  equations  yields  the  following  expression  for  the 
frictional  pressure  loss  across  a  sudden  contraction. 
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Then  the  change  in  static  pressure  at  a  sudden  contraction  can  be  given  by  the  sum  of  the 
frictional  dissipation  and  the  theoretical  kinetic  energy  change.  For  a  homogeneous 
model  [Collier,  1972] 
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(2-27) 


where,  v/,  v/g,  and  x  represent  the  specific  volume  of  liquid  phase,  specific  volume 
difference  between  vapor  and  liquid  phases,  and  vapor  quality,  respectively.  In  the  above 


equation,  the  coefficient  of  contraction,  Cc ,  is  a  function  of  a 


'.A' 

V 


Table  2.3  shows 


the  relation  between  Cc  and  <rfor  turbulent  single-phase  flow. 

Collier  said  that  Ferrel  and  McGee  [Ferre]  and  McGee,  1966]  and  Geiger  and 
Rohrer  [Geiger  and  Rohrer,  1966]  found  that  the  homogeneous  model  (Equation  (2-27)) 
satisfactorily  represented  their  experimental  data  for  steam-water  through  sudden 
contractions  [Collier,  1972].  Again  measurements  of  void  fraction  made  by  some  other 
authors  [Richardson,  1958,  for  horizontal  flow  and  Petrick  and  Swanson,  1959,  for 
vertical  flow]  in  the  vicinity  of  sudden  contractions  indicate  that  the  assumption  implicit 
in  Equation  (2-27),  or,  =  a2  =  a,  is  only  approximately  true  [Collier,  1972]. 


2.2. 3. 1.2  Piston  Ring  Scuffing  and  Lubrication 

Piston  ring  behavior  must  be  characterized  to  properly  model  the  expansion 
device  producing  work.  If  there  is  a  separating  oil  film  between  the  rubbing  surfaces, 
there  is  no  piston  ring  contact  with  the  wall,  or  “scuffing”.  Scuffing  is  a  consequence  of 
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the  contact  pressure  exceeding  the  pressure  the  momentary  local  lubricating  conditions 
can  bear  [Aue,  1975]. 

2.2.3. 1.2.1  Piston  Ring  Scuffing 

Aue  investigated  the  piston  ring  behavior  experimentally  for  internal  combustion 
engine  applications  [Aue,  1953].  He  developed  a  test  facility  to  perform  investigations  of 
piston  rings  when  running  under  various  pressures  and  lubricating  conditions.  In  the  first 
test,  lubrication  was  adjusted  to  correspond  to  normal  combustion  engine  conditions.  The 
conclusion  of  the  test  was  that  even  an  extremely  thin  lubricating  film,  so  long  as  it  was 
still  present,  would  prevent  metal-to-metal  contact  and  the  resulting  scuffing.  Other 
testing  focused  on  a  complete  stationary  set  of  five  piston  rings  on  a  piston  with  constant 
and  pulsating  pressure  loading.  The  result  of  the  test  was  that  the  last  piston  ring  alone 
took  most  of  the  pressure  drop  across  the  piston  with  both  of  well-fitting  rings  and  poor¬ 
fitting  rings. 

Aue  also  examined  the  behavior  of  a  piston  ring  theoretically  [Aue,  1974].  He 
mentioned  that  wear  would  correct  minor  imperfections  and  a  little  wear  in  service  would 
stabilize  operation.  He  suggested  rotation  of  the  piston  about  its  axis  to  improve  the 
running  condition. 

Wakuri  and  Ono  carried  out  experimental  studies  on  the  abnormal  wear  of 
cylinder  liners  and  piston  rings  in  large-sized  marine  Diesel  engines  [Wakuri  and  Ono, 
1973].  The  authors  investigated  the  mechanism  of  abnormal  wear,  detecting  procedures, 
and  the  effect  of  some  factors  on  the  occurrence  of  abnormal  wear.  Some  results  of  their 
studies  were  similar  to  those  of  work  done  by  Aue.  Abnormal  wear  is  initiated  by  the 
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scuffing  caused  by  serious  edge-loading  of  piston  rings  and  is  accelerated  by  the  abrasive 
wear  due  to  very  large-sized  worn  particles  scraped  off  from  the  cylinder  liner.  They 
discussed  that  a  roughened  surface  had  much  more  resistance  to  scuffing  than  a  smooth 
surface  and  that  rotation  of  a  piston  was  more  effective  to  increase  the  resistance  to 
scuffing. 


2. 2. 3. 1.2. 2  Piston  Ring  Lubrication 

If  there  is  a  separating  oil  or  fluid  film  between  the  rubbing  surfaces,  it  prevents 
scuffing.  In  general,  the  starting  point  for  analysis  of  the  piston  ring  lubrication  is  the 
Reynolds  equation.  The  derivation  of  Reynolds  equation  is  given  in  Appendix  A. 

Baker  et  al.  developed  a  computer  program  to  predict  lubricant  film  thickness, 
friction  for  a  piston  ring  and  pressure  differences  across  a  piston  ring  face  over  the  crank 
angle  [Baker  et  al.,  1973].  The  authors  discussed  that  serious  reduction  in  film  thickness 
may  be  expected  if  the  radius  of  curvature  of  the  piston  ring  falls  much  below  10'2  inches. 

Rohde  and  Whitaker  presented  a  mathematical  model  of  the  hydrodynamic 
lubrication  of  piston  rings  [Rohde  and  Whitaker,  1979].  The  authors  examined  the 
effects  of  engine  speed,  engine  geometry,  lubricant  viscosity,  piston  ring  face  profile, 
tension,  and  thickness.  They  calculated  the  power  loss  by  multiplying  friction  force  by 
piston  velocity.  They  did  not  consider  the  term  of  film  squeeze  in  the  power  loss 
calculation. 

Sun  developed  a  subroutine  package  to  solve  the  inverse  problem  of  slider 
lubrication,  i.e.,  load  being  given  and  the  minimum  film  thickness  being  sought  [Sun, 
1990].  The  author  mentioned  that  the  problem  might  also  contain  all  or  any  combination 
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of  finite-width,  squeeze  action,  and  oil  starvation.  The  author  calculated  the  power  loss 
by  considering  both  the  friction  obtained  by  integrating  the  shear  stress  and  squeeze  of 
film. 

Jeng  performed  theoretical  analysis  for  piston  ring  lubrication  by  considering 
fully  flooded  lubrication  and  starved  lubrication  [Jeng,  1992].  He  developed  a  system  of 
nonlinear  equations  formulated  based  on  the  hydrodynamic  lubrication  theory.  He  solved 
the  system  of  equations  using  a  modified  Newton’s  method.  He  demonstrated  the 
capability  of  his  computer  model  by  investigating  the  effect  of  operating  conditions,  ring 
tensions,  and  ring  face  profile  for  the  fully  flooded  and  starved  inlet  conditions.  He 
showed  that  the  starved  lubrication  occurred  when  the  lubricant  supply  was  not  adequate 
for  fully  flooded  lubrication. 

Radakovic  and  Khonsari  investigated  the  effects  of  viscosity  variations  with 
temperature,  in  conjunction  with  the  non-Newtonian  shear  thinning  behavior  of 
multigrade  oils  on  the  hydrodynamic  lubrication  of  the  piston  ring  [Radakovic  and 
Khonsari,  1997].  Their  simulation  results  showed  that  shear  thinning  could  have  a 
significant  effect  on  parameters  such  as  film  thickness,  viscous  drag  force,  and  power 
loss.  The  thermal  effects  from  viscous  dissipation  in  the  clearance  space  along  the  piston 
ring  also  influenced  these  parameters,  but  to  a  lesser  degree. 


2232  Fluid  Mechanic  Model  for  the  ED-WOW  Device 

A  theoretical  model  was  developed  to  predict  the  efficiency  of  the  ED-WOW 
based  on  a  piston-cylinder  type  device.  The  overall  model  is  composed  of  several  sub¬ 
models.  The  main  modeling  effort  focused  on  frictional  losses  between  the  piston  rings 
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and  cylinder  walls  and  on  loss  due  to  squeeze  of  film.  The  model  is  based  on  the 
lubrication  theory  and  assumes  that  the  piston  ring  in  the  device  is  a  slipper-pad  bearing. 
The  model  predicts  the  film  thickness  of  the  lubricant,  in  this  case  either  CO2  or  CO2  and 
oil.  The  pressure  distribution  under  the  piston  ring,  power  loss,  and  work  done  by  the 
CO2  are  also  calculated. 

Figure  2.19  shows  the  flow  diagram  of  the  model.  As  the  first  step,  the 
displacement  and  velocity  of  the  piston,  crank  angle,  and  the  work  done  by  the  CO2  are 
calculated  for  a  given  time  step.  The  power  loss  and  frictional  force  are  calculated  based 
on  lubrication  theory.  In  addition,  the  CO2  mass  flow  into  or  out  of  the  device  (based  on 
the  process)  is  predicted.  Finally,  the  program  calculates  the  pressure  and  temperature  of 
the  CO2  and  stores  the  results  for  the  given  time  step.  The  same  calculations  are  repeated 
for  the  next  time  step  and  so  on  until  the  device  goes  through  two  full  revolutions  (720°). 
Once  the  program  concludes  the  calculations  for  the  two  revolutions,  global  parameters, 
such  as  total  power  output,  power  loss,  and  efficiency  are  calculated.  This  model  is 
referred  to  as  ‘ED-WOW’  (Expansion  Device  With  Output  Work  for  CO2  cycle). 


2.2. 3. 2.1  Modeling  for  Mass  Flow  into  Chamber 

The  theoretical  considerations  of  the  various  minor  losses  are  quite  complicated; 
therefore,  minor  losses  are  usually  evaluated  by  empirical  methods.  Figure  2.10  shows 
configuration  of  the  connection  of  a  solenoid  valve  to  fittings  and  chamber  of  cylinder. 
The  mass  flow  into  the  chamber  of  the  cylinder  is  calculated  based  on  the  configuration 
of  Figure  2.10.  Figure  2.20  shows  geometries  causing  minor  losses.  The  assumptions  for 


this  model  are  as  follows. 
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•  Pressure  losses  due  to  friction  (other  than  minor  losses)  are  negligible. 

•  Potential  energy  change  is  negligible. 

•  If  two-phase  flow  occurs,  it  is  assumed  to  be  homogeneous  flow.  Therefore, 
the  single-phase  formula  can  be  applied  with  the  averaged  properties. 

•  The  chamber  is  large  compared  to  the  flow  path  dimension  so  that  it  is 
considered  as  a  reservoir. 

•  The  flow  is  isothermal. 


Four  major  contributions  for  the  pressure  losses  are  considered.  They  are  pressure  losses 
due  to  sudden  contraction,  sudden  enlargement,  miter  bend,  and  pressure  loss  at  exit. 


2. 2. 3.2. 1 . 1  Pressure  Loss  due  to  Sudden  Contraction 

The  pressure  loss  due  to  sudden  contraction  of  flow  area  is  given  as  [White, 


1986] 


4Pl-»2  K 


pv 


(2-28) 


The  equation  above  can  be  expressed  as 


Pi  P  2 
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PlV2 


(2-29) 


CD  is  a  discharge  coefficient  and  is  determined  experimentally  by  comparing  the  filling 
behavior  predicted  by  computer  model  to  the  experimental  data.  It  is  discussed  in 
Section  2.2.3.2.3. 


2.2.3.2.1.2  Pressure  Loss  due  to  Miter  Bend 
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It  is  assumed  that  the  flow  path  is  bent  by  90°  without  change  of  flow  area  as 
shown  in  Figure  2.20.  This  type  of  bend  is  known  as  ‘miter  bend’.  Pressure  loss  due  to 
bending  is  described  in  terms  of  friction  loss  as  follow  [Fox  and  McDonald,  1992]: 


De  p2v\ 

APw  =  f  P2~P^ 


(2-30) 


The  value  of  is  60  for  90°  miter  bends  as  shown  in  Figure  2.16.  The  friction  loss 


e  . 


coefficient  is  determined  by  assuming  single-phase  flow.  —  is  about  0.0032  for 


commercial  steel  for  D  =  — .  The  Reynolds  number  is  given  as 
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Re  = 


pvD 


(2-31) 


where  D  is  diameter  of  flow  channel  and  p,  v,  and  p  represent  the  density,  velocity,  and 
viscosity  of  the  fluid  in  the  flow  channel,  respectively.  In  the  case  of  CO2  between  0°C 

N  •  s 

and  100°C,  the  value  of  the  viscosity  is  about  1.5x10  5  — 5-  [Fox  and  McDonald,  1992]. 

m 

When  the  flow  is  laminar  (Re  <  2000),  the  friction  coefficient  is  calculated  as 


Re 


(2-32) 


If  the  flow  is  turbulent  (Re  >  2000),  it  is  determined  from  Moody  diagram  and  is  about 


0.025. 


2.2. 3. 2. 1.3  Pressure  Loss  due  to  Sudden  Enlargement 
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The  pressure  loss  due  to  sudden  enlargement  of  flow  area  is  given  as  [White, 


1986] 


=  K< 
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(2-33) 


The  equation  above  can  be  expressed  as 


Pl~P4  =  CD,e 
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CD,e  is  a  discharge  coefficient  for  sudden  enlargement  and  is  generally  given  as  1.0. 
Therefore,  the  pressure  loss  due  to  sudden  expansion  is 

(  .aY 

V  Ay 


Pi  ~  Pa  = 


PiV3 


(2-35) 


2.2.3.2.1.4  Pressure  Loss  due  to  Exit 

Pressure  loss  at  exit  to  a  reservoir  is  expressed  as 

(2-36) 

where,  the  loss  coefficient,  Kexih  is  given  as  1.0  for  all  shape  of  exit  [Evett  and  Liu,  1987], 


2.2.3. 2.2  Calculation  Method  for  Mass  Flow  into  Chamber 

Total  pressure  loss  from  the  outlet  of  the  solenoid  valve  to  the  chamber  consists 
of  all  of  the  losses  given  above.  It  means  that  total  pressure  loss  is  given  as 


A ptol  =  A p^2  +  Apbend  +  Ap^4  +  Ap4_,5 . 


(2-37) 


If  the  total  pressure  loss  is  known,  the  mass  flow  rate  can  be  calculated  using  Equation 
(2-37).  When  p\  and  p5  are  given,  the  calculation  procedure  is  as  follows: 
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Estimate  p\  =  £l±£l  . 
2  2 


•  Find  pi  with  assumption  of  isothermal  process  ( T2  =  7] ). 

•  Calculate  vi  by  using  Equation  (2-29).  Co  will  be  determined  later  based  on 
experiment. 

•  Calculate  pi  by  using  Equation  (2-30). 

•  Calculate  p 4  by  using  Equation  (2-35). 

•  Calculate  /94  and  v4  by  assuming  isothermal  process  ( T4  =  7] ) . 

Check  if  Equation  (2-37)  is  satisfied  with  the  given  and  predicted  values.  If  it  is  not 
satisfied,  pi  (or  p^1)  for  the  next  trial  is  determined  by  using  ‘Bisection  method’ 
[Hoffman,  1992]  as  follows: 


_n+i .  pI  +  pT 
Pi  ■  2 


(2-38) 


Repeat  the  calculation  until  Equations  (2-29),  (2-30),  (2-34),  (2-35)  and  (2-36)  are  all 
satisfied. 

Mass  flow  rate  can  be  calculated  by  using  mass  conservation  as  follow: 


™in  =  or  min  =  P3V3A3  = ... 


(2-39) 


Figure  2.21  shows  the  flow  diagram  for  the  mass  calculation. 


2. 2. 3. 2. 3  Determination  of  Discharge  Coefficient  of  Sudden  Contraction 
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As  mentioned  before,  the  discharge  coefficient  for  sudden  contraction  is 
determined  experimentally.  Experiments  were  performed  to  investigate  the 
characteristics  of  the  solenoid  valve  as  mentioned  in  Chapter  3.  The  results  of  the 
experiments  were  used  to  determine  the  coefficient.  The  discharge  coefficient  is  obtained 
based  on  Equation  (2-37)  but  not  determined  exactly  and  directly  from  a  separate 
experiment  because  no  experiment  was  conducted  for  that. 

Figure  2.22  shows  the  comparison  of  predicted  filling  behaviors  of  C02  to  the 
experimental  data.  The  y-axis  represents  dimensionless  pressure.  The  measured  pressure 
is  dynamic  pressure  but  the  predicted  pressure  is  absolute  pressure.  The  experimental 
curve  shows  the  decrease  of  pressure  after  a  peak  of  pressure  because  of  its  dynamic 
pressure.  It  is  considered  that  there  is  pressure  equilibrium  in  the  chamber  at  the  peak 
pressure  point.  When  CD  is  1.0,  the  time  at  which  the  predicted  pressure  reaches  a  peak 
is  far  away  from  the  experimental  data.  Trial-error  method  was  used  to  obtain  Co  of  6.5 
to  give  the  same  time  for  pressure  equilibrium  in  the  chamber,  as  shown  in  Figure  2.22. 


2.2.3. 2.4  Modeling  for  Mass  Flow  Out  from  Chamber 

At  the  end  of  the  expansion  process  the  C02  pressure  is  expected  to  be  about 
3.4MPa  (500  psig)  and  the  state  will  be  two-phase.  In  order  to  predict  the  mass  flow  out 
from  the  chamber,  a  two-phase  flow  model  is  required  for  sudden  contraction.  Flow  area 
is  assumed  to  decrease  suddenly  when  C02  is  discharged  from  the  chamber  to  the  inlet  of 
the  evaporator.  Figure  2.18  shows  the  sudden  contraction  in  cross-sectional  area  of  the 
pipe  with  two-phase  flow.  A  homogeneous  two-phase  flow  model  for  sudden  contraction 
is  expressed  [Collier,  1972]  as  follows: 
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(2-40) 


A  detailed  explanation  of  the  above  equation  is  given  in  Section  2. 2. 3. 1.1. 2. 

If  the  pressures  in  and  downstream  of  the  chamber  are  known,  the  mass  velocity, 
Gi,  can  be  calculated.  As  shown  in  Equation  (2-23),  the  mass  velocity  is  expressed  as 


^  _  W  _  „  _u2 

G2  —  —  p2u2  — 

A  v2 


Therefore,  the  mass  flow  rate  from  the  chamber  is  calculated  as 

t^out  =  PiuiA  —  G2A2. 


(2-41) 


(2-42) 


Subscript  ‘2’  represents  the  downstream  of  the  chamber  shown  in  Figure  2.18. 


2.2.3. 3  Modeling  for  Piston  Ring  Lubrication 

As  a  piston  moves  up  and  down,  there  is  loss  of  energy  due  to  the  friction, 
squeeze  of  the  piston  ring,  heat  transfer,  blowby,  etc.  In  this  section  the  losses  due  to  the 
friction  and  squeeze  motion  of  the  piston  ring  are  modeled.  Figure  2.23  shows  a 
schematic  of  the  piston  ring  and  cylinder  wall  conjunction. 


2. 2. 3. 3.1  Theoretical  Model 


2. 2. 3. 3. 1.1  Governing  Equation 
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The  piston  ring  can  be  treated  as  a  reciprocating,  dynamically  loaded  bearing  with 
a  combined  sliding  and  squeeze  motion  [Jeng,  1992].  As  discussed  in  Appendix  A,  the 
Reynolds  equation  for  an  incompressible  lubricant  can  be  written  as: 

l±(hidp)  =  6U^  +  12V  (t)  (2-43) 

pi  dx  {  dxj  ox 

where  pi,  h,  and  U  stand  for  the  viscosity,  film  thickness  of  the  lubricant,  and  the  piston 
speed,  respectively,  and  V (t)  represents  the  velocity  of  piston  ring  in  the  y-coordinate  (h- 
coordinate  in  this  case)  and  can  be  expressed  as: 


(2-44) 


2.2.33.2  Geometry  of  Piston  Ring 

The  film  thickness  can  be  expressed  as: 


h{x,t)  =  hm(t)+hs{x)  (2-45) 

where  hm  is  the  time  varying  minimum  film  thickness  and  hs  is  the  ring  face  profile.  A 
parabolic  form  for  the  ring  face  profile  is  sufficiently  accurate  [Jeng,  1992].  Therefore, 
the  ring  face  profile  can  be  written  with  the  assumption  of  no  offset  of  the  piston  ring  as. 


K(x)  = 


(2-46) 


where  CR  is  the  crown  height  and  b  is  the  width  of  the  ring  as  shown  in  Figure  2.23. 


2.2.3.3.3  Boundary  Conditions 
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A  fully  flooded  inlet  condition  is  assumed,  that  is,  the  lubricant  supply  is  always 
adequate  at  the  inlet.  The  inlet  and  exit  boundary  conditions  are  expressed  as: 


V  2  . 


=  Pbo, 


if) 


(2-47) 


(2-48) 


where  pb0,  and  ptop  stand  for  the  pressure  at  the  trailing  and  leading  edges  of  the  ring, 
respectively.  The  values  of  pboi  and  ptop  depend  on  the  direction  of  motion  and  piston 
position.  For  example,  when  the  piston  moves  to  the  direction  of  the  right  hand  side  in 
Figure  2.23  (during  the  down  stroke,  intake  and  expansion  stroke),  pb0t  is  the  CO2 
pressure  in  the  cylinder  chamber  and  ptop  the  pressure  between  the  first  and  second 
compression  rings. 


2. 2. 3. 3. 4  Initial  Condition 

An  initial  condition  is  needed  to  solve  the  problem  of  dynamically  loaded 
bearings.  It  is  obtained  arbitrarily  by  assuming  that 


=  0. 


(2-49) 


2. 2. 3.3. 5  Piston  Ring  Elastic  Characteristics 

The  piston  ring  is  assumed  to  be  ideal,  that  is,  the  ring  is  truly  circular  and  exerts 
uniform  pressure  all  around  the  cylinder  wall.  It  is  also  assumed  that  the  ring  only  moves 
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radially  as  a  unit  under  loading.  Therefore,  the  piston  ring  can  be  treated  as  a  thin  walled 
cylinder.  Then,  the  radial  deflection  of  the  ring  is 

C  =  (Pm -pY  (2-50) 

Etr 

where  pm  is  the  external  mean  hydrodynamic  pressure,  p'  is  the  uniform  internal  gas 
pressure,  E  is  the  Young’s  modulus  of  the  ring,  r  is  radius  of  the  ring,  and  tr  is  the  ring 
thickness.  The  radius  of  the  ring  is  assumed  to  be  half  of  the  diameter  of  the  cylinder 

wall  f  =  —  I .  The  radial  deflection  of  the  ring  is  equivalent  to  the  minimum  film 

l  2 ) 


thickness,  hm,  and  the  minimum  film  thickness  is  written  as: 

h  =  ( Pm  ~  P'V  (2-51) 

m  Etr 

Equation  (2-50)  can  be  expressed  differently  after  mathematical  manipulation  with  some 
definitions  as  follows: 


(2-52) 


where  pei  is  the  ring  elastic  pressure,  T  is  the  tangential  ring  tension  and  B  is  the  cylinder 


bore  diameter. 


2.2.3.3.6  Method  of  Solution 


2. 2. 3. 3.6.1  Pressure  Distribution 

Integration  of  Equation  (2-43)  with  respect  to  x  gives 
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h3(x,t)^-  =  6pU(t)h(x,t)  +  12ju  dh(x,t.l x  +  C, . 
dx  dt 

where  C\  is  a  constant.  The  time  derivative  of  h  can  be  expressed  as: 

dh(x,t)  _  d[hm(t)+h5{x)]  _  dhm 
dt  dt  dt 

Integration  of  Equation  (2-53)  with  respect  to  x  gives 


pM=12Afc(b) 


U(t)h(s,t)  ,  dh 


H - —  s  +  C? 

dt  2 


Ids  +  C, 


(2-53) 


(2-54) 


(2-55) 


where  C2  is  a  constant.  By  applying  the  boundary  conditions  the  pressure  distribution  is 
obtained  as: 


p{x-')=UM\i?h)i 


U{t)h(s,t)  dhi 


H - —  S  +  C2 

dt  2 


^  +  Pbo 


(2-56) 


and  the  integral  constant  C2  is 


"*rop  Pbol 

12/7 


_T_J_ 


U(t)h(s,t)  |  „ 

2  dt 


ds 


(2-57) 


f— 


As  shown  in  Figure  2.23,  the  profile  of  the  ring  face  is  assumed  to  be  parabolic 
and  there  are  convergent  and  divergent  portions  of  the  gap.  Cavitation  may  occur  in  the 
divergent  portion.  Therefore,  the  pressure  is  assumed  to  be  constant  in  the  lubricant  in 
the  divergent  portion  of  the  gap  if  a  cavitation  occurs  there.  Figure  2.24  shows  the 
pressure  distribution  under  the  piston  ring  with  this  assumption.  The  pressure  in  the 
divergent  portion  is  cut-off  and  assumed  to  be  pbot  during  the  intake  and  expansion  stroke. 
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The  x-coordinate  was  normalized  by  dividing  the  distance,  x,  by  the  half  thickness  of  the 
piston  ring,  b/2,  as  shown  in  Figure  2.23. 


2.2.3.3.6.2  Friction  Force 

Velocity  profile  in  the  y-coordinate  is  given  by  Equation  (A- 15)  in  Appendix  A 


(A- 15) 


The  shear  stress  acting  on  the  piston  ring  is  derived  using  the  above  velocity  profile  as: 

(2-58) 


du 


y=h 


h  dp  juU 
2  dx  h 


The  friction  force  for  hydrodynamic  lubrication  can  be  written  as: 

hdp 


F=  fr<M  =  f-^-  — 
J  Jl2  3bc  h 


dA 


(2-59) 


where  A  represents  the  area  of  the  lubricant  region. 

In  the  case  of  boundary  lubrication,  the  friction  force  on  the  ring  is  calculated  by: 

F  =  pc  {nB  W  )sign  ( U )  (2-60) 

where  /4,  B  and  Ware  the  coefficient  of  friction,  the  cylinder  bore,  and  the  load  imposed 

on  the  ring,  respectively.  W  is  expressed  as: 

W=b[pel+pJ  (2-61) 

where  pei  and  pbk  are  the  ring  elastic  pressure  of  Equation  (2-52)  and  the  pressure  at  the 
back  side  of  the  ring,  respectively. 
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2.2.33.6.3  Power  Loss 

When  hydrodynamic  lubrication  is  present  between  the  piston  ring  and  cylinder 
wall,  the  power  loss  is  [Jeng,  1992] 


P=  \^-U2dA  +  J- 


h 3  (dp''2 


12  fJL 


\dxj 


dA. 


(2-62) 


This  equation  automatically  includes  the  power  loss  due  to  squeeze  motion. 
In  the  case  of  boundary  lubrication,  the  power  loss  is  calculated  by: 


(2-63) 


2.2.33.7  Numerical  Scheme 


Equations  (2-51)  and  (2-56),  the  film  thickness  and  pressure  distribution  in  the 
lubricant,  are  coupled  to  each  other  and  are  a  system  of  nonlinear  equations.  The  system 
is  solved  iteratively  by  using  the  internal  halving  method  [Hoffman,  1992],  First,  the 
minimum  film  thickness,  hm,  is  determined  by  using  the  internal  halving  method.  Let  this 
minimum  film  thickness  be  hm\.  Then,  the  time  derivative  of  hm,  the  squeeze  velocity, is 
calculated  by  using  the  values  of  the  current  and  previous  steps  as: 


dh"  hnm  -  hn~l 

m  m  m 

At 


dt 


(2-64) 


where  superscripts  n  and  n- 1  represent  current  and  previous  time  steps,  respectively,  and 
At  stands  for  the  time  interval.  Now  Equation  (2-56)  can  be  solved  because  there  are  no 
unknown  variables  in  that  equation.  The  pressure  in  the  lubricant  can  be  determined  after 
double  integrations  as  shown  in  Equations  (2-56)  and  (2-57).  Analytical  expression  for 
each  term  of  the  integrations  in  Equations  (2-56)  and  (2-57)  have  been  developed  and 
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shown  in  Appendix  B.  The  external  mean  hydrodynamic  pressure,  pM,  is  calculated 
numerically  by  using  Simpson’s  1/3  rule  [Hoffman,  1992]  as: 

b 

p“=\\pdx  =  \ 

2 

where  Ax  and  subscript  n  stand  for  the  mesh  size  and  the  number  of  meshes,  respectively. 
Then  hm  is  calculated  again  using  Equation  (2-51)  and  let  this  be  hm 2.  The  two  minimum 
film  thicknesses,  hm\  and  hm 2,  are  compared  to  check  if  the  Equations  (2-51)  and  (2-56) 
are  related  correctly.  This  procedure  is  terminated  when  the  periodic  condition  is 
satisfied,  i.e.,  when: 

Ki-hJ<£  (2'66) 

where  £  is  the  convergence  tolerance.  If  the  solved  hm  is  less  than  the  composite  rms 
roughness,  cr,  it  is  assumed  that  the  lubrication  mode  is  boundary  lubrication  and  hm  is 
equal  to  the  composite  rms  roughness.  In  the  case  of  hydrodynamic  lubrication,  friction 
force  and  power  loss  are  calculated  numerically  by  using  Simpson’s  1/3  rule. 

2.2.3. 3. 8  Validation  of  the  Model  for  the  Piston  Ring  Lubrication 

This  model  for  the  piston  ring  lubrication  has  been  tested  by  simulating  an 
internal  combustion  engine  of  2000  rpm  revolutionary  speed.  Its  results  have  been 
compared  to  the  results  of  [Jeng,  1992]  and  [Radakovic  and  Khonsari,  1997]  who 
performed  lubrication  analysis  for  the  same  engine.  Table  2.4  shows  the  values  of  the 
parameters  of  the  piston  ring  and  operating  conditions  for  the  simulated  engine.  Figures 


\>Ax{p0  +4/?,  +  2p2  +  4/?3  +  •  •  •  +  4pn_1  +p„) 


(2-65) 
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2.25  and  2.26  illustrate  the  piston  speed  of  the  engine  as  a  function  of  crank  angle  and  the 
combustion  pressure  of  the  engine  versus  crank  angle,  respectively. 

Figure  2.27  shows  the  minimum  film  thickness  predicted  by  the  models  of  Jeng 
and  Radakovic  and  Khonsari.  In  the  legend  of  the  figure,  ‘current  analysis’  notation 
stands  for  the  prediction  by  Radakovic  and  Khonsari  and  ‘Y.  R.  Jeng’  indicates  the 
prediction  by  the  model  of  Jeng.  Figure  2.28  illustrates  the  power  losses  predicted  by  the 
model  of  Radakovic  and  Khonsari.  They  considered  four  kinds  of  lubricants  in  their 
analysis  as  shown  in  the  legend  of  the  figure.  Among  them  ‘Newtonian/Isothermal’ 
curve  is  compared  to  the  power  loss  predicted  by  this  model  because  of  the  same 
characteristics  of  the  lubricant  simulated  in  this  analysis. 

Figures  2.29  and  2.30  show  the  predicted  minimum  film  thickness  and  power  loss 
by  this  model,  respectively.  As  shown  in  the  figures,  in  general,  the  trends  and 
magnitude  are  very  similar  to  those  of  Jeng,  and  Radakovic  and  Khonsari.  The 
differences  between  exact  magnitudes  of  the  minimum  film  thickness  and  power  loss 
seem  to  be  due  to  the  difference  in  the  treatment  of  the  piston  ring  and  handling  some 
parameters  such  as  the  pressure  distribution  under  the  piston  ring.  Jeng  used  a  value  for 
tangential  ring  tension  to  define  the  elastic  deflection  of  the  ring  but  the  present  analysis 
treats  the  piston  ring  as  a  thin  ring  subject  to  internal  and  external  pressure  which  will 
deflect  according  to  its  elastic  modulus.  The  present  analysis  assumes  that  the  pressure  in 
the  divergence  section  of  the  gap  is  the  same  as  the  pbot,  pressure  at  the  trailing  edge  of 
the  ring.  Jeng  considered  no  pressure  at  all  if  a  cavitation  occurred  under  ring  and 
Radakovic  and  Khonsari  assumed  that  the  pressure  in  the  divergence  section  of  the  gap 
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was  the  same  as  the  pressure  at  the  throat  (the  minimum  flow  area).  As  mentioned 
before,  however,  the  overall  trends  are  in  a  good  agreement. 


2.2.3.4  Modeling  for  Thermodynamic  Processes 

There  are  three  processes  within  one  revolution  of  crankshaft:  intake,  expansion, 
and  exhaust.  As  discussed  in  Chapter  3,  the  intake  valve  opens  and  the  supercritical  high 
pressure  C02  flows  into  the  cylinder  and  pushes  the  piston  down  from  the  Top  Dead 
Center  (TDC)  to  the  location  of  expansion  initiation  (LOE)  during  the  intake  process.  All 
valves  close  as  the  piston  reaches  the  LOE  and  the  expansion  process  starts  with  constant 
C02  mass.  When  the  piston  reaches  the  Bottom  Dead  Center  (BDC),  the  exhaust  valve 
opens  and  C02  is  swept  out  as  the  piston  moves  up  during  the  exhaust  stroke.  This 
exhaust  process  continues  until  the  piston  reaches  TDC.  The  same  processes  are  repeated 
during  the  next  revolution  of  the  crankshaft  of  the  device.  This  calculation  continues 
until  the  crank  angle  becomes  720°  (two  revolutions  of  the  crankshaft). 


2.2.3.4.1  Intake  Process 

The  intake  valve  opens  during  the  intake  process  and  C02  flows  into  the  chamber 
as  shown  in  Figure  2.31.  In  the  modeling  of  the  intake  process,  the  kinetic  and  potential 
energies  are  neglected  and  the  thermodynamic  properties  are  assumed  to  be  uniform  in 
the  chamber.  No  heat  transfer  is  assumed  (adiabatic  flow)  but  the  power  loss  due  to  the 
friction  and  squeeze  motion  of  the  piston  ring  is  considered  as  heat  loss.  The  last 
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assumption  in  the  intake  model  is  that  the  pressure  in  the  chamber  is  constant  at  each  time 
step  so  that  the  work  done  by  CO2  [Wark,  1995]  can  be  calculated  as 

V2 

w  =  -  \pdv  =  -Pl  (V2  -  Vt )  =  -Pl  A  V .  (2-67) 

v, 


2. 2. 3. 4.2  Calculation  Procedure 

There  are  two  steps  for  the  calculation  of  mass  flow  rate  and  thermodynamic 
properties  at  each  time  step  in  this  process.  First,  the  mass  flow  rate  is  calculated  by 
using  the  method  in  Section  2. 2.3. 2.1  with  the  pressure  difference  between  upstream  of 
the  intake  valve  and  chamber.  In  this  calculation  the  piston  does  not  move  i.e., 
V,  =  constant .  The  amount  of  C02  flowing  into  the  chamber  becomes 

min  =  ™in  ■ At  •  (2-68) 

In  the  second  step,  the  piston  moves  down  with  the  new  mass  in  the  chamber 
(m  =  constant)  for  the  given  At  and  V2  is  determined.  Then,  after  calculating  the  work 
done  by  the  carbon  dioxide  by  using  Equation  (2-67)  and  power  loss  by  using  Equation 
(2-62)  or  (2-63),  the  thermodynamic  properties  are  determined.  If  the  predicted  pressure 
is  higher  than  or  equal  to  that  upstream  of  the  valve  (p>2  >  Pin ) ,  the  pressure  in  the 
chamber  is  assumed  to  be  the  same  as  the  inlet  pressure  and  constant  during  the  rest  of 
intake  process. 


2. 2. 3.4. 3  Determination  of  Thermodynamic  Properties 
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When  the  pressure  in  the  chamber  is  lower  than  the  inlet  pressure,  the  1st  law  of 
thermodynamics  is  expressed  as: 

^isL  =  Q  +  W  +  mlnh„.  (2-69) 

dt 

The  above  equation  can  be  written  by  multiplying  both  of  the  left  hand  side  (LHS)  and 
right  hand  side  (RHS)  of  the  equation  by  dt  as: 

A E  =  A U  =  m2u2  -  mxux  =  -P  +  W  +  minhin  (2-70) 

where  the  subscripts  ‘1’,  ‘2’,  and  ‘in’  represent  the  states  before  and  after  the  piston 
movement,  and  inlet  state,  respectively.  /nin,  P,  and  W  stand  for  the  mass  flowed  into  the 
chamber,  power  loss,  and  work  done  by  C02,  respectively.  The  specific  internal  energy 
and  density  of  C02  at  state  ‘2’  can  be  obtained  as: 

u2  =  —  [mxux  -  P  +  W  +  minhin )  and  (2-7 1 ) 

m2 


The  pressure  and  temperature  in  the  case  of  single-phase  or  the  pressure  and  quality  in 
the  case  of  two-phase  are  determined  by  using  the  thermodynamic  property  data  [Span 
and  Wagner,  1996]  for  C02  with  Equations  (2-71)  and  (2-72). 

When  the  pressure  in  the  chamber  becomes  equal  to  the  inlet  pressure,  the  method 
in  Section  2.2.3.2.1  can’t  be  applied  to  predict  the  amount  of  mass  flowed  into  the 
chamber  because  the  pressure  in  the  chamber  is  assumed  to  be  constant  and  there  is  no 
pressure  difference  between  the  upstream  of  the  intake  valve  and  the  chamber  during  the 
rest  of  the  intake  process.  The  1st  law  of  thermodynamics  in  this  case  is 
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A U  =  m2u2  -  »!]«,  =  -P  +  W  +  (m2  -  ml  )hin .  (2-73) 

Unknowns  are  m2  and  U2  in  Equation  (2-73)  but  m2  can  be  expressed  as: 

m2  ~  V2 Pi  ■  {2-14) 

The  Secant  method  is  used  to  find  the  temperature  that  is  satisfied  with  ui  and  pi  for  the 
given  constant  pressure  ( p2  =  pin  =  constant)  .  Let’s  define  a  function  F{T)  from 
Equation  (2-73)  as  follows: 

F(T)=V2p2(Tlu2(T)-h, .]-[«,(«,  -A,„)-P  +  H']  =  0  (2-75) 


where  T  stands  for  temperature  and  the  unknown  variable  is  temperature.  The  slop  of  the 
secant  passing  through  two  points  is 


r  ft.  '  F(T, )  -  F(T„_, )  F(T„, )  -  F(T„ )  _  -  F(r, ) 

\  n  /  rrr  rj-t  rri  rj-t  m  rfi  '  ^ 


T  -T 


n- 1 


T  -  T 

1  n+l  1  - 


T  -T 

In+ 1  1  ” 


where  subscripts  n  and  n+l  represent  the  iteration  steps  and  F(Tn+])  is  the  same  as  F{a). 
a  is  the  solution  of  the  above  equation  and  F(Tn+ 0  becomes  zero  by  the  definition  of  the 
function.  Hence,  the  temperature  at  n+l  iteration  step  is 


1  n+l 


(2-77) 


This  calculation  procedure  is  repeated  until  the  convergence  tolerance  is  satisfied  as 
follows: 


(2-78) 


2.2.3.4.4  Expansion  Process 
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All  the  solenoid  valves  close  during  the  expansion  process  so  that  mass  is 
constant  but  volume  increases  in  the  chamber  as  shown  in  Figure  2.32.  The  pressure  and 
temperature  are  changed  due  to  the  volume  increase  and  some  energy  losses.  As  done  in 
the  intake  process,  the  thermodynamic  properties  are  predicted  by  considering  power  loss 
as  well  as  volume  increase.  No  heat  transfer  is  assumed  (adiabatic  flow)  but  the  power 
loss  is  considered  as  heat  loss.  The  pressure  is  assumed  to  be  constant  in  the  chamber  at 
each  time  step  so  that  the  work  done  by  CO2  can  be  calculated  by  using  Equation  (2-67). 

The  1st  law  of  thermodynamics  in  this  case  is 

^L  =  Q+W.  (2-79) 

dt 

By  multiplying  both  of  the  LHS  and  RHS  of  the  equation  by  dt,  the  above  equation  can 
be  written  as: 

m(u2  —  u, )  =  -P  +  W  (2-80) 

where  m  is  the  C02  mass  in  the  chamber.  Unknown  variable  is  only  u2  and  it  can  be 


expressed  as: 


w2  =  w,  + 


-P  +  W 
m 


(2-81) 


The  density  of  C02  at  state  ‘2’  becomes 


(2-82) 


The  pressure  and  temperature  (in  the  case  of  single-phase)  or  pressure  and  quality  (in  the 
case  of  two-phase)  at  state  ‘2’  are  determined  with  the  specific  internal  energy  and 
density.  The  bisection  method  is  used  to  find  them. 
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2.2.3. 4.5  Exhaust  Process 

The  exhaust  valve  opens  during  the  exhaust  period  and  mass  of  CO2  is  discharged 
to  the  inlet  to  the  evaporator  as  shown  in  Figure  2.33.  The  state  of  CO2  is  two-phase  in 
this  process  and  the  mass  flow  rate  is  calculated  by  using  the  method  in  Section  2.2. 3. 2.4. 
The  assumptions  are  the  same  as  in  the  intake  process.  The  kinetic  and  potential  energies 
are  neglected.  Even  though  the  process  is  an  adiabatic  flow,  the  power  loss  is  considered 
as  heat  loss.  The  pressure  is  assumed  to  be  constant  at  each  time  step  so  that  the  work 
done  by  CO2  is  calculated  by  Equation  (2-67).  Additional  assumptions  are  that  the  flow 
is  quasi-steady-state  and  homogeneous. 

There  are  two  steps  for  the  calculation  of  mass  flow  rate  and  thermodynamic 
properties  at  each  time  step  like  in  the  intake  process.  In  the  first  step,  the  mass  amount 
flowed  out  from  the  chamber  is  calculated  with  the  pressure  difference  between  the 
chamber  and  downstream  of  the  exhaust  valve.  There  is  no  piston  movement  at  this  time, 
i.e.,  V,  =  constant .  The  amount  of  CO2  flowed  out  from  the  chamber  becomes 

mou,=™ou<  'At-  (2-83) 

Then,  the  piston  moves  up  and  the  volume  of  the  chamber  decreases  with  the  new  mass 
of  CO2  in  the  second  step. 

The  thermodynamic  properties  are  calculated  with  new  specific  internal  energy 
and  density  by  using  the  bisection  method.  The  1st  law  of  thermodynamics  in  this  case  is 

dF 

—rr~  =  Q  +  W  -  ihouthOM .  (2-84) 

at 

After  some  mathematical  manipulations  with  thermodynamic  definitions,  the  above 
equation  can  be  written  as: 
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„2  =  — (m,B,  -  P  +  W  -  rnmhm )  (2-85) 

m2 

and  the  density  becomes  p2=  —  .  Then  the  pressure  and  quality  of  the  C02  are  found 

^2 

by  using  the  bisection  method. 

2.2.3.4.6  Total  Power  Loss  and  Net  Work  Output 

Total  power  loss  for  one  revolution  or  during  unit  time  scale  can  be  obtained  by 

summing  the  power  loss  at  each  time  step  as  follow: 

P--ZP.  <2-86) 

Net  work  output  is  also  obtained  by  adding  the  work  at  each  time  step  as. 

W„=  £w  +  2>+  (2-87) 

intake  expansion  exhaust 

During  the  intake  and  expansion  processes,  the  work  is  generated  by  C02  but  during  the 
exhaust  process,  the  work  is  needed  to  sweep  out  the  C02  from  the  chamber  to  the  inlet  to 
the  evaporator. 


2.3  Experiments  for  Cycle  with  ED-WOW  Device 
To  evaluate  the  effects  on  system  performance  the  prototype  expansion  device 
was  installed  in  a  working  transcritical  C02  refrigeration  cycle.  Figure  2.34  shows  a 
schematic  of  the  C02  based  environment  control  unit  (ECU)  of  a  unit  designed  for  the 
United  States  Army  [Li  and  Groll,  2001].  This  system  was  built  to  test  a  carbon  dioxide 
based  air-to-air  breadboard  system.  The  ED-WOW  will  be  incorporated  into  this 
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experimental  setup  after  modification.  Figure  2.35  illustrates  the  schematic  of  the  test 
setup  with  the  expansion  device.  As  shown  in  this  figure,  the  ED-WOW  does  not  replace 
the  expansion  valve  entirely.  Depending  on  the  final  capacity  of  the  prototype  expander, 
a  portion  of  total  working  fluid  may  be  bypassed  through  a  normal  expansion  valve. 


2.3.1  Description  of  Testing  Procedure 

Carbon  dioxide  is  pressurized  by  a  compressor,  achieving  a  supercritical  state  at 
the  outlet  of  the  compressor.  Heat  is  removed  through  the  gas  cooler  to  the  atmosphere 
and  the  thermodynamic  properties  of  carbon  dioxide  change  from  ‘vapor  like’  to  ‘liquid 
like’  fluid.  Then,  C02  expands  through  the  ED-WOW  from  high  pressure  to  low 
pressure  and  work  is  produced  during  this  process.  The  bypass  valve  allows  a  portion  of 
the  flow  to  pass  through  an  ordinary  expansion  valve.  The  carbon  dioxide  becomes  two- 
phase  at  the  outlet  of  the  ED-WOW.  Heat  is  then  transferred  to  the  C02  through  the 
evaporator  and  it  becomes  saturated  or  superheated  vapor  at  the  outlet  of  the  evaporator. 

The  system  of  Figure  2.35  is  installed  in  two  temperature  and  humidity  controlled 
psychrometric  rooms  so  that  independent  control  of  the  indoor  and  outdoor  environment 
can  be  achieved.  The  gas  cooler  is  located  in  the  outdoor  psychrometric  room  and  the 
other  parts  of  the  system  in  the  indoor  psychrometric  room. 

To  provide  a  load  to  absorb  extracted  work,  a  hydraulic  pump  with  water  as  the 
working  fluid  is  connected  to  the  output  shaft  of  the  ED-WOW  device.  The 
revolutionary  speed  of  the  ED-WOW  is  controlled  by  increasing  or  decreasing  the  load  of 
the  hydraulic  pump,  and  the  work  produced  by  the  ED-WOW  is  measured  indirectly  by 
measuring  the  performance  of  the  hydraulic  pump.  Figure  2.36  shows  the  schematic  of 
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the  connection  of  the  hydraulic  pump  system  to  the  ED-WOW.  Static  pressures  are 
measured  before  and  after  the  pump,  and  the  volumetric  flow  rate  of  the  water  (the 
working  fluid)  is  measured  using  a  rotameter.  The  ideal  pump  work  is  then  calculated 
using  the  pressure  difference  across  the  pump  and  the  water  flow  rate,  as  well  as  the 
pump  efficiency.  As  discussed  above,  this  hydraulic  pump  work  will  be  used  to  measure 
the  work  produced  through  the  ED-WOW  with  the  efficiency  of  the  pump. 
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Table  2.1:  Operating  conditions  for  a  typical  expansion  device  in  the  transcritical  CO2 


cycle. 


Parameter 

Operating  value 

High-side  pressure,  pi,igh 

10.2MPa 

High-side  temperature,  Thigh 

49°C 

Low-side  pressure,  piow 

3.4MPa 

Low-side  temperature,  Tiow 

Saturation  temperature  at  3.4MPa. 

CO2  flow  rate 

0.04kg/s 

Table  2.2:  Loss  coefficients  for  gradual  contractions:  Round  and  rectangular  ducts 

[Fox  and  McDonald,  1992]. 


A2 

Ai 

Included  angle,  0,  degree 

10 

15-40 

50-60 

90 

120 

150 

180 

0.50 

0.05 

0.05 

0.06 

0.12 

0.18 

0.24 

0.26 

0.25 

0.05 

0.04 

0.07 

0.17 

0.27 

0.35 

0.41 

0.10 

0.05 

0.05 

0.08 

0.19 

0.29 

0.37 

0.43 

Note:  Coefficients  are  based  on  hlm  =  K — .  Subscript  2  stands  for  small  pipe. 
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Table  2.3:  Relation  between  the  coefficient  of  contraction  Cc  and  area  ratio  a 


[Collier,  1972; 


1 

a 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

Cc 

0.586 

0.598 

0.625 

0.686 

0.790 

1.0 

U  J 

0.5 

0.45 

0.36 

0.21 

0.07 

0 

Table  2.4:  Re 

levant  engine  parameters. 

Parameters  used  by  Jeng,  and 

Radakovic  and  Khonsari 

Values 

Unit 

Engine  speed 

S  =  2000 

rpm 

Cylinder  bore 

B  =  88.9 

mm 

Bore  radius 

r  =  44.45 

mm 

Crank  radius 

R  =  40.0 

mm 

Connecting  rod  length 

L=  141.9 

mm 

Composite  roughness 

<7=  0.5 

pm 

Piston  ring  width 

b  =  1.475 

mm 

Ring  thickness 

tr  =  3.8 

mm 

Ring  modulus 

E  =  70.0 

GN/rrr 

Ring  crown  height 

CR=  14.9 

pm 

Friction  coefficient 

/=  0.08 

“ 

Lubricant  viscosity 

H  =  0.00689 

Pa-s 

Tangential  ring  tension 

T=  22.38 

N 

Figure  2.1:  Basic  transcritical  CO2  cycle  (Base  cycle)  -  Case  1. 


Evapaorator 

Figure  2.2:  Base  cycle  with  ED-WOW  -  Case  2. 


Evapaorator 

Figure  2.3:  Base  cycle  with  internal  heat  exchanger  (IHX  cycle)  -  Case  3 


Figure  2.4:  IHX  cycle  with  ED-WOW  -  Case  4. 
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Evapaorator 

Figure  2.5:  Intercooler  cycle  -  Case  5  and  Case  8. 


Evapaorator 

Figure  2.6:  Intercooler  cycle  with  ED-WOW  -  Case  6  and  Case  9. 
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1st  stage 
compressor 


2nd  stage 
compressor 


Evapaorator 

Figure  2.7:  Economizer  cycle  -  Case  7  and  Case  10. 


Figure  2.8:  Piston-cylinder  device  with  the  related  components  and  piston  positions. 
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Volume  (cm3) 


Figure  2.9:  Piston  processes  with  assumption  of  ideal  process. 
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Figure  2.10:  Schematic  of  connection  of  the  solenoid  valve  to  the  fitting  and  cylinder 

chamber. 


coefficients  for  flow  through  sudden  area  changes  [Fox 
McDonald,  1992]. 


Expansion  loss  coefficient,  Kt 


Reentrant: 

K=0.78 


Sharp-edged: 

K=0.4-0.5 


Slightly  rounded: 
K=0.2-0.25 


Well-rounded: 

K=0.05 


Figure  2.13:  Entrance  loss  coefficients  [Evett  and  Liu,  1987], 


Figure  2.14:  Exit  loss. 
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Fully  developed  inlet  flow 


Thin  inlet  boundary  layer 


Total  included  diffuser  angle  20,  degrees 


Figure  2.15:  Flow  losses  in  a  gradual  conical  expansion  region  [Evett  and  Liu,  1987] 
Note  that  the  loss  is  based  on  velocity  head  in  the  small  pipe. 


%  20 


Deflection  angle,  $  (degrees) 

Figure  2.16:  Representative  total  resistance  (L</D)  for  miter  bends 
[Fox  and  McDonald,  1992]. 
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Figure  2.19:  Flow  diameter  of  the  computer  program  modeling. 


(a)  Sudden  contraction 


(b)  Miter  bend 


(c)  Sudden  enlargement  (d)  Exit  loss 


Figure  2.20:  Pressure  losses. 
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Crank  Angla  (dagraaa) 


Figure  2.28:  Power  losses  predicted  by  model  of  Radkovic  and  Khonsari  [Radkovic  and 

Khonsari,  1997]. 


Figure  2.29:  Minimum  film  thickness  versus  crank  angle. 


Legend  of  sensors  :  T-temperature,  P-pressure,  V-volumetric  flow  rate,  M-mass  flow  rate 
Db-dry  bulb  temperature,  Wb-wet  bulb  temperature,  W-power, 
DP-differential  pressure. 


Figure  2.34:  Schematic  of  CO2  based  ECU  test  setup. 


Legend  of  sensors  :  T-temperature,  P-pressure,  M-mass  flow  rate,  W-power, 
DP-differential  pressure. 

Figure  2.35:  Schematic  of  experimental  setup  for  transcritical  C02  cycle  with  ED-WOW 


T ranscritical  C02  cycle 


Figure  2.30: 


Schematic  of  the  connection  of  the  hydraulic  pump  system  to  the  expansion 
device  producing  work. 
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3  DESIGN  OF  ED-WOW  EXPANSION  DEVICE 

As  discussed  in  Chapter  2,  a  piston-cylinder  assembly  device  was  chosen  as  the 
appropriate  expansion  device  for  the  prototype  application.  A  representative  expansion 
device  was  designed  and  fabricated.  The  parameters  for  a  typical  expansion  device  for 
prototype  transcritical  CO2  cycle  are  given  in  Table  2.1. 

The  details  of  the  design  of  the  expansion  device  with  output  of  work  (ED-WOW) 
are  given  in  this  chapter.  The  basic  design  is  based  on  a  piston  expander  controlled  by 
electrically  actuated  valves.  The  crankshaft  of  the  piston  expander  is  connected  to  a 
mechanical  loading  device,  a  hydraulic  pump,  which  absorbs  the  work  produced  during 
the  expansion  of  C02. 


3.1  Design  of  Control  Valves 

The  expansion  process  is  controlled  using  fast-acting  solenoid  valves  as  intake 
and  exhaust  valves.  As  smaller  valves  exhibit  more  rapid  cycling  response,  the  valves 
must  be  sized  for  the  specific  application.  Two  valve  sizes,  indicated  by  orifice  diameter, 
are  used  for  the  solenoid  valves.  Solenoid  valves  with  3/64  in  orifice  diameter  are  used 
as  the  intake  valve.  Solenoid  valves  with  1/16  in  orifice  diameter  are  the  exhaust  valves. 
There  is  some  time  lag  for  the  solenoid  valves  to  response  to  an  input  signal.  In  order  to 
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control  the  expansion  process  as  desired,  the  characteristics  of  the  valve  should  be 
known. 


3.1.1  Experiments  for  the  Characteristics  of  the  Solenoid  Valve 

Experiments  have  been  performed  to  characterize  the  filling  behavior  of  C02  into 
the  ED-WOW.  For  this  purpose,  a  fixed-volume  chamber  was  designed  that  represents 
the  volume  of  the  ED-WOW  at  the  instant  of  time  on  the  expansion  stroke  when  the 
intake  valve  is  required  to  close.  In  order  to  guarantee  minimum  throttling  in  the 
expansion  device,  the  time  that  it  takes  for  the  C02  to  fill  the  chamber  volume  and  to 
pressurize  the  chamber  up  to  the  pressure  upstream  of  the  intake  valve  needs  to  be 
known.  Figure  3.1  shows  the  schematic  of  the  experimental  test  stand.  The  electric 
heater  is  used  to  increase  the  C02  pressure  inside  the  small  tank  up  to  the  desired 
conditions. 

The  solenoid  valve  with  3/64  in  orifice  diameter  was  tested.  The  pressure 
transducer  installed  in  the  chamber  is  used  to  measure  the  dynamic  pressure  inside  the 
chamber  after  energizing  the  solenoid  valve  (opening  of  the  valve).  This  transducer  is 
AC  coupled,  and  thus  only  provides  the  unsteady  portion  of  the  pressure  signal.  Thus, 
the  peak  of  the  unsteady  signal  indicates  when  the  volume  can  be  considered  ‘filled’. 
Figures  3.2  and  3.3  show  the  pressure  signal  (in  Volts)  as  C02  fills  the  chamber.  In  the 
legend  of  the  figures,  p\  and  p2  stand  for  the  initial  pressures  upstream  of  the  valve  and  in 
the  chamber,  respectively.  As  shown  in  these  figures,  the  pressure  starts  to  increase 
approximately  0.01  seconds  after  opening  of  the  valve.  This  time  is  the  time  lag  when 
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the  intake  valve  opens.  This  time  will  be  taken  as  that  representative  of  the  time  required 
to  open  the  intake  valves. 

It  can  be  seen  from  Figure  3.2  that  the  time  that  it  takes  to  reach  the  maximum 
pressure  increases  slightly  as  the  inlet  pressure  increases.  Figure  3.3  shows  the 
comparison  of  the  pressure  increase  for  similar  absolute  pressure  differences  across  the 
solenoid  valve,  but  with  different  inlet  pressures.  From  Figure  3.3  it  can  be  seen  that  the 
time  that  it  takes  to  reach  the  maximum  pressure  decreases  as  the  inlet  pressure  increases. 
As  shown  in  Figures  3.2  and  3.3,  it  takes  more  than  0.075  seconds  for  the  pressures  to 
reach  maximum  values.  The  design  calculation  shows  that  the  possible  valve  opening 
time  is  approximately  0.072  seconds  for  the  revolution  speed  of  the  expansion  device  of 
120  rpm  (See  Section  3.7).  Therefore,  it  was  determined  that  two  intake  valves  and  one 
exhaust  valve  must  be  used  in  order  to  give  large  enough  flow  area  into  the  piston- 
cylinder  chamber. 

Other  experiments  have  been  conducted  without  CO2  flow,  i.e.,  no  pressure 
difference  across  the  valve  to  measure  the  valve  closing  time.  Table  3.1  shows  the 
response  time  to  the  signal  to  close  the  solenoid  valve  with  3/64  in  orifice  diameter.  The 
average  time  is  0.05323  s  and  its  95%  confidence  interval  is  0.05253-0.05393.  The 
average  time  will  be  taken  as  that  representative  of  the  time  required  to  close  the  intake 
and  exhaust  valves. 

Table  3.1  also  shows  the  response  time  to  the  signal  to  open  the  valve  when  there 
is  no  pressure  difference  across  the  valve.  Its  average  time  and  95%  confidence  interval 
are  0.0104  s  and  0.009994-0.010806  s,  respectively.  This  average  time  will  be  taken  as 
that  representative  of  the  time  required  to  open  the  exhaust  valve. 


91 


3.1.2  Example  of  Valve  Control 

Figure  3.4  shows  the  desired  piston  displacement  from  the  top  ceiling  of  the 
piston-cylinder  device  over  time  as  determined  from  the  expansion  device  producing 
work  design.  The  assumed  revolutionary  speed  is  120  rpm  and  the  lengths  of  the 
connecting  rod  and  crank  are  38.8  mm  and  11.0  mm.  The  detailed  description  about  the 
ED-WOW  is  given  in  the  next  sections.  Also  indicated  in  Figure  3.4  are  the  timings  for 
the  intake  and  exhaust  valves.  As  discussed  before,  the  design  calculation  indicates  that 
the  intake  valve  should  open  at  the  TDC  at  time  zero  and  close  when  the  piston  reaches 
5.21  mm  from  TDC  (at  the  displacement  of  8.01  mm)  at  0.072  seconds.  The  exhaust 
valve  should  open  at  the  displacement  of  24.80  mm  (Bottom  Dead  Center,  BDC)  at  0.25 
seconds  and  close  when  the  piston  reaches  the  TDC  again.  The  lead  times  for  the  control 
signals  required  for  the  valves  to  open  and  close  fully  are  indicated  in  the  figure.  The 
circle  and  square  stand  for  the  timings  of  input  signal  to  open  and  close  the  intake  and 

exhaust  valves  in  Figure  3.4,  respectively. 

An  encoder  is  attached  to  the  ED-WOW  to  read  the  exact  angle  of  the  crankshaft 
of  the  device.  The  encoder  generates  a  gate  signal  (Z  signal)  and  an  A  signal  of  2048 
pulses  per  revolution.  This  encoder  is  connected  to  a  data  acquisition  system  and  a 
program  was  written  to  control  the  solenoid  valves  based  on  the  reading  of  the  angle  of 
the  crankshaft. 


3.2  Design  of  Exnander  Cylinder,  Piston,  and  Crank 
The  prototype  piston  type  ED-WOW  is  based  on  a  highly  modified  small  four¬ 
cycle,  two-piston  engine  of  displacement  13.26  cc  x  2  that  is  commercially  available.  In 
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the  original  design,  both  pistons  are  on  the  same  firing  order,  meaning  that  both  pistons 
achieve  TDC  at  the  same  instant.  In  the  CO2  expander  application,  the  goal  is  to  reduce 
the  need  for  mechanical  inertia  (e.g.  a  flywheel)  by  using  an  out-of-phase  firing  order, 
with  the  piston  on  the  expansion  stroke  providing  the  force  required  to  drive  the  opposing 
piston  through  the  evacuation  stroke.  To  achieve  this,  the  crankshaft  of  the  original 
device  was  modified,  and  is  shown  in  Figure  3.5. 

An  assembly  of  the  piston,  piston  ring,  and  connecting  rod  is  used  for  the 
expansion  device.  Figure  3.6  shows  the  assembly  that  is  commercially  available.  The 
assembly  is  used  to  deliver  the  force  of  CO2  inside  the  cylinder  chamber  and  rotate  the 
crank  of  the  device  described  above.  Table  3.2  shows  the  dimensions  of  the  cylinder, 
connecting  rod,  crank,  etc. 


3.3  Design  of  Top  of  Expander  Cylinder 

The  cylinder  head  connects  the  dual  intake  and  single  exhaust  solenoid  valves  to 
the  cylinder.  The  major  design  constraint  for  the  device  was  achieving  the  absolute 
minimum  dead  volume  in  the  system  (i.e.  the  volume  downstream  of  the  intake  valve  at 
piston  top  dead  center).  This  minimal  dead  volume  is  critical  to  achieving  an  efficient 
expansion  process,  as  it  is  directly  proportional  to  the  throttling  losses  at  the  beginning  of 
the  expansion  stroke.  Figure  3.7  presents  several  views  of  the  cylinder  head  design, 
showing  the  three  ports  that  allow  for  connection  of  the  valves.  Each  of  the  two  cylinder 
heads  is  bolted  to  the  crank  housing  with  four  machine  screws.  The  three  valves  for  each 
head  connect  via  1/8  in  NPT  fittings,  and  in  turn  are  attached  to  a  retention  plate  (not 
shown)  with  machine  screws.  Tie  rods  span  between  the  cylinder  head  retention  plates 
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on  each  of  the  two  opposed  cylinders,  allowing  for  a  rigid  structure  resistant  to  vibration- 
induced  leakage  in  the  head  gasket.  Figures  3.8  and  3.9  illustrate  the  dimensions  of  the 
head  (in  inches). 


3.4  Design  of  Enclosure 

The  prototype  ED- WOW  has  been  designed  to  operate  at  an  inlet  pressure  of 
1500  psig  and  an  outlet  pressure  of  500  psig.  Due  to  these  high  pressures,  it  is 
anticipated  that  some  C02  will  leak  from  the  expansion  volume  to  the  crank  housing 
along  the  piston  rings,  and  from  the  crank  housing  to  the  surrounding  environment 
through  bearings  and  seals.  Therefore,  the  whole  ED-WOW  is  placed  in  an  enclosure  to 
prevent  the  refrigerant  from  leaking  out  of  the  system.  The  enclosure  is  under  discharge 
pressure  as  it  is  intended  that  the  C02  exiting  the  discharge  valve  enters  the  enclosure. 
Thus,  the  enclosure  needs  to  be  strong  enough  to  withstand  the  discharge  pressure  of 
about  500  psig  during  normal  operation. 

Based  on  current  design  of  the  ED-WOW,  the  inner  diameter  of  the  enclosure 
needed  to  be  at  least  16  inches  to  fit  all  valves  and  mountings  into  the  enclosure.  Figure 
3.10  shows  a  cut  view  of  the  enclosure  with  the  expansion  device  drawn  to  scale.  As 
shown  in  this  figure,  the  expansion  device  can  be  located  in  two  domed  pipe  end  caps 
with  an  18  inch  outer  diameter.  Flanges  are  welded  to  each  domed  end  cap  to  allow  the 
opening  and  closing  of  the  enclosure  using  bolted  joints.  The  material  of  the  domed  end 
caps  is  schedule  10,  304-stainless  steel,  the  inner  diameter  is  17.25  inch,  and  the  wall 
thickness  is  3/8  inch.  Figure  3.11  shows  several  views  of  the  enclosure.  The  flanges  of 
the  enclosure  were  sealed  using  an  O-ring. 
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The  crankshaft  of  the  ED-WOW  must  penetrate  through  the  enclosure  so  that  it 
can  be  connected  to  the  loading  device  which  absorbs  the  extracted  energy  when  the  CO2 
expands  from  high  to  low  pressure.  A  penetration  is  made  at  the  center  of  one  of  the 
domed  end  caps  for  the  crankshaft.  In  order  to  prevent  CO2  leakage  through  the  shaft 
penetration,  a  high-pressure  rotary  shaft  seal  was  installed  in  the  penetration. 

The  enclosure  has  two  additional  penetrations:  one  for  the  high-pressure  inlet 
piping  to  the  expansion  device  and  the  other  one  for  the  extraction  of  the  expanded  CO2 
at  the  bottom  of  the  enclosure  which  is  fed  to  the  inlet  of  the  evaporator.  In  addition,  a 
safety  relief  valve  was  installed  at  the  top  of  the  enclosure. 

Several  wires  needed  to  be  connected  to  the  encoder  and  the  solenoid  valves 
inside  the  enclosure  in  order  to  read  the  angle  of  the  crankshaft  and  to  open  and  close  the 
valves  to  control  the  CO2  expansion  process.  A  commercially  available  “electrical 
feedthrough”  was  plugged  into  the  wall  of  one  of  the  domed  end  caps.  Since  the 
electrical  feedthrough  is  designed  for  high  pressure  operation,  no  CO2  leakage  is 
expected.  The  outside  and  inside  wires  are  connected  on  both  sides  of  the  electrical 
feedthrough. 

A  flat  plate  is  welded  to  the  inside  of  the  enclosure  to  provide  the  base  for  the 
expansion  device.  Figure  3.12  shows  the  fabricated  enclosure.  Figures  3.12  (a)  and  (b) 
illustrate  the  front  half  of  the  enclosure  with  the  ED-WOW  mounted  on  the  plate  and  the 
rear  half  of  the  enclosure,  respectively.  Table  3.3  shows  the  summary  of  the  design  of 


the  enclosure. 
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3.5  Design  of  Loading  Device 

A  mechanical  loading  device  is  connected  to  the  ED-WOW.  A  hydraulic  pump  is 
used  as  the  mechanical  loading  device.  Figure  3.13  shows  the  ED-WOW  with  the 
hydraulic  pump  before  the  expansion  device  is  placed  inside  the  enclosure.  The  load  is 
increased  by  closing  the  hand  valves  at  the  inlet  to  the  hydraulic  pump  and  thus, 
increasing  the  flow  resistance  to  the  fluid  in  the  hydraulic  pump  system.  There  are  two 
kinds  of  valves,  a  ball  valve  and  a  needle  valve.  The  ball  valve  is  used  to  adjust  the 
coarse  load  amount  and  the  needle  valve  is  for  the  fine  adjustment. 


3.5.1  Calibration  of  Hydraulic  Pump 

The  hydraulic  pump  has  been  calibrated  to  calculate  its  efficiency.  This  pump 
efficiency  is  used  to  measure  the  work  output  through  the  expansion  device.  During 
normal  operation  of  the  transcritical  CO2  system,  the  power  is  measured  with  the 
hydraulic  pump.  Then,  the  extracted  work  through  the  ED-WOW  is  calculated  using  the 
measured  power  and  pump  efficiency. 

Figure  3.14  shows  the  schematic  of  the  set-up  for  the  measurement  of  torque 
using  a  spring  scale  and  string.  A  string  is  tied  up  to  a  spring  scale,  wound  on  the  large 
rod  several  times  and  a  mass  is  given  at  the  end  of  the  other  side  of  the  string,  as  shown 
in  Figure  3.14.  As  the  crankshaft  of  the  expansion  device  rotates,  the  string  becomes  taut 
due  to  the  friction  between  the  surface  of  the  rod  and  the  string.  At  a  certain  point  the 
system  reaches  an  equilibrium  state  owing  to  the  torque  balance  between  the  ED-WOW 
and  the  system  of  the  spring  scale,  string,  and  mass.  At  that  point  the  rotational  speed  of 
the  expansion  device  becomes  constant  and  the  string  delivers  forces  to  the  spring  scale 


96 


and  mass.  The  spring  force  and  rotation  speed  of  the  ED- WOW  are  measured.  The 

radius  of  the  large  rod  is  4.37  cm.  The  torque  is  calculated  as: 

T,=(Fsa,„+Fm)xr  0-1) 

where  Tq,  Fscaie,  F,mss,  and  r  are  the  torque,  the  force  measured  by  spring  scale,  the  force 
due  to  the  mass,  and  radius  of  the  rod.  Then,  the  power  becomes 

P„=Ttxa>  (3-2) 

Pss  and  co  stands  for  the  power  calculated  using  the  torque  of  the  ED-WOW  and  the 
rotational  speed  of  the  expansion  device  in  radian/s,  respectively. 

Then,  the  pump  power  is  calculated.  The  pump  connection  to  the  ED-WOW  is 
shown  in  Figure  3.15.  The  pressures  are  measured  before  and  after  the  hydraulic  pump 
and  volumetric  water  flow  rate  is  also  measured  using  a  rotameter.  The  power  of  the 
hydraulic  pump  is  calculated  as  follows. 

^=Apx(2  (3-3) 

where  Ppp,  Ap,  and  Q  are  the  power  of  the  hydraulic  pump,  pressure  difference  across  the 

hydraulic  pump,  and  the  volumetric  water  flow  rate,  respectively. 

The  hydraulic  pump  efficiency  can  be  calculated  by  comparing  Pss  and  Ppp  for  the 
similar  test  conditions.  Table  3.4  and  Table  3.5  show  the  measured  values  of  the  forces 
and  pressures,  the  calculated  torque  and  powers,  and  the  pump  efficiencies.  The 
description  for  the  parameters  is  given  in  the  tables.  As  shown  in  the  tables,  the  pump 
efficiency  becomes  68.14%  when  p;n,  the  inlet  pressure  of  the  ED-WOW,  is  about  4.3 
MPa  and  the  rotational  speed  is  about  300  RPM.  It  also  becomes  47.48%  when  pjn  and  CD 


97 


are  about  3.0  MPa  and  294  RPM.  As  discussed  before,  this  efficiency  of  the  hydraulic 
pump  is  used  to  measure  the  work  produced  through  the  ED- WOW. 

The  pump  was  calibrated  at  two  points  only  due  to  the  difficulty  of  the 
experiment.  The  initial  desire  of  this  calibration  was  to  obtain  enough  data  and  to  make  a 
calibration  curve  for  the  pump  efficiency.  However,  the  system  was  too  sensitive  to  the 
friction  of  the  string.  The  ED-WOW  did  not  work  once  the  friction  was  further  increased 
to  decrease  the  rotational  speed  of  the  ED-WOW. 


3.6  Establishment  of  Input  Timings  for  the  Intake  and  Exhaust  Valves 
Another  experiment  has  been  performed  to  establish  correct  input  timings  for  the 
solenoid  valves  to  control  the  expansion  process  of  C02  through  the  expansion  device 
producing  work.  Figure  3.16  shows  the  schematic  of  the  test  stand. 

As  mentioned  in  the  previous  sections,  the  input  signals  should  be  given  prior  to 
certain  locations  of  the  crankshaft  due  to  the  time  lags  of  the  solenoid  valves.  However, 
the  expansion  device  was  not  working  when  the  originally  predicted  time  lags  were 
considered.  It  seemed  that  immediately  after  the  signal  to  open  the  intake  valve  #1  is 
given,  high-pressure  gas  flows  into  cylinder  #1  even  though  piston  #1  has  not  reached 
Top  Dead  Center  (TDC).  At  the  same  time,  the  high-pressure  gas  inside  cylinder  #2 
pushes  piston  #2  as  it  expands  from  the  high  pressure  to  the  low  pressure.  Number  1  (#1) 
and  number  2  (#2)  illustrate  the  one  side  of  two  piston-cylinders  and  the  other  side, 
respectively.  If  the  pressure  in  cylinder  #1  increases  rapidly  and  the  device  has  not 
enough  momentum  to  move  piston  #1  through  the  TDC  position  (against  the  force  of  the 
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incoming  high-pressure  gas,  since  the  intake  valve  opened  too  early),  piston  #1  will  move 
backwards  until  the  pressures  in  both  cylinders  have  equalized  and  the  device  is  stuck. 

Thus,  the  timing  table  shown  in  Figure  3.4  has  been  modified.  The  signal  to  open 
the  intake  solenoid  valve  was  delayed  until  just  after  the  piston  reached  TDC.  However, 
the  signal  to  close  the  intake  valve  remained  the  same  in  order  to  achieve  the  desired 
volume  expansion  ratio.  Figures  3.17  and  3.18  show  the  new  timing  tables  for  the  intake 
valves  #1  and  the  exhaust  valve  #2,  and  for  the  intake  valves  #2  and  the  exhaust  valve  #1, 
respectively.  Figure  3.19  shows  the  timing  table  for  all  solenoid  valves.  The  x-axis 
stands  for  the  number  of  pulses  of  the  encoder  connected  to  the  expansion  device  for  one 
revolution  of  the  crankshaft  (2048  pulses/revolution),  and  the  y-axis  stands  for  on  (=1, 
open)  and  off  (=0,  closed)  of  the  valve  actuation.  As  shown  in  the  figures,  the  signal  to 
open  the  intake  valves  is  given  10  pulses  after  Top  Dead  Center  (TDC)  and  the  valves 
remain  open  for  70  pulses.  The  signal  to  open  the  exhaust  valve  is  given  10  pulses  before 
BDC  and  the  signal  to  close  the  exhaust  valve  is  given  10  pulses  after  TDC.  This  timing 
table  is  based  on  a  crankshaft  speed  of  120  rpm. 

Appendix  C  shows  the  program  written  in  ‘Visual  C++’  to  control  the  solenoid 
valves.  As  mentioned  above,  this  program  is  written  based  on  the  rotational  speed  of  120 
RPM. 


3.7  Determination  of  Revolution  Speed  of  the  Expansion  Device 
As  mentioned  above,  the  revolution  speed  of  the  ED-WOW  is  determined  as  120 
RPM.  Based  on  the  desired  operating  condition  shown  in  Table  2.1  and  the  dimensions 
of  the  ED-WOW  shown  in  Table  3.2,  the  LOE  (location  of  expansion  initiation)  must  be 
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5.21  mm  from  the  location  of  the  TDC.  Table  3.6  shows  the  revolution  speed  of  the 
expansion  device  versus  the  time  that  it  takes  for  the  piston  to  reach  5.21  mm  from  the 

TDC  during  the  down  stroke  (intake  stroke). 

As  mentioned  above,  the  time  lag  of  the  intake  valve  to  the  input  signal  for 
closing  is  0.053s.  The  input  signal  must  be  given  0.053s  prior  to  the  piston  reaching  the 
LOE.  The  time  to  reach  the  LOE  should  be  larger  than  0.053s.  As  shown  in  Table  3.6, 
when  the  revolution  speed  is  equal  to  or  larger  than  180  rpm,  the  intake  valve  can  t  be 
closed  at  the  exact  position  of  the  piston  of  the  LOE  because  of  the  short  time  duration  of 
the  intake  process  from  the  TDC  to  the  LOE.  The  intake  valve  will  be  closed  after  the 
LOE.  In  this  case,  C02  will  be  under-expanded  due  to  the  larger  amount  of  mass  at  the 
initiation  of  the  expansion.  Therefore,  the  revolution  speed  should  be  slower  than  180 
rpm. 

When  the  speed  of  the  device  is  too  low,  the  momentum  becomes  too  low  to  keep 
running  the  device  during  normal  operation.  Therefore,  the  speed  of  120  rpm  is 
determined  rather  than  60  rpm  as  the  operating  speed  of  the  expansion.  In  this  case,  there 
is  enough  time  to  close  the  intake  valve  at  the  LOE  to  get  the  desired  low-side  pressure 
and  there  is  also  comparatively  large  momentum  to  run  the  ED-WOW  continuously. 
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Table  3.1:  Valve  fully  opening  and  closing  timings  for  the  solenoid  valve  with  3/64  in 

orifice. 


Case 

Time  (s) 

Average  time 
(s) 

Standard 
deviation  (s) 

95%  confidence 
interval  (s) 

Measurement 
of  closing  time 

0.05323 

0.000667 

0.05253 

0.05393 

Measurement 
of  opening  time 

0.0106 

0.0105 

0.0106 

0.0096 

0.0108 

0.0103 

0.0104 

0.000387 

0.009994 

0.010806 

Table  3.2:  Dimensions  of  the  expansion  device  producing  work. 


Diameter  of  cylinder 

27.7  mm 

Displacement 

13.26  cc 

Stroke 

22.0  mm 

Length  of  crank 

11.0  mm 

Length  of  connecting  rod 

38.75  mm 

Length  from  the  top  of  the  cylinder  to  the  head 
of  the  piston  at  TDC 

2.8  mm 

Volume  of  cylinder  at  TDC 

1.685  cm3 

Volume  of  cylinder  at  BTC 

14.95  cm3 

Volume  of  fittings  connected  to  the  device 

0.813  cm3 
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Table  3.3:  Summary  of  the  design  of  the  enclosure. 


Dimension 

1.  Outer  diameter:  18” 

2.  Inner  diameter:  17.624” 

3.  Wall  thickness:  3/16” 

4.  Total  length:  16” 

5.  Length  of  straight  section:  about  6.6” 

Material 

1.  Schedule  10  304-stainless  steel 

Penetrations 

1.  One  at  the  top  for  the  crankshaft  that  will  be  connected  to  the 
loading  device. 

2.  One  for  high-pressure  inlet  piping  to  the  expansion  device. 

3.  One  at  the  top  for  a  safety  relief  valve  (1/4”  male  NPT). 

4.  One  at  the  top  for  a  pressure  gage. 

5.  One  at  the  top  for  temperature  measurement. 

6.  One  at  the  bottom  of  the  enclosure  to  extract  the  expanded 
CO2  to  an  evaporator. 

7.  One  at  the  bottom  to  release  the  leaked  CO2  to  the  inlet  of  the 
evaporator  through  a  check  valve. 

8.  Two  at  the  top  for  the  electrical  feedthroughs  to  connect  the 
inside  and  outside  wires 

Welding 

1.  One  flat  plate  to  support  the  expansion  device  inside  the 
enclosure. 

2.  Two  flanges  to  fasten  the  domed  caps  with  some  bolts. 

Total  weight 

1.  About  235  lb. 

Table  3.4:  Power  measurement  using  torque  of  the  expansion  device  producing  work. 


Test  no. 

1 

V  * - 

2 

3 

4 

Pm  (kPa), 

STD 

4323.9, 

48.4 

2974.6, 

46.0 

2990.1, 

38.9 

3038.5, 

60.1 

Pout  (kPa), 

STD 

161.9, 

83.2 

137.6, 

43.6 

N/A 

N/A 

Tin  (°C), 

STD 

16.38, 

0.33 

16.78, 

0.42 

18.01, 

0.29 

17.43, 

0.39 

Ton,  (°C), 

STD 

3.40, 

0.72 

3.64, 

0.84 

12.27, 

0.72 

10.73, 

0.40 

F-scale  (kgf) 

6.03  ±0.23 

4.9010.23 

4.9010.23 

4.9010.23 

F-mass  (kgf) 

0.05 

0.15 

0.055 

0.055 

Total  F  (N) 

6.0810.23 

5.0510.23 

4.9610.23 

4.9610.23 

(0  (RPM) 

300112 

282112 

288112 

294112 

Torque  (N-m) 

2.61 

2.16 

2.12 

2.12 

Power  (W) 

Pss 

85.19 

66.63 

64.0 

64.71 

102 


Table  3.5:  Power  measurement  of  the  hydraulic  pump  and  its  efficiency 


Test  no. 

5 

6 

7 

8 

9 

Pin  (kPa), 

4436.4, 

2980.4, 

3055.1, 

3065.5, 

N/A 

STD 

36.7 

26.6 

64.0 

55.7 

Poul  (kPa), 
STD 

139.4, 

48.2 

130.1, 

34.2 

133.2, 

36.1 

125.4, 

38.9 

N/A 

Ti„  (°C), 

20.37, 

16.54, 

16.42, 

N/A 

STD 

0.27 

0.22 

mBRM' 

0.15 

Ton.  (°C), 

STD 

15.08, 

0.98 

8.21, 

0.20 

7.75, 

0.21 

7.80, 

0.21 

N/A 

PpJn  (Psig) 

-2.5  ±0.1 

-  2.1  ±0.1 

-2.1  ±0.1 

-2.1±0.1 

—  2. 1  ±  0. 1 

Pp.out  (psig) 

23.0  ±0.2 

12.0  ±0.2 

11.0  ±0.2 

10.0  ±0.2 

12.0  ±0.2 

Q  (GPM) 

3.3±0.1 

3.3±0.1 

3.3±0.1 

3.3±0.1 

3.3±0.1 

CO  (RPM) 

294  ±12 

294  ±12 

294  ±12 

294  ±12 

294  ±12 

Power  (W) 

p 

r  DP 

58.04 

32.10 

29.82 

27.54 

32.10 

Pump  Eff. 

(%) 

68.14 

50.15 

46.59 

43.04 

50.15 

Descriptions  of  parameters:  STD:  standard  deviation,  Pin:  inlet  pressure  of  the  ED- 
WOW,  Pout:  outlet  pressure  of  the  ED-WOW,  Tj„:  inlet  temperature  of  the  ED-WOW, 
Tout:  outlet  temperature  of  the  ED-WOW,  F-scale:  force  measured  by  the  scale,  F-mass: 
force  due  to  the  mass,  Pppn:  inlet  pressure  of  the  hydraulic  pump,  Pp,out:  inlet  pressure  of 
the  hydraulic  pump,  Q :  volumetric  flow  rate  of  water. 


Table  3.6:  Revolution  speed  of  the  expansion  device  versus  the  time  to  reach  LOE. 


Speed  (rpm) 

Angular  velocity  (rad/s) 

Time  (s) 

60 

6.28 

0.1442 

120 

12.57 

0.0721 

180 

18.85 

0.0481 

240 

25.13 

0.0360 

Voltage  (V) 


Figure  3.1: 


10: 


Pressure  Pressure 


Schematic  of  test  stand  to  determine  the  characteristics  of  the  solenoid  valve. 


Figure  3.2:  Filling  behaviors  of  C02  in  the  chamber  with  the  intake  valve 

(Dorifice==3/64  in). 
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Figure  3.3:  Comparison  of  filling  behaviors  for  similar  pressure  differences,  but  with 

different  p\ s. 


Figure  3.4:  Timings  for  input  signals  to  the  solenoid  valves. 
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Figure  3.5:  Crankshaft  of  the  expansion  device  producing  work. 


Figure  3.6: 


Assembly  of  a  piston,  piston  ring  and  connecting  rod. 


■  •  • 


6 


7 
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Figure  3.10:  Cut  view  of  the  designed  enclosure  with  expansion  device  in  scale. 


Figure  3.12:  Fabricated  enclosure. 

Front  half  of  the  enclosure  with  the  ED- WOW. 
(b):  Rear  half  of  the  enclosure. 
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device 


Motor 


Hydraulic 

pump 


Ball  valve 


Ball  valve 


Hydraulic 

pump 


Expansion 

device 


Motor 


Figure  3.13:  Front  and  side  views  of  the  expansion  device  and  hydraulic  pump 
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Figure  3.14:  Schematic  of  set-up  for  torque  measurement  using  spring  scale  and  string. 


On/Off  of  valves  (-)  C  On/Off  of  valves  (-) 


1.2 


Timing  table  for  the  1 81  intake  and  2nd  exhaust  valves. 


1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

-0.2 

-256  0  256  512  768  1024  1280  1536  1792  2048  2304 

Number  of  pulses  (-) 


- 1*  Intake  valve 

- 2nd  Exhaust  valve 


TDC  <1"  piston)  8DC  o'*  piston)  TDC  (1“  piston) 

BOC  (2“  piston)  TDC  <2"  Pte,on>  BDC  <2”  P^'0") 

J _ I _ l _ I _ l _ I _ l _ I _ l _ I _ l _ I _ l _ I _ l _ I _ i _ I _ l _ 


3.17:  Valve  on/off  timing  of  1st  intake  and  2nd  exhaust  valves. 


1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 


-0.2 


Timing  table  for  the  2nd  intake  and  1 81  exhaust  valves. 


- 1*  Exhaust  valve 

-  -  -  -  2“  Intake  valve 


1  piston) 

BDC  {2  piston) 

TDC  (2nd  piston) 

piston) 

TDC  (1*  piston) 

BDC  (1“  piston) 

i  .  i  .  i 


-256  0  256  512  768  1024  1280  1536  1792  2048  2304 


Number  of  pulses  (-) 


Figure  3.18:  Valve  on/off  timing  of  2nd  intake  and  1st  exhaust  valves 


On/Off  of  Valve  (-) 


Timing  table  for  the  solenoid  valves  to  control  the  expansion  process. 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

-0.2 

-256  0  256  512  768  1024  1280  1536  1792  2048  2304 

Number  of  pulses  (-) 


- 1  Intake  valve 
- 18*  Exhaust  valve 
2nd  Intake  Valve 
-  2nd  Exhaust  Valve 


Figure  3.19:  Desired  timing  table  to  control  all  solenoid  valves 
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4  RESULTS 


4.1  Results  of  Analysis  for  Cycles  with  an  Economizer  or  Intercooler 
The  transcritical  CO2  cycle  was  analyzed  with  the  effect  of  an  intercooler  or 
economizer  incorporated  into  the  base  cycle. 


4.1.1  Comparison  of  COPs  of  Base,  Economizer,  and  Intercooler  Cycles 


4. 1 . 1 . 1  Economizer  and  Base  Cycles 

Figure  4.1  shows  the  comparison  of  COPs  of  Economizer  cycle  (Figure  1.7)  and 
Base  cycle  (Figure  1.1).  For  this  analysis,  the  temperature  and  pressure  at  the  inlet  of  the 
lst-stage  compressor,  TE\  and  pE\,  and  the  air  temperature,  Tair,  are  assumed  to  be  25°C, 
3MPa  and  35°C,  respectively.  The  same  conditions  are  assumed  for  the  Base  cycle.  As 
shown  in  this  figure,  the  Economizer  cycle  has  larger  COPs  than  Base  cycle  for  the 
pressure  ratio  larger  than  1.6.  The  maximum  COPs  of  Economizer  and  Base  cycles  are 
3.564  and  3.043,  respectively.  The  maximum  COP  increases  by  17.1%  by  using  the 
Economizer  cycle. 

The  same  trend  has  been  observed  when  the  pressure  at  the  inlet  to  the  lst-stage 
compressor,  pE  1,  is  changed  but  the  other  input  variables  are  fixed.  It  can  be  seen  in 
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Figure  4.1  that  the  COP  of  the  Economizer  cycle  is  zero  from  values  of  pressure  ratio 
from  1.0  to  1.5.  This  means  that  the  COP  was  not  calculated  for  that  pressure  ratio  range. 
When  the  pressure  ratio  is  small  and  the  air  temperature  (heat  sink  in  the  gas  cooler)  is 
high,  the  enthalpy  at  the  outlet  of  the  gas  cooler  is  higher  than  that  of  the  vapor  phase  at 
the  intermediate  pressure.  Therefore,  the  CO2  becomes  the  superheated  vapor  and  there 
is  no  phase  separation  in  the  accumulator  after  the  1st  expansion.  From  this  figure  it  can 
be  seen  that  the  C02  becomes  two-phase  in  the  accumulator  for  the  pressure  ratios  larger 
than  1.5. 

It  can  also  be  seen  from  Figure  4.1  that  the  COP  of  the  Base  cycle  is  zero  for 
pressure  ratios  from  1.0  to  1.5.  This  means  that  the  COP  of  the  Base  cycle  was  not 
calculated  because  the  Economizer  cycle  is  not  applicable  at  this  range  of  the  pressure 
ratios. 


4, 1 . 1 .2  Intercooler  and  Base  Cycles 

Figure  4.2  shows  the  comparison  of  the  COP  of  the  Intercooler  cycle  (Figure 
1.10)  with  that  of  the  Base  cycle.  It  is  assumed  that  the  temperature  and  pressure  at  the 
inlet  of  the  lst-stage  compressor,  Tn,  pn,  and  the  air  temperature,  Tair,  are  constant  and 
those  values  are  25°C,  3.5MPa  and  35°C,  respectively.  Again,  the  same  conditions  are 
given  for  the  Base  cycle.  Identical  pressure  ratios  are  used  for  both  of  the  compressors 
and  varied  from  1.0  to  2.4  to  calculate  the  COP.  The  solid  line  in  the  legend  stands  for 
the  Base  cycle  and  the  dashed  line  the  Intercooler  cycle  in  the  figure.  It  is  observed  that 
the  COP  of  the  Intercooler  cycle  is  higher  than  that  of  the  Base  cycle  and  the  maximum 
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COP  of  4.018  occurs  at  the  pressure  ratio  of  1.58  in  the  case  of  the  Intercooler  cycle.  The 
maximum  COP  of  Base  cycle  is  3.618  at  a  pressure  ratio  of  1.56.  The  maximum  COP 
increases  by  1 1.1%  by  using  Intercooler  cycle. 

As  shown  in  the  figure,  the  COP  of  the  Intercooler  cycle  is  zero  from  values  of 
pressure  ratio  from  1.0  to  1.12.  This  is  due  to  the  fact  that  when  the  intermediate 
temperature,  7/2,  of  the  CO2  after  the  lst-stage  compressor  is  lower  than  Tair,  there  will  be 
no  heat  transfer  through  the  intercooler  and  no  calculation  of  the  COP  has  been 
performed.  From  this  figure  it  can  be  seen  that  7>  begins  to  exceed  Tair  for  pressure 
ratios  larger  than  1.14. 

In  addition,  the  COP  of  the  Base  cycle  is  zero  for  pressure  ratios  from  1.0  to  1.12. 
This  means  that  the  COP  of  the  Base  cycle  was  not  calculated  because  the  intercooler 
cycle  is  not  applicable  at  this  range  of  the  pressure  ratios. 

It  can  also  be  observed  in  Figure  4.2  that  the  COP  of  the  Intercooler  cycle 
increases  again  at  a  pressure  ratio  of  approximately  2.2,  then  decreases.  That  means  there 
is  the  2nd  peak  of  COP  in  the  Intercooler  cycle.  This  is  caused  by  the  decrease  of  the 
total  combined  work  of  the  1st-  and  2nd-stage  compressors.  Figure  4.3  shows  the 
comparison  of  the  COP  and  the  compression  work.  From  the  pressure  ratio  of  2.2  the 
total  combined  work  of  the  1st-  and  2nd-stage  compressors  becomes  less  and  this  results  in 
increase  of  COP  as  shown  in  Figure  4.3.  Starting  at  a  pressure  ratio  of  2.36,  the  COP 
decreases  again  as  the  compression  work  increases. 

As  discussed  above,  the  system  performance  can  be  improved  by  modifying  the 
Base  cycle  into  the  Economizer  cycle  or  the  Intercooler  cycle.  In  addition,  the  system 
efficiencies  of  the  Economizer  and  Intercooler  cycles  can  be  increased  much  further  if  the 
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optimum  operating  condition  is  selected.  In  the  next  sections,  the  analyses  are  discussed 
in  detail  about  the  effects  of  some  important  parameters  on  COPs  of  Economizer  and 
Intercooler  cycles  for  the  optimum  operating  condition. 

4.1.2  Effects  of  Pressure  Ratios  on  COP 

In  order  to  give  the  minimum  compression  work  in  the  Intercooler  cycle,  it  is 
common  to  give  the  same  pressure  ratio  across  the  two  compressors  [Wark,  1988].  In 
addition,  the  temperature  of  the  C02  at  the  outlet  of  the  lst-stage  compressor  (TIV)  is 
cooled  down  through  the  intercooler  to  the  temperature  of  C02  at  the  inlet  to  the  lsl-stage 
compressor  ( Tn ).  During  warm  weather  periods,  however,  the  intercooler  may  not  be 
able  to  reject  enough  heat  for  the  C02  temperature  to  be  reduced  to  Tn.  In  this  case, 
different  pressure  ratios  for  each  of  the  two  compressors  will  be  required  to  give  the 
minimum  compression  work.  Therefore,  it  is  necessary  to  evaluate  what  combination  of 
pressure  ratios  across  the  compressors  give  the  maximum  COP  for  a  given  operating 
condition. 

Figure  4.4  shows  the  behavior  of  the  COP  for  the  Intercooler  cycle.  The  C02 
pressure,  pn,  and  temperature,  Tn,  at  the  inlet  to  the  lst-stage  compressor  are  assumed  to 
be  3.5MPa  and  25°C,  respectively.  The  air  temperature,  Tair,  is  a  heat  sink  temperature 
and  assumed  to  be  35°C.  The  pressure,  pn,  at  the  outlet  of  the  2nd-stage  compressor,  is 
taken  to  be  lOMPa.  If  the  same  pressure  ratio  is  assumed  across  the  compressors,  the 
maximum  COP  of  3.842  occurs  at  pressure  ratio  of  1.69  and  intermediate  pressure  of 
5.92MPa.  The  COP  of  the  Base  cycle  for  this  condition  becomes  3.361.  As  shown  in 
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Figure  4.4,  however,  the  COP  continues  to  increase  if  the  intermediate  pressure  is 
increased  with  an  optimum  COP  of  3.956  occurring  at  intermediate  pressure  of  8.12MPa. 
In  this  case,  the  pressure  ratios  across  the  lst-stage  compressor  and  2nd-stage  compressor 
are  2.32  and  1.23,  respectively.  From  this  result  it  could  be  deduced  that  improvements 
of  COP  could  be  obtained  if  the  pressure  ratios  are  chosen  appropriately. 

As  shown  in  Figure  4.4,  there  is  sudden  increase  of  COP  at  around  the  pressure 
ratio  of  2.25  and  the  COP  reaches  its  maximum  point  at  the  pressure  ratio  of  2.32.  This 
behavior  of  COP  is  due  to  the  sudden  decrease  of  total  compression  work.  Figure  4.5 
shows  the  comparison  of  the  compression  works  for  this  case.  In  the  legend  of  the  figure, 
‘lst-stage  compression  work’,  ‘2nd-stage  compression  work’,  and  ‘Total  compression 
work’  stand  for  the  compression  work  of  the  lst-stage  compressor,  the  compression  work 
of  the  2nd-stage  compressor,  and  the  summation  of  the  1st-  and  2nd-stage  compression 
works,  respectively.  As  shown  in  this  figure,  the  compression  work  of  the  2nd-stage 
compressor  decreases  much  at  around  the  pressure  ratio  of  2.25  and  the  total  compression 
work  becomes  the  minimum  at  the  pressure  ratio  of  2.32  at  which  the  maximum  COP 
occurs.  At  this  point,  the  pressure  ratio  across  the  2nd-stage  compressor  is  1.23. 

The  behavior  of  the  compression  work  over  the  pressure  ratio  across  the 
compressors  is  due  to  the  thermodynamic  characteristics  of  the  carbon  dioxide.  Figure 
4.6  shows  the  processes  of  the  Intercooler  cycle  with  the  fixed  pressure  ratios  and  the 
different  pressure  ratios  on  the  pressure-enthalpy  diagram.  When  the  fixed  pressure 
ratios  of  1.69  are  applied,  the  intermediate  pressure  at  the  intercooler  becomes  5.92MPa. 
In  this  case,  the  paths  of  the  refrigeration  processes  are  11-11  ’j-I2  yi2 j-13 -14-11 .  The 
subscript  T  represents  the  case  of  the  fixed  pressure  ratio.  If  the  different  pressure 
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ratios,  2.32  for  the  lst-stage  compressor  and  1.23  for  the  2nd-stage  compressor,  are 
considered,  the  intermediate  pressure  becomes  8.12MPa.  Hence  the  paths  of  the 
processes  become  11-11  'd-12  ’ d-12d-13-14-U .  The  subscript  ‘d’  represents  the  case  of  the 
different  pressure  ratios.  The  total  compression  works  for  both  cases  are  as  follows. 

wf  ~  (^/r,  —  hn )+  {^i2,  ~  hir, ) 
and 

Wd  ~  ~  )+  l^/2rf  _  klTj  )’ 

Equations  (4-1)  and  (4-2)  represent  the  total  compression  work  of  the  cases  of  the 
fixed  pressure  ratio  and  different  pressure  ratio,  respectively.  The  first  and  second  terms 
in  the  right  hand  sides  of  Equations  (4-1)  and  (4-2)  stand  for  the  compression  works  of 
the  lst-stage  and  2nd-stage  compressors,  respectively.  As  discussed  above,  the  2  -stage 
compression  work  decreases  suddenly  at  the  pressure  ratio  of  2.32  across  the  1  -stage 
compressor.  As  shown  in  Figure  4.6,  the  enthalpy  increase  through  the  2  -stage 
compressor  becomes  very  small  at  that  point.  It  results  in  the  minimum  total 
compression  work.  Therefore,  the  COP  of  the  Intercooler  cycle  becomes  the  maximum 
with  the  constant  capacity  of  the  evaporator.  It  can  be  also  seen  in  Figure  4.6  that  there  is 
a  large  enthalpy  decrease  through  the  intercooler  at  the  intermediate  pressure  of  8.12MPa 
due  to  the  almost  zero  slope  of  the  temperature  curve  of  carbon  dioxide  with  temperature 
of  35°C  at  that  pressure,  even  though  the  figure  does  not  show  it  clearly  because  of  the 


(4-1) 


(4-2) 
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coarse  picture.  Therefore,  in  this  case,  most  of  the  heat  of  C02  in  the  system  is  removed 
through  the  intercooler  cycle  and  a  small  portion  of  the  heat  through  the  gas  cooler. 


4.1.3  Effect  of  Low-Side  Pressure  on  COP 
In  the  previous  section,  the  pressure  ratio  effect  on  the  COP  was  investigated.  In 
this  section,  identical  pressure  ratios  across  both  of  the  two  compressors  are  assumed  to 
illuminate  the  specific  thermodynamic  performance  difference  between  the  proposed 
transcritical  C02  cycles.  For  the  Base  cycle,  the  pressure  ratio  across  the  compressor  is 
square  of  the  pressure  ratio  of  the  modified  cycles  to  give  the  same  high-side  pressure 
( PE2  or  pn). 

4. 1.3.1  Economizer  Cycle 

Figure  4.7  shows  the  effects  of  the  C02  pressure  at  the  inlet  to  the  lst-stage 
compressor  (low-side  pressure),  pe\,  on  the  COP  of  the  Economizer  cycle.  It  can  be 
observed  from  this  figure  that  as  the  pe\  increases,  the  COP  increases  and  the  pressure 
ratio  across  the  lst-stage  compressor  at  which  the  maximum  COP  occurs  decreases. 
When  pe\  is  4.5MPa,  the  maximum  COP  of  5.835  occurs  at  the  pressure  ratio  of  1.36. 
When  pe\  is  2.5MPa,  the  maximum  COP  becomes  3.017  at  the  pressure  ratio  of  1.84. 
The  COP  can  increase  by  93.4%  by  increasing  the  C02  pressure  at  the  inlet  to  the  1st- 
stage  compressor  from  2.5MPa  to  4.5MPa. 

Actually,  the  control  of  the  low-side  pressure  as  a  function  of  operating 
conditions  is  a  difficult  task  to  achieve  because  the  low-side  pressure  strongly  depends  on 
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the  capacity,  heat  source  temperature,  and  geometry  of  the  evaporator.  However,  this  is  a 
parametric  study,  in  which  it  is  assumed  that  depending  on  the  evaporator  design,  a 
higher  the  low-side  pressure  may  be  achievable. 

Figure  4.8  shows  a  comparison  of  the  heat  removal  capacity  of  the  evaporator  in 
the  Economizer  cycle.  The  capacity  is  1 18.08kJ/kg  at  the  pressure  ratio  of  1.36  when  pE\ 
is  4.5MPa.  When  pE\  is  2.5MPa,  it  becomes  145.98kJ/kg  at  a  pressure  ratio  of  1.84.  The 
evaporator  heat  removal  capacity  is  decreased  by  19.1%  by  increasing  the  pressure  at  the 
inlet  to  the  lst-stage  compressor  from  2.5MPa  to  4.5MPa.  Table  4.1  shows  the 
compression  work  at  the  pressure  ratio  of  each  low-side  pressure  at  which  the  maximum 
COP  occurs.  As  shown  in  the  table,  the  compression  work  decreases  considerably,  by 
53.3%,  when  the  low-side  pressure  increases  from  2.5Mpa  to  4.5MPa.  Therefore,  the 
effect  of  increasing  the  pressure  at  the  inlet  to  the  evaporator  is  a  net  increase  in  the  COP, 
due  to  a  larger  decrease  of  the  compression  work  and  a  smaller  decrease  of  the  heat 
removal  capacity  of  the  evaporator. 

4. 1.3.2  Intercooler  Cycle 

Figure  4.9  shows  the  effect  of  the  pressure  at  the  inlet  to  the  lst-stage  compressor, 
pn,  on  the  COP  of  the  Intercooler  cycle.  The  effect  of  the  pressure  is  the  same  as  in  the 
case  of  Economizer  cycle.  As  shown  in  the  figure,  the  maximum  COP  increases  as  pn 
increases.  Also  the  pressure  ratio  at  which  this  maximum  COP  occurs  decreases  as  the 
pn  increases.  When  pn  is  4.5MPa,  the  maximum  COP  of  5.612  occurs  at  a  pressure  ratio 
of  1.38.  When  pn  is  2.5MPa,  the  maximum  COP  becomes  2.872  at  the  pressure  ratio  of 
1.88.  The  COP  is  thus  increased  by  100.7%  by  increasing  the  C02  pressure  at  the  inlet  to 
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the  lst-stage  compressor  from  2.5MPa  to  4.5MPa.  The  reason  for  the  increase  of  COP 
with  increase  of  Pn  is  the  same  as  in  the  case  of  the  Economizer  cycle.  It  is  due  to  a 
larger  decrease  of  the  compression  work  and  a  smaller  decrease  of  the  heat  removal 
capacity  of  the  evaporator. 


4.1.4  Effect  of  High-Side  Pressure  on  COP 

As  mentioned  above,  the  COP  of  the  transcritical  CO2  cycle  depends  strongly  on 
the  low-side  pressure  if  the  evaporator  is  selected  correctly.  In  addition,  the  COP  is 
dependent  on  pressure  ratio  split  between  the  two  compressors  in  case  of  the  two-stage 
compression  cycles.  In  this  section,  the  optimum  condition  for  pressure  at  the  gas  cooler 
is  investigated  for  the  Intercooler  cycle  at  a  given  air  temperature  and  CO2  pressure  at  the 
evaporator. 

The  COP  of  the  Intercooler  cycle  has  been  calculated  by  changing  the  high-side 
pressure  for  given  CO2  temperature,  Tn,  and  pressure,  pn,  at  the  inlet  to  the  lst-stage 
compressor  with  several  air  temperature,  ra,>.  The  values  of  Tn  and  pn  are  26.7°C  and 
3.5Mpa,  respectively.  Figures  4.10  through  4.14  show  the  change  of  the  maximum  COPs 
over  high-side  pressure  for  the  given  air  temperatures  of  27.8°C,  35°C,  40.6°C,  46.1°C 
and  51.7°C,  respectively.  In  the  legend  of  the  figures,  ‘Intercooler-l’  and  ‘Intercooler-2’ 
stand  for  the  Intercooler  cycle  with  identical  pressure  ratios  for  both  compressors  and 
different  pressure  ratios,  respectively.  The  calculation  has  been  performed  at  the  high- 
side  pressures  (Pi 2)  of  8MPa,  lOMPa  and  12MPa.  It  can  be  observed  from  these  figures 
that  the  maximum  COP  does  not  occur  at  the  same  high-side  pressure.  In  other  words, 
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there  is  a  specific  high-side  pressure  at  which  the  maximum  COP  occurs  for  a  given  air 
temperature.  When  the  air  temperature  is  low,  7’a,r=27.8°C,  the  maximum  COP  occurs  at 
the  high-side  pressure  of  8MPa.  At  an  intermediate  air  temperature,  35°C  <  Tair< 40.6°C, 
the  maximum  COP  occurs  at  lO.OMPa.  When  the  air  temperature  is  high,  it  occurs  at 
12MPa. 

This  phenomenon  happens  due  to  the  significant  change  of  the  enthalpy  of  the 
C02  over  temperature  drop  at  a  specific  pressure  after  the  critical  point  as  shown  in 
Figure  4.15.  Figure  4.15  shows  the  thermodynamic  characteristics  of  the  C02  on  the  p-h 
diagram.  One  of  the  assumptions  is  that  the  C02  temperature  drops  down  to  the  air 
temperature  due  to  the  heat  rejection  in  the  gas  cooler.  When  the  air  temperature  is 
27.8°C,  the  ratio  of  the  decrease  of  the  C02  enthalpy  through  the  gas  cooler  to  the 
compression  work  becomes  larger  at  the  high  pressure  of  8.0MPa  than  at  lO.OMPa  or  at 
12.0MPa.  Therefore,  when  the  air  temperature  is  27.8°C,  the  COP  becomes  the  largest 
at  the  high-side  pressure  of  8MPa  by  assuming  the  isenthalpic  expansion.  However, 
when  Tair  is  35°C  or  40.6°C,  the  largest  ratio  of  the  C02  enthalpy  change  to  the 
compression  work  occurs  at  lO.OMPa.  When  the  air  temperature  is  very  high,  46.1°C  or 

51.7°C,  it  occurs  at  12.0MPa  of  the  high-side  pressure. 

Some  other  cases  are  analyzed  under  similar  conditions  to  the  previous  cases, 
except  for  the  C02  temperature  (770  of  32.2°C  at  the  inlet  to  the  lst-stage  compressor. 
The  values  of  Tn  and pn  are  32.2°C  and  3.5Mpa,  respectively.  Figures  4.16  through  4.19 
show  the  change  of  the  maximum  COPs  over  high-side  pressure  for  the  given  air 
temperatures  of  35°C,  40.6°C,  46.1°C  and  51.7°C,  respectively.  As  shown  in  these 


128 


figures,  the  results  of  these  cases  demonstrate  similar  behavior  of  the  COP  to  the  previous 
cases. 

In  order  to  obtain  the  high-side  pressure  at  which  the  maximum  COP  occurs  at  a 
given  air  temperature,  the  calculation  of  the  COP  has  been  performed  again  by  increasing 
the  high-side  pressure  from  8MPa  to  12MPa  in  0.2MPa  increments.  The  different 
pressure  ratios  for  the  compressors  are  considered  by  variation  of  the  pressure  ratio 
across  the  lst-stage  compressor  from  1.0  to  3.5  by  0.02.  Figures  4.20  through  4.24  show 
the  behavior  of  the  COP  in  the  three-dimensional  axes.  The  x-axis,  y-axis  and  z-axis 
represent  the  high-side  pressure,  pressure  ratio,  and  COP,  respectively.  From  these 
figures,  it  can  be  seen  that  there  is  a  specific  high-side  pressure  at  which  the  maximum 
COP  of  the  cycle  occurs.  It  can  also  be  observed  that  there  is  a  pressure  ratio  across  the 
lst-stage  compressor  at  which  the  COP  becomes  the  maximum.  Table  4.2  shows  the 
maximum  COP  values  for  a  given  condition,  the  high-side  pressure  and  pressure  ratio 
across  the  lst-stage  compressor  at  which  the  maximum  COP  occurs. 

It  can  be  deduced  from  these  results  that  if  the  pressure  at  the  gas  cooler  is  chosen 
carefully  for  a  given  CO2  temperature  at  the  outlet  of  the  gas  cooler,  high  thermodynamic 
performance  of  the  cycle  can  be  obtained.  In  addition,  the  different  pressure  ratio  across 
the  lst-stage  compressor  should  be  also  decided  for  optimum  thermodynamic 
performance. 


4. 1 .5  Summary  for  Cycles  with  an  Economizer  or  Intercooler 
The  transcritical  CO2  cycles  with  an  economizer  or  an  intercooler  (Economizer 
cycle  and  Intercooler  cycle)  have  been  considered  for  analysis  of  their  thermodynamic 
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performances.  The  results  show  that  there  is  significant  improvement  of  the  COP  in  the 
Economizer  cycle  and  Intercooler  cycle.  It  is  also  observed  that  the  higher  the  low-side 
pressure  is,  the  higher  the  COP  of  the  C02  cycle.  The  effect  of  the  pressure  ratios  for  the 
compressors  has  also  been  investigated.  When  the  cycles  are  optimized,  different 
pressure  ratios  for  each  compressor  give  better  performance  than  using  identical  pressure 
ratios  across  the  both  compressors  during  warm  weather  period.  The  coefficient  of 
performance  of  the  transcritical  C02  cycle  with  two  compressors  including  intercooling 
(Intercooler  cycle)  depends  on  the  high-side  pressure  for  a  given  air  temperature  as  a  heat 
sink.  It  is  also  dependent  of  the  pressure  ratio  across  the  lsl-stage  compressor  for  a  given 
pressure  difference  between  the  low-side  pressure  and  the  high-side  pressure. 

4.2  Results  of  Analysis  for  Cycles  with  Work  Extraction 

4.2.1  Constant  Pressure  Ratios  across  Compressors 
Table  4.3  shows  the  results  of  the  analysis  for  various  C02  cycles  with  and 
without  the  ED-WOW.  In  the  analysis  shown  in  Table  4.3,  constant  pressure  ratios 
across  the  compressors  were  considered  for  the  Intercooler  cycle  (Case  5  and  6)  and  the 
Economizer  cycle  (Case  7).  The  percentage  of  the  pressure  increase  across  the  1st-  and 
2nd-stage  compressors  of  the  Intercooler  cycle  were  85.22%  and  14.78%,  respectively, 
which  equals  to  pressure  ratios  of  1.88  (lst-stage)  and  1.08  (2nd-stage)  at  an  evaporator 
pressure  of  5.087  MPa  and  a  C02  gas  cooler  pressure  of  10.34  MPa.  These  pressure 
ratios  were  chosen  because  they  gave  maximum  COP  at  upper  limit  design  point 
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conditions  of  an  evaporator  pressure  of  5.087  MPa  (equal  to  an  evaporation  temperature 
of  15°C),  a  CO2  gas  cooler  outlet  temperature  of  55°C,  and  a  CO2  gas  cooler  pressure  of 
10.34  MPa. 

For  the  Economizer  cycle,  equal  pressure  ratios  across  each  compressor  stage 
were  used.  This  is  due  to  the  fact  that  if  the  pressure  ratios  of  1.88  and  1.08  (as  used  in 
the  Intercooler  cycle)  are  used,  the  intermediate  pressure  becomes  higher  than  the  critical 
pressure  at  the  design  point  conditions.  Therefore,  a  single-phase  fluid  will  enter  the 
accumulator,  which  is  at  the  intermediate  pressure,  and  no  phase  separation  will  occur. 
Thus,  the  cycle  is  not  applicable  (N/A).  Even  by  assuming  equal  pressure  ratios  across 
the  first  and  second  stage  compressors,  for  some  operating  conditions  the  state  of  carbon 
dioxide  at  the  accumulator  inlet  is  compressed  liquid  (“Comp.  Liq.”)  or  superheated 
vapor  (“Sup.  Vapor”),  which  makes  the  cycle  not  applicable,  as  shown  in  Table  4.3.  In 
these  cases,  an  analysis  of  the  compression  work,  heat  removal  capacity  of  the 
evaporator,  and  COP  was  omitted  and  labeled  “N/A”. 

There  are  eight  kinds  of  operating  condition  labeled  “OC-1”  through  “OC-8”  for 
each  CO2  cycle.  “OC”  stands  for  the  operating  condition.  The  operating  condition  OC-8 
is  the  upper  limit  design  point  condition  (C02  evaporation  temperature  (TeVap)  =  15°C, 
C02  temperature  at  the  gas  cooler  outlet  (Tgas.ou,)  =  55°C,  and  high-side  pressure  (phigh)  = 
10.34  MPa).  For  the  first  four  operating  conditions,  OC-1  through  OC-4,  the  evaporator 
temperature  is  changed  while  the  other  parameters  are  fixed.  As  shown  in  the  table,  the 
COP  of  the  Base  cycle  (Case  1)  increases  with  an  increase  in  evaporation  temperature  (or 
pressure).  This  is  due  to  the  small  decrease  of  the  heat  removal  capacity  of  the 
evaporator  compared  to  the  large  decrease  of  the  compression  work  as  Tevap  increases. 
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For  the  next  five  conditions  (OC-4  through  OC-8),  Tgas.out  is  changed  while  the 
other  parameters  are  fixed.  As  shown  in  Table  4.3,  the  COP  of  the  Base  cycle  decreases 
with  the  increase  of  Tgas.out.  The  COP  becomes  0.49  at  the  upper  limit  design  point 
condition,  OC-8.  This  is  caused  by  the  significant  decrease  of  the  heat  removal  capacity 
of  the  evaporator. 

The  same  trend  is  observed  in  the  other  CO2  cycles  (Case  2  through  Case  7). 
When  the  COPs  of  the  CO2  cycles  are  compared  at  the  operating  condition  OC-8,  the 
IHX  cycle  with  the  ED-WOW  (Case  4)  gives  the  highest  COP,  1.798,  by  using  the  work 
output  of  the  ED-WOW  to  assist  the  compression  work.  This  is  due  to  the  relatively 
large  heat  removal  capacity  of  the  evaporator  caused  by  the  heat  transfer  from  hot  CO2  to 
cold  C02  through  the  internal  heat  exchanger.  As  shown  in  Table  4.3,  Case  4  gives  the 
highest  COP  when  Tgas.ou,  is  very  high  (50°C  and  55°C).  However,  when  Tgas.out  is  not  as 
high  (less  than  or  equal  to  45°C),  the  Intercooler  cycle  with  the  ED-WOW  (Case  6)  gives 
the  highest  COP.  This  is  due  to  relatively  small  compression  work  compared  to  that  of 
Case  4. 

Two  cases  are  considered  for  the  ED-WOW.  The  first  case  is  that  only  the  lower 
enthalpy  effect  of  the  ED-WOW  at  the  inlet  to  the  evaporator  is  considered.  Then,  the 
work  extraction  effect  as  well  as  the  lower  enthalpy  effect  has  been  considered  and  the 
extracted  work  has  been  assumed  to  be  used  to  run  the  compressor.  As  shown  in  Table 
4.3,  the  COPs  of  the  C02  cycle  are  increased  a  little  in  the  first  case.  However,  the 
system  performance  is  improved  significantly  when  both  effects  of  the  expansion  device 
are  considered  as  shown  in  the  table. 
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Based  on  this  analysis,  it  can  be  concluded  that  at  high  Tgas_0ut,  the  IHX  cycle  with 
the  ED-WOW  gives  the  highest  system  performance,  and  at  intermediate  Tgas-out,  the 
intercooler  cycle  with  the  ED-WOW  shows  the  best  system  performance.  It  can  also  be 
concluded  that  when  the  ED-WOW  is  used  instead  of  a  throttling  expansion  valve,  the 
system  performances  of  the  CO2  cycles  increase  significantly.  For  instance,  the  COP  of 
the  Base  cycle  is  increased  by  up  to  112%  with  the  ED-WOW. 


4.2.2  Different  Pressure  Ratios  across  Compressors 

Table  4.4  shows  the  results  of  the  analysis  for  the  Intercooler  cycle  (Case  8), 
Intercooler  cycle  with  the  ED-WOW  (Case  9),  and  Economizer  cycle  (Case  10)  with 
varying  pressure  ratios  across  the  first  and  second  stage  compressors.  For  each  operating 
condition,  the  pressure  ratios  across  the  first  and  second  stage  compressors  were  adjusted 
to  give  the  maximum  system  COP.  Table  4.4  indicates  the  optimum  pressure  ratios  and 
the  maximum  COPs  for  the  cycles. 

By  comparing  the  COPs  in  Table  4.3  and  Table  4.4,  it  can  be  observed  that  the 
COPs  of  the  Intercooler  and  Economizer  cycles  can  be  significantly  increased  by 
selecting  the  optimum  pressure  ratios  for  the  given  operating  conditions.  However,  to 
independently  control  the  intermediate  pressure  in  two-stage  compression  systems  as  a 
function  of  operating  conditions  may  be  a  difficult  task  to  achieve  in  actual  equipment. 

As  discussed  in  the  previous  section,  the  system  performances  of  the  transcritical 
CO2  cycles  are  increased  considerably  with  the  ED-WOW  as  shown  in  Table  4.4. 
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4.3  Results  of  Theoretical  Model  for  ED-WOW  Device 

Analysis  of  the  ED-WOW  has  been  performed  using  the  developed  simulation 
model.  As  discussed  before,  the  ED-WOW  consists  of  flat  twin-cylinder  in  which  the 
pistons  are  on  out-of-phase  firing  order.  In  the  analysis  it  is  assumed  that  there  is  no 
time-delay  of  the  valves  to  respond  to  input  signals.  This  assumption  gives  conservatism 
in  terms  of  the  CO2  mass  flow  into  the  cylinder  of  the  ED-WOW. 

Four  revolution  speeds  of  crankshaft  are  considered.  Those  are  120  rpm,  240 
rpm,  360  rpm,  and  600  rpm.  The  one  revolution  of  the  crankshaft  takes  0.5  s  for  120 
rpm,  0.25  s  for  240  rpm,  0.167  s  for  360  rpm,  and  0.1  s  for  600  rpm,  respectively.  The 
input  values  for  the  parameters  are  shown  in  Table  4.5.  Among  the  input  values,  the 
crown  height  of  the  piston  ring  and  friction  coefficient  of  the  cylinder  are  assumed  to  be 
14.9  pm  and  0.08.  As  discussed  before,  it  is  considered  that  the  valves  open  and  close 
instantaneously  and  the  LOE  is  fixed  in  this  analysis.  Figure  4.25  shows  the  piston 
speeds  of  the  ED-WOW  over  its  crank  angle.  As  expected,  the  piston  speed  becomes 
faster  as  the  revolution  speed  increases. 

The  analysis  has  been  done  for  several  revolutions  of  the  crankshaft  to  remove 
any  uncertainties  which  may  exist  in  the  input  values.  As  shown  in  figure  4.26,  the  same 
trends  of  the  pressure  behavior  are  repeated  for  120  RPM  and  240  RPM  during  the  2  , 
3rd,  4th,  and  5th  revolutions.  Therefore,  the  second  revolution  (360°  -720°)  of  the 
crankshaft  is  the  focused  of  the  current  analysis. 


134 


4.3.1  Chamber  Pressure 

Figure  4.27  shows  the  behavior  of  the  CO2  pressure  in  the  chamber  over  the  crank 
angle.  It  can  be  seen  in  this  figure  that  the  low  revolution  speeds  of  this  expansion 
device,  120  rpm  and  240  rpm,  are  better  than  other  faster  speeds  to  get  the  inlet  high 
pressure  inside  the  cylinder.  When  the  speed  becomes  equal  to  or  faster  than  360  rpm, 
the  valve  opening  time  is  too  short  for  the  carbon  dioxide  to  fill  the  chamber  and  to  reach 
the  pressure  upstream  of  the  inlet  valve  because  of  the  characteristics  of  the  valves  as 
described  in  Chapter  3.  In  addition,  the  power  loss  becomes  larger  due  to  the  increase  of 
friction  force  with  the  faster  speed  and  it  reduces  the  pressure  and  temperature  inside  the 
chamber  of  the  ED-WOW. 

When  the  speed  of  device  is  240  rpm,  there  is  an  increase  of  pressure  in  the 
cylinder  during  the  exhaust  process  as  shown  in  Figure  4.27.  This  is  because  of  small 
orifice  diameter  of  the  exhaust  valve  and  relatively  short  time  period  to  sweep  out  the 
expanded  CO2  from  the  cylinder  compared  to  the  case  of  120  rpm.  Therefore,  the 
revolution  speed  of  the  device  of  120  rpm  is  considered  as  the  optimum  operating  speed 
for  our  application.  As  discussed  in  Chapter  3,  the  speed  of  120  rpm  was  also  considered 
as  the  optimum  operating  speed  based  on  the  experiments.  Only  the  case  of  120  rpm  is 
analyzed  and  discussed  in  detail  in  the  rest  of  this  section. 

In  the  case  of  120  rpm  the  pressure  in  the  cylinder  reaches  the  intake  pressure 
very  quickly  and  remains  constant  during  the  rest  of  the  intake  period  as  shown  in  Figure 
4.27.  Then,  it  decreases  suddenly  as  C02  expands  in  the  expansion  process.  At  the  end 
of  the  expansion,  the  CO2  pressure  becomes  about  3.4MPa  that  is  the  desired  operating 
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low-side  pressure.  This  means  that  the  location  of  the  expansion  initiation  of  8.01  mm 
from  the  top  of  the  chamber  (5.21  mm  from  TDC)  was  determined  reasonably  well. 

The  pressure  oscillates  at  the  end  of  the  exhaust  period.  This  is  due  to  the 
dynamic  behavior  of  the  C02  flow.  The  pressure  in  the  chamber  reaches  the  pressure  at 
the  inlet  of  the  evaporator  (or  in  the  crank  housing)  at  the  angle  of  the  crankshaft  of 
around  680  degree.  At  the  next  time  step,  the  chamber  pressure  becomes  less  than  that  of 
crank  housing  due  to  the  C02  discharge  so  that  there  is  no  mass  discharge  because  there 
is  no  driving  force  to  sweep  out  the  expanded  C02  from  the  cylinder.  As  the  piston 
moves  up  further,  however,  the  chamber  volume  decreases  and  it  leads  to  increase  of  the 
chamber  pressure  and  C02  mass  is  discharged.  This  behavior  continues  until  the  end  of 
the  exhaust  process. 


4.3.2  Work  Done  by  C02 

This  simulation  model  predicts  the  work  done  by  C02  as  sown  in  Figure  4.28. 
Positive  work  stands  for  the  work  done  into  the  C02  and  negative  work  means  the  work 
done  by  C02.  It  can  be  seen  in  this  figure  that  the  work  is  produced  during  the  intake  and 
expansion  processes  but  work  is  required  to  push  the  piston  up  and  to  sweep  out  C02 
from  the  device  to  the  inlet  of  the  evaporator  through  the  exhaust  valve  during  the 
exhaust  period,  as  expected.  However,  the  work  output  is  much  larger  than  the  work 
required  so  that  work  is  generated  through  the  device.  The  work  generated  is  40.6 
J/piston-revolution  when  the  device  is  being  operated  at  120  rpm  with  the  high-side 
pressure  of  10.3  MPa  and  the  low-side  pressure  of  3.4  MPa. 
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4.3.3  Minimum  Film  Thickness  Under  the  Piston  Ring 
Figure  4.29  shows  the  minimum  film  thickness  under  the  piston  ring.  As  shown 
in  the  figure,  the  film  thickness  is  constant  over  the  whole  processes.  Its  value  is  too 
small  and  the  same  as  the  composite  roughness  of  the  cylinder.  This  means  that  there  is 
always  the  boundary  lubrication  but  no  hydrodynamic  lubrication.  There  is  always  the 
piston  scuffing.  This  boundary  lubrication  occurs  due  to  the  very  low  piston  speed  in  this 
case,  as  shown  in  Figure  4.25.  The  low  piston  speed  leads  to  low  pressure  under  the 
piston  ring  so  that  the  internal  gas  pressure  is  always  high  enough  to  push  the  piston  ring 
radially  and  outwardly  and  the  ring  touches  the  surface  of  the  cylinder  wall  always. 


4.3.4  Power  Loss 

As  discussed  in  the  previous  section,  the  piston  ring  scuffing  occurs  during  the 
whole  processes.  Therefore,  the  power  loss  is  due  to  only  the  friction  of  the  piston  ring 
as  shown  in  Equation  (2-63).  There  is  no  power  loss  due  to  the  squeeze  motion  of  the 
piston  ring  because  the  ring  always  touches  the  surface  of  the  cylinder  wall.  The  power 
loss  becomes  small  as  shown  in  Figure  4.30.  The  total  power  loss  is  about  2.1  J/piston- 
revolution  in  this  case. 

There  is  a  sharp  decrease  of  the  power  loss  as  the  expansion  process  begins.  This 
is  due  to  the  sudden  decrease  of  the  pressure  inside  the  device  at  that  point  as  shown  in 
Figure  4.27.  The  power  loss  depends  on  the  pressure  inside  the  expansion  device  as 
described  in  Equations  (2-61)  through  (2-63). 
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4.4  Experimental  Results  for  Cycle  with  ED-WOW  Device 
The  ED-WOW  has  been  implemented  into  a  prototype  transcritical  C02  cycle  as 
shown  in  Figure  2.35.  As  shown  in  the  figure,  the  system  consists  of  an  evaporator,  a 
compressor,  a  gas  cooler,  an  internal  heat  exchanger,  and  an  ED-WOW.  In  addition, 
there  is  an  ordinary  expansion  valve.  Experiments  have  been  performed  to  investigate 
the  system  performance  and  the  concept  of  the  ED-WOW  has  been  proved  by  the 
improvement  of  the  system  efficiency. 

4.4.1  Experimental  Test  Conditions 

Experiments  have  been  performed  for  three  different  operating  conditions.  Table 
4.6  shows  the  experimental  conditions  for  the  cases.  In  the  table,  Tindoor  and  Toutdoor  stand 
for  the  indoor  temperature  and  outdoor  temperature,  respectively.  The  same  indoor  and 
outdoor  temperatures  have  been  considered  for  all  cases  and  they  are  20  C  and  35  C, 
respectively.  Pc>0ut,  Pcjn,  Tc,ou„  and  Tc,in  indicate  the  compressor  discharge  pressure, 
compressor  suction  pressure,  compressor  discharge  temperature,  and  compressor  suction 
temperature,  respectively.  The  compressor  discharge  pressure  considered  is  6620  kPa  to 
7971  kPa  and  the  compressor  suction  pressure  changes  from  2528  kPa  to  2771  kPa.  The 
compressor  discharge  temperatures  considered  are  from  88°C  to  118°C  and  the 
compressor  suction  temperature  is  approximately  20°C  for  all  cases.  tnCo2  represents  the 

C02  mass  flow  rate  in  the  system  and  the  range  of  9.7  kg/s  to  13.5  g/s  has  been 
determined  based  on  the  compressor  capacity.  g)ed  2nd  Atjntake  are  the  rotational  speed 
and  the  intake  valve  opening  time  of  the  ED-WOW. 
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Pgc.ini  Pgc.outj  PED,in>  PED,ouh  2nd  Pevap.out  are  the  pressures  at  the  inlet  and  outlet  of 
the  gas  cooler,  at  the  inlet  and  outlet  of  the  ED-WOW,  and  at  the  outlet  of  the  evaporator, 
respectively.  Tgc,in,  Tgc>out,  TED,in,  and  TeVap,out  are  the  temperature  at  the  inlet  and  outlet  of 
the  gas  cooler,  at  the  inlet  of  the  ED-WOW,  and  at  the  outlet  of  the  evaporator, 
respectively. 


4.4.2  Cycle  Performance 

Figure  4.31  shows  the  p-h  diagram  of  the  processes  of  the  system  for  all  cases. 
The  numbers  in  the  figure  show  the  state  at  each  point.  Those  states  are  also  shown  in 
Figure  2.35.  As  shown  in  Table  4.6  and  Figure  4.31,  the  pressure  at  the  outlet  of  the 
compressor,  7971  kPa,  is  higher  than  the  critical  pressure  of  C02  (7377  kPa)  in  Case  3. 

There  is  a  relatively  large  pressure  drop  from  state  4,  the  outlet  of  the  gas  cooler, 
to  state  5,  the  inlet  of  the  ED-WOW.  This  is  due  to  the  many  valves  and  fittings  and  the 
internal  heat  exchanger  (IHX)  between  these  stations  as  shown  in  Figures  2.35  and  2.36. 
These  valves  have  been  arranged  to  measure  the  mass  flow  rate  of  C02  during  the 
experiment.  It  can  also  be  seen  in  the  figure  that  pressure  decreases  relatively  much  from 
the  outlet  of  the  evaporator  to  the  inlet  of  the  compressor.  This  is  caused  by  the  internal 
heat  exchanger  and  the  accumulator  between  the  evaporator  and  compressor. 

The  process  from  state  5  to  state  7  stands  for  the  case  of  the  isenthalpic  expansion 
process  (assuming  the  expansion  through  the  ordinary  expansion  valve)  in  Figure  4.31. 
As  expected,  the  enthalpy  decreases  at  the  inlet  of  the  evaporator  when  C02  expands 
through  the  ED-WOW  even  though  the  amount  of  the  enthalpy  decrease  does  not  seem 


great 
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As  previously  mentioned,  the  work  generated  through  the  ED-WOW  has  been 
measured  indirectly  by  measuring  the  power  of  the  hydraulic  pump.  The  efficiency  of 
the  pump  was  measured  at  approximately  50%  at  an  ED-WOW  speed  of  300  RPM.  As 
the  design  speed  of  the  ED-WOW  in  the  current  tests  is  120  RPM,  the  pump  efficiency 
can  be  expected  to  be  less  than  this  value.  As  the  actual  pump  efficiency  could  not  be 
measured  at  this  low  speed,  the  current  assumption  of  the  pump  efficiency  will  be  taken 
at  50%  which  will  thus  provide  a  conservative  estimate  of  the  work  extracted  from  the 
ED-WOW. 

The  measured  works  produced  by  the  ED-WOW  are  23.2  W,  28.7  W,  and  34.8  W 
for  Cases  1,  2,  and  3,  respectively.  As  pressure  difference  increases  across  the  ED- 
WOW,  the  work  produced  increases.  The  measured  isentropic  efficiency  of  the  ED- 
WOW  (8S, ed)  is  about  10%  as  shown  in  Table  4.7. 

Table  4.7  and  Table  4.8  show  the  results  of  the  experiments.  PED,in,  Ped.oui,  and 
Ted  in  stand  for  the  pressures  at  the  inlet  and  outlet  of  the  ED-WOW  and  the  temperature 
at  the  inlet  of  the  ED-WOW,  respectively.  Qisemh  and  QED  represent  the  evaporator 

capacity  when  C02  expands  through  the  ordinary  expansion  valve  and  the  ED-WOW, 
respectively.  As  discussed  above,  the  evaporator  capacity  increases  with  the  C02 
expansion  through  the  ED-WOW.  Wcomp  and  WED  illustrate  the  compressor  work  and  the 

work  produced  through  the  ED-WOW,  respectively.  COPexv  is  the  COP  based  on  the 
isenthalpic  expansion.  COPed-i  and  COPed-2  are  the  COPs  based  on  the  expansion 
through  the  ED-WOW.  COPEd-i  is  the  COP  of  the  system  when  only  the  lower  enthalpy 
effect  is  considered.  When  not  only  the  lower  enthalpy  effect  but  also  work  output  effect 
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is  considered,  COPed-2  is  calculated  by  reducing  the  compressor  work  input.  This  means 
that  the  work  produced  through  the  ED-WOW  is  considered  to  use  to  run  the  compressor, 
some  part  of  the  compression  work.  So,  it  saves  some  compression  work.  Those  COPs 
are  calculated  as  follows: 
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The  concept  of  the  ED-WOW  to  improve  system  performance  has  been  proven 
experimentally.  As  shown  in  Tables  4.7  and  4.8,  the  COP  increases  with  the  ED-WOW. 
When  only  the  effect  of  lower  enthalpy  is  considered,  the  COP  increases  by  3.3%  to  5.4% 
compared  to  the  COP  based  on  the  isenthalpic  expansion.  If  both  of  the  lower  enthalpy 
and  work  producing  effects  are  considered,  COP  are  increased  much  by  6.7%  to  9.9%. 
Again,  as  mentioned  previously,  the  hydraulic  pump  efficiency  that  is  used  to  calculate 
the  work  output  of  the  ED-WOW  is  a  conservative  estimate.  If  a  more  accurate 
efficiency  at  the  given  condition  is  used,  it  is  likely  the  work  output  through  the  ED- 
WOW  will  be  larger  than  the  value  resulting  in  more  improvement  in  the  system 
performance. 

As  shown  in  Table  4.7,  the  capacity  of  the  evaporator  increases  with  increase  of 
the  high-side  pressure  (compressor  discharge  pressure).  In  addition,  the  work  produced 
through  the  ED-WOW  is  increased  when  Pc>out  increases.  However,  the  percentage  of  the 
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increase  of  the  evaporator  and  COPed-i  decreases  with  the  increase  of  the  Pc,out-  The 
largest  percentage  for  the  increase  of  Qed  occurs  in  Case  1  and  for  the  increase  of 
COPed-i  also  in  Case  1.  When  both  of  the  lower  enthalpy  and  work  output  effects  are 
considered,  there  is  the  same  tendency  that  the  system  performance  improvement 
decreases  with  the  increase  of  Pc,out-  The  largest  improvement  of  the  system  performance 
occurs  in  Case  1  (lowest  compressor  discharge  pressure,  6621  kPa)  and  it  is  9.9%.  When 
the  compressor  discharge  pressure  is  highest  (7971  kPa,  Case  3),  the  system  efficiency 
improvement  becomes  the  smallest  and  it  is  6.7%  even  though  the  absolute  amounts  of 
enthalpy  decrease  at  the  inlet  of  the  evaporator  and  work  output  through  the  ED-WOW 
are  the  largest. 


4.4.3  Uncertainty  Analysis 

An  uncertainty  analysis  was  conducted  for  the  results  of  the  experiments.  The 
Kline  and  McClintock  method  [Kline  and  McClintock,  1953]  was  used  which  sums  the 
square  of  errors: 

(4-6) 
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where  n^z^and  wz  are  the  total  uncertainty  associated  with  the  dependent  variable  A, 
the  independent  variable  which  affects  the  dependent  variable  A,  and  the  uncertainty  of 
the  variable  z{ ,  respectively.  Equation  (4-6)  is  referred  to  as  the  second-power  equation. 

The  uncertainties  of  the  measurements  are  given  in  Table  4.9. 

The  work  generated  through  the  ED-WOW  is  calculated  as: 
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W  = 

VY  ED 


V&P„  V(p 

pump,o  P pump,i ) 


where  V ,  APpump ,  r)pump ,  Ppumpn ,  and  Ppumpi  are  volumetric  flow  rate  of  water,  pressure 

difference  across  the  hydraulic  pump,  the  efficiency  of  the  pump,  the  static  pressure  at  the 
outlet  of  the  pump,  and  the  static  pressure  at  the  inlet  of  the  pump,  respectively. 
Actually,  there  is  a  large  uncertainty  in  the  efficiency  of  the  pump.  However,  the 
uncertainty  analysis  has  been  conducted  for  the  given  pump  efficiency  (fixed  pump 
efficiency). 

The  total  uncertainty  for  the  generated  work  becomes: 
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where  ww  ,wv,wp  ,  and  wp  represent  the  uncertainties  for  the  work  produced 


through  the  ED-WOW,  volumetric  water  flow  rate,  pressure  at  the  outlet  of  the  pump, 
and  pressure  at  the  inlet  of  the  pump,  respectively. 

The  ideal  work  output  in  an  isentropic  expansion  process  is: 

Ws,ed  =  m{hED  i  -  hs  ED  o )  (4-9) 

where  m,hEDi,  and  hsEDo  stand  for  CO2  mass  flow  rate,  enthalpy  at  the  inlet  of  the  ED- 

WOW,  and  enthalpy  at  the  outlet  of  the  ED-WOW  for  the  isentropic  expansion  process, 
respectively.  The  enthalpy  at  the  outlet  of  the  ED-WOW  is  a  function  of  the  pressure  at 
the  outlet  of  the  ED-WOW  and  the  entropy  at  the  inlet  of  the  ED-WOW: 


hs,ED,o  ~  f  ( PeD.o  >  SED,i  ) 


(4-10) 
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where  pEDo  and  sED ,  are  the  pressure  at  the  outlet  of  the  ED- WOW  and  the  entropy  at 
the  inlet  of  the  ED-WOW,  respectively.  The  entropy  at  the  inlet  of  the  ED-WOW  is 


given  by: 


SED,i  ~  f  {PED,i  ’  ^ED,i  ) 


(4-11) 


where  pFD , ,  and  TED  .  are  the  pressure  and  temperature  at  the  inlet  of  the  ED-WOW, 
respectively.  Also,  the  enthalpy  at  the  inlet  of  the  ED-WOW  is  a  function  of  pressure  and 


temperature  is  given  by: 


h-EDJ  ~  /{Ped.i^edj)- 
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The  uncertainty  of  the  ideal  work  output  through  the  ED-WOW  becomes: 
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where  w-  ,  w . ,  w.  ,  and  w  .  are  the  uncertainties  of  the  ideal  work  produced 

through  the  ED-WOW,  C02  mass  flow  rate,  enthalpy  at  the  inlet  of  the  ED-WOW,  and 
enthalpy  at  the  outlet  of  the  ED-WOW  for  the  isentropic  expansion  process,  respectively. 
The  uncertainty  for  the  enthalpy  at  the  inlet  of  the  ED-WOW  is  calculated  as  follows: 
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where  w  ,  and  wT  are  the  uncertainties  for  the  pressure  and  temperature  at  the  inlet 

PED.i  1  ED  j 


of  the  ED-WOW,  respectively.  The  enthalpy,  hEDi ,  is  a  relative  value  to  a  reference 
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point.  For  conservatism  and  convenience,  the  enthalpy  difference  of  hED  i  -  hsEDo  is  used 
instead  of  hED  i .  Hence,  the  uncertainty  becomes: 


The  uncertainty  for  the  enthalpy  at  the  outlet  of  the  ED- WOW  is: 


where  w  ,  and  w  are  the  uncertainties  for  the  pressure  at  the  outlet  of  the  ED-WOW 

PED.o  *ED.i  r 

and  entropy  at  the  inlet  of  the  ED-WOW,  respectively.  Again,  the  enthalpy  difference  of 
hEDi  -  hs  ED  o  is  used  instead  of  hs  ED  o .  The  uncertainty  is  written  as: 


(4-17) 


The  uncertainty  at  the  inlet  of  the  ED-WOW  becomes: 


i 


(4-18) 


The  total  uncertainty  for  the  ideal  work  produced  through  the  ED-OW  can  be  obtained  by 
substituting  the  uncertainties  of  Equations  (4-15),  (4-17)  and  (4-18)  in  Equation  (4-13). 
The  isentropic  efficiency  of  the  ED-WOW  is  calculated  as  follows: 

(4-19) 

Ws.ED 

The  uncertainty  of  the  isentropic  efficiency  becomes: 


145 


(4-20) 


The  heat  removal  capacity  of  the  evaporator  is  expressed  as: 

QeD  =  ™{hevap,o  -  Kvap.i  )  (4'2 1 } 

where  QED ,  hemp  0 ,  and  hevap  i  are  the  heat  capacity  of  the  evaporator,  enthalpy  at  the  outlet 
of  the  evaporator,  and  enthalpy  at  the  inlet  of  the  evaporator,  respectively.  The  enthalpy 
at  the  inlet  of  the  evaporator  is  given  by: 

™Kap.i  =  ™hED,<  ~  Wed  (4'22) 


Hence,  Equation  (4-22)  can  be  written  as: 

Qed  =  ™(Kap,o  -  hED.i  )+WED.  (4-23) 

The  uncertainty  of  the  heat  capacity  of  the  evaporator  becomes: 

wQa)  =  \(hevap<0  ~  hEDJ  K  )2  +  )  +  (wwED  }  f2  (4"24) 

where  w  and  w,  are  the  uncertainties  of  the  heat  capacity  of  the  evaporator  and  the 

Qed  nevap,o 

enthalpy  at  the  outlet  of  the  evaporator,  respectively.  The  enthalpy  at  the  outlet  of  the 
evaporator  is  a  function  of  pressure  and  temperature  at  the  outlet  of  the  evaporator  as 


follows: 


^ evap,o  =  f(p  evap,o 


(4-25) 


where  pevap0  and  Tevap0  are  the  pressure  and  temperature  at  the  outlet  of  the  evaporator, 
respectively.  Therefore,  the  uncertainty  of  the  enthalpy  at  the  outlet  of  the  evaporator  can 


be  obtained  as  follows: 
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(4-26) 


where  wh  ,  wn  ,  and  wT  are  the  uncertainty  of  the  enthalpy,  pressure,  and 

nevap,o  Pevap.o  ievap,n 

temperature  at  the  outlet  of  the  evaporator,  respectively.  Again,  for  conservatism  and 
convenience,  the  enthalpy  difference  of  hevapo  -  hED  i  is  used  instead  of  hevap  o .  Hence, 


the  uncertainty  becomes: 


(4-27) 


The  uncertainty  of  the  heat  capacity  of  the  evaporator  can  be  obtained  by  substituting  the 
correct  uncertainty  values  into  Equation  (4-24). 

The  compression  work  is  expressed  as  follows: 

™  (4-28) 

where  Wcomp ,  t,  and  xn  are  the  compression  work  rate,  torque  of  the  compressor,  and 

rotational  speed  of  the  compressor,  respectively.  The  uncertainty  of  the  compression 
work  becomes: 


where  w ^  ,  wr,  and  ww  stand  for  the  uncertainties  of  the  torque  and  rotational  speed  of 

”comp 


the  compressor,  respectively. 
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The  COP  is  calculated  by  dividing  the  heat  capacity  of  the  evaporator  by  the 
compression  work  as  follows: 

COP  =  .  ®ED  - — .  (4-30) 

Wcomp  ~WED 

The  uncertainty  of  the  COP  is  obtained  as: 

1 

2 

(4-31) 


w, 


COP 


COP 


'w.  ' 2 

r  o 

icevaj 

^  Qevap  J 


W, 


W„ 


W  -W 

^  ''  comp  '  ED  J 


W  ■ 
WtD 


\Wcomp  WED  J 


where  wCOP  is  the  uncertainty  of  the  COP. 

Table  4.10  shows  the  calculated  uncertainties  for  Case  3.  As  shown  in  the  table, 
the  uncertainty  of  the  work  generated  through  the  ED-WOW  is  relatively  large.  This  is 
due  to  the  fact  that  the  instrumentations  that  were  used  to  measure  the  pressures  and 
water  flow  rate  were  oversized  for  the  recorded  measurements.  In  order  to  reduce  the 
uncertainties,  instrumentations  with  appropriate  measurements  ranges  need  to  be  used. 
The  large  uncertainty  value  of  the  isentropic  efficiency  of  the  ED-WOW  is  caused  by  the 
large  uncertainty  of  the  work  produced  through  the  ED-WOW. 


4.5  Analysis  of  the  ED-WOW  as  Installed  in  Cycle 
The  ED-WOW  has  been  analyzed  with  the  theoretical  model  for  the  experiments 
(Case  1  through  Case  3).  For  more  realistic  simulation  the  time-delay  of  0.053  s  is 
assumed  for  the  intake  valve  to  respond  to  the  input  signal  for  closing  in  this  analysis.  As 
discussed  in  Chapter  3,  the  time-delays  of  the  valves  were  measured  experimentally. 
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Figure  4.32  shows  the  results  of  the  simulation  for  the  pressure  behavior  inside 
the  cylinder  of  the  ED-WOW.  Table  4.11  also  indicates  some  results  of  the  simulation. 
In  Table  4.11  the  first  row  represents  the  case  number  and  the  next  four  rows  stand  for 
the  ED-WOW  inlet  pressure,  ED-WOW  outlet  pressure,  CO2  mass  flow  rate,  work 
produced  through  the  ED-WOW  in  the  experiments.  The  rest  of  the  rows  in  Table  4.11 
indicate  the  predicted  values. 

As  shown  in  Figure  4.32,  the  predicted  pressures  at  the  end  of  the  expansion  are 
higher  than  the  experiments  for  all  cases.  This  means  that  the  predicted  density  inside  the 
cylinder  at  that  point  is  higher  than  the  experiment.  The  predicted  density  should  be 
reduced  to  decrease  the  pressure  at  the  end  of  the  expansion.  This  means  that  there  is 
some  leakage  of  C02  from  the  cylinder  to  the  crank  housing  through  the  piston  ring 
during  the  experiments.  This  fact  can  be  ascertained  by  comparing  the  C02  mass  flow 
rates  of  the  experiments  and  simulations  shown  in  the  4th  and  7th  rows  in  Table  4.1 1.  The 
value  of  mass  flow  rate  in  the  7th  row  has  been  predicted  by  considering  zero  leakage  of 
C02  through  the  piston  ring.  As  shown  in  the  table,  the  measured  mass  flow  rates  are 
larger  than  the  predicted  ones. 

The  simulations  have  been  repeated  by  considering  some  leakage  of  C02  through 
the  piston  ring  during  the  intake  and  expansion  processes  as  shown  in  Table  4.11.  The 
values  of  the  10th  row  through  15th  row  are  predicted  with  the  assumption  of  the  leakage 
of  C02.  The  rates  of  the  leakage  have  been  determined  by  the  trial  and  error  method  in 
order  to  make  the  similar  mass  flow  rates  with  the  experiments  as  shown  in  the  4th  and 
1 1th  rows.  With  some  leakage  of  C02,  the  predicted  pressures  at  the  end  of  expansion  are 
decreased  much  compared  to  the  values  of  the  zero  leakage  assumption  and  those  values 
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are  very  similar  to  the  measured  ones  especially  for  Case  1  and  Case  2  as  shown  in  the  3 
and  10th  rows  in  Table  4.11.  The  predicted  pressure  is  still  a  bit  higher  than  the  measured 
one  for  Case  3.  The  general  results  and  trends  of  the  simulation  with  the  assumption  of 
the  leakage  are  in  a  good  agreement  with  the  experiments.  Therefore,  it  can  be  deduced 
that  there  is  a  leakage  of  C02  from  the  cylinder  of  the  ED-WOW  to  the  crank  housing 
through  the  piston  ring  during  the  intake  and  expansion  processes.  The  predicted  C02 
leakage  rate  and  its  percentage  to  the  total  C02  flow  rate  are  shown  in  the  13  and  14 
rows  in  Table  4.11.  The  range  of  the  leakage  rate  is  about  31.3%  to  39.6%  of  the  C02 
mass  flow  rate.  Figure  4.33  shows  the  pressure  behavior  inside  the  cylinder  with  the 
assumption  of  some  leakage  of  C02. 

The  leakage  rate  is  4.21  g/s  and  its  percentage  is  31.3%  of  the  C02  flow  rate  for 
Case  3.  The  percentage  seems  to  be  high  and  thus,  it  was  checked  if  the  predicted 
leakage  rate  is  reasonable  by  examining  the  leakage  area.  The  leakage  area  is  calculated 
based  on  the  leakage  flow  rate. 

A  two-phase  state  is  assumed  in  the  cylinder  of  the  ED-WOW  and  two-phase 
flow  is  considered  for  the  leakage  through  the  piston  ring  to  the  crank-housing.  Equation 
(2-40)  is  used  for  the  leakage  calculation.  The  pressure  of  the  crank-housing  is  assumed 
to  be  3.4  MPa.  Two  cylinder  pressures  are  considered.  They  are  7  MPa  and  4  MPa.  In 
addition,  two  vapor  qualities,  x  =  0.9  and  x  =  0.5 ,  are  examined.  These  conditions  will 
confine  the  possible  operating  condition  of  the  ED-WOW  during  the  experiments.  Table 
4.12  shows  the  calculation  conditions  and  the  results  for  Case  3. 

The  mass  velocities  are  calculated  by  using  Equation  (2-40)  as  mentioned.  Then, 
the  flow  area  for  the  mass  flow  rate  of  4.21  g/s  is  calculated  and  is  given  in  the  3rd  row  of 
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Table  4.12.  The  ED-WOW  consists  of  two  cylinders  and  its  rotational  speed  was  about 
120  RPM  (2  Hz)  during  the  experiments.  In  addition,  it  is  assumed  that  the  CO2  leakage 
occurred  only  during  the  intake  and  expansion  processes  because  the  cylinder  pressure 
needs  to  be  higher  than  the  pressure  in  the  crank-housing.  Therefore,  the  leakage  area 
becomes  half  of  the  calculated  flow  rate  as  shown  in  the  4th  row  of  Table  4.12. 

If  it  is  assumed  that  the  end  tips  of  the  piston  ring  does  not  touch  and  there  is  a 
gap  between  the  ends  of  piston  ring,  the  size  of  the  gap  becomes  0.0685  mm  when  pi  is  7 
MPa  and  vapor  quality  is  0.9.  The  measured  width  of  the  piston  ring  is  1.067  mm  as 
shown  in  Table  4.5.  The  calculated  sizes  of  the  leakage  gap  are  shown  in  the  5th  row  of 
Table  4.12.  As  shown  in  the  table,  the  range  of  the  size  of  the  leakage  gap  is  0.0685  mm 
to  0.2158  mm. 

In  addition,  another  situation  is  considered.  If  the  piston  ring  does  not  touch  the 
cylinder  wall,  there  is  a  gap  between  the  cylinder  wall  and  the  piston  ring.  This 
assumption  is  the  opposite  of  the  results  of  the  simulation  because  the  piston  ring  was 
predicted  to  always  touch  the  cylinder  wall  based  on  the  simulation.  However,  this 
assumption  is  considered  in  order  to  get  some  idea  about  the  leakage  area.  The  calculated 
sizes  of  the  gap  are  shown  in  6th  row  of  Table  4.12.  The  range  of  the  size  of  the  gap 
between  the  cylinder  wall  and  the  piston  is  8.4  x  1(T4  mm  to  2.65  x  10 3  mm. 

The  gap  sizes  of  the  leakage  path  seem  to  be  possible  values  as  shown  in  5th  and 
6th  row  of  Table  12.  Therefore,  it  can  be  deduced  that  the  predicted  leakage  flow  rate  of 
CO2  through  the  piston  ring,  31.3%  to  39.6%  of  CO2  flow  rate,  is  reasonable. 

As  shown  in  Table  4.11,  the  work  produced  through  the  ED-WOW  increases  with 
the  leakage  of  CO2.  This  is  due  to  the  low  pressure  during  the  exhaust  process.  If  the 
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exhaust  pressure  is  small,  there  is  less  CO2  mass  remained  inside  the  cylinder  and  less 
work  is  required  to  sweep  out  the  expanded  CO2  from  the  cylinder  to  the  inlet  of  the 
evaporator.  This  can  be  seen  in  Figure  4.34  and  Figure  4.35.  Figure  4.34  shows  the  work 
produced  through  the  ED-WOW  with  the  assumption  of  zero  leakage.  Figure  4.35 
illustrates  the  work  with  some  leakage  of  C02  through  the  piston  ring.  The  negative 
work  stands  for  the  work  produced  by  the  ED-WOW  and  the  positive  work  indicates  the 
work  required  to  sweep  out  the  expanded  CO2  from  the  ED-WOW.  As  shown  in  the 
figures,  the  work  produced  during  the  intake  and  expansion  processes  are  very  similar 
with  both  assumptions  but  the  works  required  during  the  exhaust  process  are  different. 
As  expected,  less  work  is  required  when  the  pressure  is  low  during  the  exhaust  process 
due  to  the  less  C02  mass.  However,  there  is  energy  loss  due  to  the  C02  leakage  and  the 
isentropic  efficiency  of  the  ED-WOW  decreases  as  shown  in  Table  4.11.  If  there  is  no 
leakage,  less  amount  of  C02  is  required  to  achieve  the  desired  low  pressure  at  the  end  of 
the  expansion  process  and  to  obtain  the  desired  work  output  through  the  ED-WOW. 
Therefore,  some  improvement  of  the  leakage  through  the  piston  ring  should  be 
considered  to  increase  the  efficiency  of  the  ED-WOW. 

It  can  also  be  seen  from  Table  4.11  that  the  predicted  works  through  the  ED- 
WOW  (12th  row)  are  much  larger  than  ones  in  the  experiments  (5th  row).  This  means  that 
the  current  method  to  measure  the  work  produced  is  very  conservative. 

Figure  4.36  shows  the  minimum  film  thickness  of  the  lubricant  under  the  piston 
ring  with  the  assumption  of  some  leakage.  As  shown  in  this  figure,  its  value  is  too  small 
and  constant  over  the  whole  processes.  This  means  that  there  is  always  boundary 
lubrication  for  all  cases.  That  is,  the  piston  ring  always  touches  the  cylinder  wall. 
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Figure  4.37  shows  the  predicted  power  loss  due  to  the  friction  between  the  piston 
ring  and  the  cylinder  wall  with  the  assumption  of  some  leakage.  As  shown  in  the  figure, 
the  power  loss  becomes  the  largest  when  the  pressure  inside  the  cylinder  is  the  highest. 
Case  3.  This  is  due  to  the  fact  that  when  the  pressure  is  high,  there  is  much  force  to  push 
the  piston  ring  in  radial  direction  to  the  cylinder  wall  and  it  results  in  increase  of  the 
friction  loss.  The  calculated  values  of  the  power  loss  are  7.16  W,  7.24  W,  and  8.44  W  for 
Case  1,  Case  2,  and  Case  3,  respectively. 
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Table  4.1:  Total  compression  works  of  the  Economizer  cycle  for  low-side  pressures. 


Low-side  pressure 
(MPa) 

Pressure  ratio  at  which  the 
maximum  COP  occurs 

Compression  work 
(kJ/kg) 

2.5 

1.84 

48.383 

3.0 

1.70 

42.017 

3.5 

1.56 

33.249 

4.0 

1.46  _ 

27.330 

4.5 

1.38 

22.575 

Table  4.2:  Maximum  COP  values,  pressure  ratio  and  high-side  pressure  at  which  the 

maximum  COP  occurs.  _ 


Case 

Tj i  and  pn 

Air 

temperature 

TairCQ 

Pressure  ratio 
across  1st 
compressor 

High-side 

pressure 

Pl4 

(MPa) 

COP 

1 

7}i=26.7°C 

/j/i=3.5MPa 

27.8 

1.64 

8.0 

5.457 

2 

35 

1.78 

8.8 

4.081 

3 

40.6 

2.06 

10.4 

3.254 

4 

46.1 

2.32 

12.0 

2.702 

5 

51.7 

2.34 

12.0 

2.187 

6 

7/,=:  32.2°C 
p/i=3.5MPa 

35 

1.72 

8.8 

4.146 

7 

40.6 

2.0 

10.4 

3.307 

8 

46.1 

2.26 

12.0 

2.749 

9 

51.7 

2.26 

12.0 

2.246 
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Table  4.3:  Results  of  thermodynamic  analysis  for  several  CO2  cycles 


perating  condition 


'  at  evaporator  (°C) 


1  at  outlet  of  gas  cooler  (°C) 


'  at  inlet  of  comp.  (°C)  (Cases  3  &  4) 


P  at  evaporator  (MPa) 


P  at  gas  cooler  (MPa) 


P  intermediate  (MPa)  (Cases  5  &  6) 


esults 


ase  cycle  (Case  1) 


fork  of  comp,  (kj/kg) 


Heat  cap.  of  evap  (kJ/kg) 


30P 


ase  +  ED- WOW  (Case  2) 


Heat  cap.  of  evap  (kJ/kg) 


fork  output  throu.  ED-WOW  (kJ/kg) 


30P  without  using  W  of  WEOD 


^OP  using  W  of  WEOD  to  run  comp. 


IIX  Cycle  (Case  3) 


fork  of  comp,  (kj/kg) 


Heat  cap.  of  evap  (kJ/kg) 


30P 


HX  +  ED-WOW  (Case  4) 


Heat  cap.  of  evap  (kJ/kg) 


ork  output  throu.  ED-WOW  (kJ/kg) 


ntcrcooler  cycle  (Case  5) 


fork  of  1-st  comp.  (kJ/kg) 


fork  of  2-nd  comp.  (kJ/kg) 


fork  of  total  comp  (kJ/kg) 


Heat  cap.  of  evap  (kJ/kg) 


30P 


ntercooler  +  ED-WOW  (Case  6) 


eat  cap.  of  evap  (kJ/kg) 


fork  output  throu.  ED-WOW  (kJ/kg) 


20P  without  using  W  of  WEOD 


30P  using  W  of  WEOD  to  run  comp. 


'  at  outlet  of  intercooler  (°C) 


iconomizer  cycle  (Case  7) 


fork  of  1  -st  comp.  (kJ/kg) 


fork  of  2-nd  comp.  (kJ/kg) 


fork  of  total  comp  (kJ/kg) 


Heat  cap.  of  evap  (kJ/kg) 


:op 


uality  at  accumulator 


20P  without  using  W  of  WEOD 

2.34 

30P  using  W  of  WEOD  to  run  comp. 

2.49 

56.058 


2. 


58.161 


143.518; 


2.46 


OC-3 

BISS 

10 

15 

35 

35 

27.5 

29 

4.502 

5.087 

150.04 

3 

145.742 

6.528 

5.636 

2.391 

2.66 

43.738 


135.508 


3.09 


140.325: 


4.817 


3.20 


3.605 


54.393 


168.52 


3.09 


49.017 


140.1 


2.85 


38.69 


1.772 


40.463 


135.508 


3.34 


OC-7 

OC-8 

15 

15 

50 

55 

39.5 

43 

5.087 

5.087 

35.098 


108.115 


35.098 

35.098 

43.619 

17.208 

1.243 

0.49 

43.552 

45.17 

160.767 

145.597 

3.691 

3.223 

48.142 

49.526 

92.114 

70.861 

1.913 

1.431 

149.195 


3.598 


3.303 


3.58 


126.062 


4.325 


>  30.883 

30.883 

30.883 

30.883 

!  1.588 

1.867 

2.556 

3.223 

I  32.471 

32.749 

33.438 

34.106 

1  129.27 

108.115 

43.619 

•  3.981 

3.301 

2.351 

1.279 

133.339 

51.52 

4.069 

4.839 

6.13 

7.901 

35.02 

21.021 

106.772 

47.782 

3.049 

2.273 

0.025 

0.563 

Table  4.4:  Results  of  thermodynamic  analysis  for  intercooler  and  economizer  cycles  to 


find  pressure  ratios  across  com] 

pressors  that  give  maximum  system  performanct 

Operating  condition 

OC-1 

OC-2 

OC-3 

OC-4 

OC-5 

OC-6 

OC-7 

OC-8 

T  at  evaporator  ( C) 

0 

5 

10 

15 

15 

15 

15 

15 

T  at  outlet  of  gas  cooler  ( C) 

35 

35 

35 

ME 

40 

45 

50 

55 

P  at  evaporator  (MPa) 

3.485 

MmMk 

5.087 

5.087 

5.087 

5.087 

5.087 

P  at  gas  cooler  (MPa) 

Results 


Untercooler  cycle 


p-ratio  across  1st  compressor 

2.34 

2.06 

1.82 

1.61 

1.74 

1.78 

1.82 

1.88 

p-ratio  across  2nd  compressor 

mme 

1.26 

W&BEL 

mmm 

Intermediate  pressure  (MPa) 

8.155 

8.177 

mmi 

8.19 

8.851 

- — - 

9.055 

22 

Work  of  1st  compressor  (kJ/kg) 

47.61 

38.773 

30.533 

22.804 

26.803 

27.991 

HBaBj 

30.882 

Work  of  2nd  compressor  (kJ/kg) 

5.707 

5.554 

5.448 

5.47 

4.613 

5.117 

ROSS 

3.701 

Total  compression  work  (kJ/kg) 

53.317 

44.327 

35.981 

28.273 

31.417 

33.108 

34.583 

Without  ED- WOW  (Case  8) 

Heat  cap.  of  evaporator  (kJ/kg) 

143.518 

140.106 

135.508 

129.27 

108.115 

78.614 

43.619 

17.208 

COP 

2.692 

3.161 

3.766 

:$EZ 

3.441 

mmSL 

M m 

|  With  ED- WOW  (Case  9) 

Heat  cap.  of  evaporator  (kJ/kg) 

150.046 

145.742 

140.325 

133.339 

112.954 

84.744 

51.52 

26.624 

Work  output  throu.  ED-WOW  (kJ/kg) 

6.528 

5.636 

4.817 

4.069 

4.839 

6.13 

HE 

COP  without  using  W  of  ED-WOW 

MSB 

3.288 

3.9 

4.716 

HSl^ 

COP  using  W  to  run  compressor 

■US 

mam* 

4.503 

5.509 

3.141 

1.058 

[Economizer  cycle  (Case  10) 


p-ratio  across  1st  compressor 
p-ratio  across  2nd  compressor 
Intermediate  pressure  (MPa) 
Work  of  1st  compressor  (kJ/kg) 
Work  of  2nd  compressor  (kJ/kg) 
Total  compression  work  (kJ/kg) 
Heat  cap.  of  evaporator  (kJ/kg) 
COP 


1.48 

2 

5.158 

15.367 

41.079 

56.447 

139.556 

2.472 


1.36 

1.92 

5.398 

11.931 

36.588 

48.519 

136.546 

2.814 


1.261 


1.82 


5.673 


8.869 


32.011 


40.88 


132.49S 


1.17 

1.74 

5.952 

5.919 

27.597 

33.516 

127.023 

3.79 


1.24 

1.64 

6.308 

7.017 

24.133 

31.15 

102.401 

3.287 


1.28 

1.59 

6.511 

5.426 

21.144 

26.57 

65.949 


1.11 

1.83 

5.647 

1.11 

28.58 

29.691 

37.877 


1.01 

2.01 

5.138 

0.043 

34.434 

34.477 

16.655 


3.241 


2.48: 


1.271 
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Table  4.5:  Input  values  for  simulation  of  the  expansion  device  producing  work. 


Parameters 

Values 

Unit 

Parameters 

Values 

Unit 

Engine  speed 

S  =  120,  400, 
or  600 

RPM 

Volume  at  TDC 

Vtdc  =  1.685 

mm3 

B  =  27.7 

mm 

Fittings’  volume 

msssm 

Bore  radius 

r=  13.85 

mm 

Location  of 
expansion 

LOE  =  8.01 
(5.21  from 
TDC) 

mm 

Crank  radius 

R  =  11.0 

mm 

Diameter  of 
fitting 

Dfit  —  2.286 

mm 

Connecting  rod 
length 

L  =  38.8 

mm 

Orifice  diameter 
of  intake  valve 

Dint-  1.191 

mm 

Composite 

roughness 

<7=0.37 

pm 

Orifice  diameter 
of  exhaust  valve 

Dexh  =  2.286 

mm 

Piston  ring  width 

b=  1.067 

mm 

Number  of  intake 
valves  per 
cylinder 

nin  =  2 

- 

Ring  thickness 

tr=  1.067 

mm 

Number  of 
exhaust  valves 
per  cylinder 

flexh  =  1 

- 

tmdmm 

E  =  70.0 

KS2SZG3I 

Pressure  at  intake 

W3B. ESEEH 

MPa 

Ring  crown 
height 

Cr  =  14.9 

pm 

Desired  pressure 
at  exhaust 

Pexh  =  3.4474 

MPa 

Friction 

coefficient 

/=  0.08 

- 

Desired  pressure 
in  crank  housing 

pcrk=  3.4474 

MPa 

Lubricant 
viscosity 
(C02  at  300  K) 

//  = 

1.50  xl0~5 

Pa-s 

Temperature  at 
intake 

Tin  =  322.0 

K 

Tangential  ring 
tension 

T=  22.38 

N 

Desired 
temperature  at 
exhaust 

Tea  =  298.25 

K 

Volume  at  BDC 

Vbdc  —  14.952 

mm3 

Time  step 

At  =  0.00125, 
0.0005,  or 
0.00025 

Sec. 
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Table  4.6:  Test  conditions  to  investigate  the  performance  of  the  ED-WOW. 


Case 

1 

2 

. _ _ _ 1 

3 

Tjndoor  (  Q 

20 

20 

20 

Toutdoor  (  C) 

35 

35 

35 

Pc, out  (kPa) 

6620.7 

7068.7 

7971.3 

Pc, in  (kPa) 

2770.9 

2637.3 

2528.2 

P gc,in  (kPa) 

6603.5 

7045.3 

7938.2 

Pgc.out  (kPa) 

6536.9 

6931.5 

7845.2" 

PED,in(kPa) 

6052.1 

6392.4 

7255.7 

PED,out  (kPa) 

2921.7 

2814.3 

2759.8 

Pevap.out  (kPa) 

2854.8 

2735.7 

2657.4 

Tc,out(°C) 

87.6 

99.4 

117.9 

Tc,in  (°C) 

20.2 

19.8 

19.6 

Tgc.in  (°C) 

68.0 

77.6 

93.7 

Tgc.out  (°C) 

36.0 

36.2 

36.7 

TED,in  (°C) 

26.0 

27.9 

31.4 

Tevap.out  (  C) 

12.5 

11.0 

10.3 

mCo2  (g/s) 

9.7 

11.6 

13.5 

Wed  (RPM) 

114 

120 

120 

Atjntake  (s) 

0.129 

0.123 

0.123 
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Table  4.7:  Results  of  the  experiments  with  ED-WOW  in  a  transcritical  CO2  cycle. 


Case 

1 

2 

3 

PED,in  (kPa) 

6052.1 

6392.4 

7255.7 

Ped.oui  (kPa) 

2921.7 

2814.3 

2759.8 

Ted, in  (°C) 

26.0 

27.9 

31.4 

Qisenlh  (W) 

425.6 

570.7 

1030 

Qeo  (W) 

448.7 

599.4 

1064 

wcomp  (W) 

567.6 

773.1 

1110 

w£D  (W) 

23.18 

28.70 

34.77 

COPexv 

0.7498 

0.7382 

0.9279 

COPed-i 

0.7905 

0.7753 

0.9586 

COPed-2 

0.8242 

0.8052 

0.9896 

£s,ed  (%) 

10.47 

9.97 

10.27 

Table  4.8:  Increases  of  the  evaporator  capacity  and  COP. 


Case 

1 

2 

3 

Qisen'h  <W) 

425.6 

570.7 

1030 

<2£D  (W) 

448.7 

599.4 

1064 

Increase  of  evaporator 
capacity  (%) 

5.43 

5.03 

3.30 

COPexv 

0.7498 

0.7382 

0.9277 

COPed- 1 

0.7905 

0.7753 

0.9586 

5.43 

5.03 

3.30 

COPed-2 

0.8242 

0.8052 

0.9896 

Increase  of  COP  with 
enthalpy  &  work  output 
effects  (%) 

9.92 

9.08 

6.67 
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Table  4.9:  Uncertainties  of  measurements. 


Variable 

Uncertainty  value 

Flow  meter 

0.5%  of  mass  flow  rate 

T-type  thermocouple 

0.5°C 

Pressure  transducer 

43.1  kPa 

Vacuum  pressure  gage 

4.1  kPa(0.6psi) 

Static  pressure  gage 

13.8  kPa  (2  psi) 

Rotameter 

0.09  nfThr  (0.4  GPM) 

Torquemeter 

0.085  N-m 

Tachometer 

50(1)  RPM 

Note:  U)  The  uncertainty  value  of  the  Tachometer  is  assumed. 


Table  4.10:  Calculated  Results  and  associating  uncertainties  (Case  3). 


Item 

Measure  or 

calculated  value 

Absolute 

uncertainty 

Relative 

uncertainty 

Work  generated 

through  the  ED-WOW 

34.77  W 

10.69  W 

30.74  % 

Isentropic  efficiency 

of  the  ED-WOW 

10.27  % 

3.16  % 

30.79  % 

Evaporator  capacity 

1030  W 

5.43  W 

0.5% 

Compression  work 

mow 

44.23  W 

3.98  % 

COPed-2 

0.9896 

0.004972 

0.5  % 

160 


Table  4.11:  Results  of  the  ED-WOW  simulation  for  the  experimental  test  conditions. 


Case 

1 

2 

3 

PED.in  (kPa) 

6052.1 

6392.4 

7255.7 

PED.out  (kPa) 

2921.7 

2814.3 

2759.8 

Experiment 

™Ca2  (g/s) 

9.7 

11.58 

13.48 

WED  (W) 

23.18 

28.70 

34.77 

Ped.oui  (kPa) 

3657.4 

3889.5 

4569.3 

™COl  (g/s) 

6.16 

6.80 

9.0 

/jvI  U 

leakage 

WED  (W) 

94.28 

98.01 

104.53 

Isentropic 
efficiency  (%) 

56.9 

53.9 

46.0 

PED.out  (kPa) 

2710 

2589.9 

3546.5 

Simulation 

™Co2  (g/s) 

9.71 

11.54 

13.45 

W£D  (W) 

135.6 

154.4 

149.2 

Some 

™C02-leak  (g/s) 

3.27 

4.57 

4.21 

Rate  of  leakage 
(%) 

33.67 

39.6 

31.28 

Isentropic 
efficiency  (%) 

33.7 

28.8 

29.2 

Table  4.12:  Calculation  of  leakage  area  of  CO2  through  the  piston  ring  for  Case  3. 


Parameters 

Vapor  quality,  x  =  0.9 

Vapor  quality,  x  =  0.5 

pi  =  7  MPa 

Pi  =  7  MPa 

Pi  =  4  MPa 

Mass  velocity,  G, 
(kg/s-m2) 

28718.3 

11724.2 

22337.4 

9119.2 

Flow  area 
(mm2) 

0.146248 

0.358233 

0.188025 

0.460566 

Leakage  area 
(mm2) 

0.073124 

0.179116 

0.094013 

0.230283 

Gap  between  the  ends  of 
the  piston  ring 
(mm) 

0.0685 

0.1679 

0.0881 

0.2158 

Gap  between  the 
cylinder  wall  and  the 
piston  ring 
_ (mm) _ 


0.00084 


0.00206 


0.00108 


0.00265 


Jj 


Figure  4.1:  Comparison  of  COPs  of  Economizer  cycle  and  Base  cycle. 


Figure  4.2:  Comparison  of  COPs  of  Intercooler  cycle  and  Base  cycle. 
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Figure  4.3:  COP  and  compression  work  of  Intercooler  cycle. 


Intermediate  pressure  (MPa) 

Pn  =  3.5  MPa,  Tn  =  25  C,  T.r  =  35  C  &  P|2  =  1 0  MPa 

Figure  4.4:  COP  of  Intercooler  cycle  with  different  pressure  ratios  across  the  lst-stage 

and  2nd-stage  compressors. 
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Figure  4.5:  Behavior  of  compression  works  of  Intercooler  cycle  over  intermediate 

pressure. 


Figure  4.6:  Processes  of  Intercooler  cycle  on  pressure-enthalpy  diagram. 


Pressure  ratio  across  compressors 


Figure  4.8:  Effect  of  pressure  at  the  inlet  to  the  evaporator  on  heat  removal  capacity  of 

Economizer  cycle. 


- -  Pn=2.5  MPa,  T=-12.0C 

- P,=3.0  MPa,  T  =-5.6C 

- P,  =3.5  MPa,  T  =0.2C 

- P„=4  0  MPa,  Ts=5.3C 

- P,'=4.5  MPa,  T  =10.0C 


Pressure  ratio 
T  =25C  &  T  =35C 


Figure  4.9:  Effect  of  pressure  at  the  inlet  to  the  evaporator  on  COP  of  Intercooler  cycle. 


■  Base  cycle 

-  Intercooler-1 

-  Intercooler-2 


High  pressure  (MPa) 

T,  =26.7  C,  Ph=3.5  MPa  &  Talr=27.8  C 

Figure  4.10:  Comparison  of  COPs  of  transcritical  C02  cycles  for  7)i=26.70C 

d  Tnir  =  27.8"  ~ 


High  side  pressure  (MPa) 

T,  =26.7  C,  P,  =3.5  MPa  &  T  =35  C 

II  II  air 


High  side  pressure  (MPa) 

T,  =26.7  C,  P„=3.5  MPa  &  T  .  =40.6  C 


Figure  4.12:  Comparison  of  COPs  of  transcritical  CO2  cycles  for  7>i=26.7°C 

and  Tair  =  40.6°C. 


High  side  pressure  (MPa) 
T,=26.7  C,  P,=3.5  MPa  &  T  =46.1  C 


Figure  4.13:  Comparison  of  COPs  of  transcritical  C02  cycles  for  7n-26.7°C 

and  Tair  =  46.1°C. 


Figure  4.14:  Comparison  of  COPs  of  transcritical  C02  cycles  for  Tn-26.1°C 

and  Tnir  =  51.7°C. 


h  [kJ/kg] 


High  side  pressure  (MPa) 

T  =32.2  C,  P  =3.5  MPa  &  T  . =35  C 

1  ’  1  air 


Figure  4.16:  Comparison  of  COPs  of  transcritical  CO2  cycles  for  7/i=32.2°C 


High  side  pressure  (MPa) 

T,=32.2  C,  P=3.5  MPa  &  T,.  =40.6  C 

Figure  4.17:  Comparison  of  COPs  of  transcritical  C02  cycles  for  Tn 

and  Tair  =  40.6°C. 


High  side  pressure  (MPa) 
T1=32.2  C,  P,=3.5  MPa  &  Talr=46.1  C 


:32.2°C 


Figure  4.18:  Comparison  of  COPs  of  transcritical  C02  cycles  for  7n=32.2°C 

and  T^r  =  46.1°C. 


Sf'°  3.5  8.0  ^ 
T,  =26.7  C.  P,  =3.5  MPa  &  T_=40.6  C 


Figure  4.22:  Behavior 


of  COP  of  Intercooler  cycle  for  Tn=26.7°C  and  ra/r-40.6°C 


3 


a 


a 


3 


[%] 


540  630 

Angle  of  crankshaft  (degree) 

Figure  4.25:  Piston  speed  versus  crank  angle  of  the  ED-WOW. 


Figure  4.26:  Repeat  of  the  same  trends  of  the  pressure  behavior  for  120RPM  and 

240RPM. 
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Intake 


Figure  4.27:  CO2  pressure  in  the  chamber  versus  crank  angle. 


Figure  4.28:  Work  done  by  CO2  versus  crank  angle  (120  RPM). 


Power  loss  (W)  £  Minimum  film  thickness  (urn) 


.29:  Minimum  film  thickness  under  the  piston  ring  (120  RPM). 


Intake 


Figure  4.30:  Power  loss  due  to  friction  of  the  piston  ring  (120  RPM). 


P  [kPa] 
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h [kJ/kg] 

Figure  4.31:  Pressure-enthalpy  diagram  of  the  processes  of  the  system. 
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Angle  of  crankshaft  (Degree) 


Figure  4.32:  Pressure  behavior  inside  the  cylinder  of  the  ED-WOW  with  the  assumption 

of  zero  leakage  of  CO2. 


Pressure  (MPa) 
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Angle  of  crankshaft  (Degree) 


Figure  4.33:  Pressure  behavior  inside  the  cylinder  of  the  ED-WOW  with  the  assumption 

of  leakage  of  CO2  through  piston  rings. 


Work  (J 


0.4 


Figure  4.34:  Predicted  work  produced  with  zero  leakage  assumption. 


Minimum  film  thickness  (|im) 


Power  loss  (W) 
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Angle  of  crankshaft  (Degree) 


5  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


A  study  to  increase  system  performance  of  a  transcritical  CO2  cycle  has  been 
conducted.  Specifically,  an  expansion  device  with  output  of  work  (ED-WOW)  has  been 
designed,  constructed  and  incorporated  into  such  a  cycle.  The  ED-WOW  extracts  energy 
from  CO2  during  the  expansion  process  in  a  transcritical  CO2  cycle  which  can  be  used  to 
drive  the  compression  process.  In  addition,  with  work  extraction  reduced  enthalpy  of  the 
CO2  can  be  achieved  at  the  inlet  of  the  evaporator.  These  two  effects  of  the  ED-WOW 
increase  the  system  performance  of  CO2  cycle  much.  This  final  chapter  will  summarize 
the  tasks  completed  to  meet  the  objectives  of  this  research,  make  conclusions  based  on 
work  done,  and  give  some  suggestions  for  future  direction  of  the  ED-WOW. 


5.1  Summary  of  Results 

Because  of  high  volumetric  heat  capacity  of  CO2,  it  offers  the  potential  for 
reduced  weight  and  volume  in  package  system,  which  is  a  major  focus  for  the  automobile 
and  military  applications.  However,  the  overall  efficiency  of  the  transcritical  C02  cycle 
decreases  as  the  weight  and  volume  of  the  system  decrease.  System  improvements  must 
be  achieved  to  meet  the  goals  of  the  reduction  of  weight  and  volume  and  while  still 
retaining  competitive  cycle  efficiency. 
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5.1.1  Cycle  Analysis  Conclusions 

A  thermodynamic  analysis  was  performed  for  several  modified  transcritical  CO2 
cycles,  with  a  goal  of  determining  the  proper  configuration  to  achieve  the  optimum 
system  performance.  The  analysis  shows  that  the  system  efficiency  of  the  transcritical 
CO2  cycle  can  be  dramatically  increased  by  modifications  of  the  system  configuration, 
and  that  the  optimum  performance  is  a  strong  function  of  the  operating  conditions. 

The  system  efficiency  of  the  transcritical  C02  cycle  is  increased  considerably  by 
nearly  100  %  when  the  low-side  pressure  increases  from  2.5  MPa  to  4.5  MPa.  There  is  a 
specific  high-side  pressure  at  which  the  maximum  COP  occurs  in  the  CO2  cycle.  There 
are  also  specific  pressure  ratios  across  the  compressors  for  the  minimum  compression 
work  and  the  maximum  COP  for  the  two-compressor  system. 

As  previously  discussed,  the  employment  of  a  work  extraction  device  in  the  role 
of  an  expander  is  expected  to  give  the  largest  improvement  of  the  system  efficiency.  This 
device  will  not  only  increase  the  heat  removal  capacity  of  the  evaporator  but  also  can 
reduce  the  compression  work.  The  thermodynamic  analysis  shows  that  system  efficiency 
of  the  transcritical  CO2  cycle  is  increased  by  16%  to  112%  with  the  work  extraction 
device  compared  to  the  cycle  without  the  device. 

5.1.2  Cycle  with  Work  Extraction 

A  major  objective  of  this  research  was  to  design,  construct,  test,  and  analyze  an 
expansion  device  which  provides  for  work  extraction  (the  ED-WOW  device).  Based  on 
the  review  of  literature  and  technologies,  a  piston-cylinder  type  device  was  chosen  as  the 
ED-WOW.  A  theoretical  model  was  developed  to  verify  the  proper  operation  of  the  ED- 
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WOW  in  prototype  application  and  to  identify  loss  mechanisms  of  the  ED-WOW.  This 
model  bases  on  the  lubrication  theory  with  assumption  that  the  piston  ring  is  a  slippery- 
pad  bearing.  This  model  calculates  the  pressure  inside  the  device,  the  work  done  by  CO2, 
the  power  loss,  etc.  After  preliminary  analysis,  the  ED-WOW  was  constructed  and 
installed  in  a  prototype  transcritical  CO2  refrigeration  system,  replacing  the  simple 
throttling  valve.  During  cycle  tests,  the  use  of  the  device  provided  7%  to  10% 
improvement  in  COP.  These  figures  are  considered  to  be  highly  conservative,  as  the 
work  extraction  measurement  was  biased  toward  a  lower  value  due  to  uncertainty  in  the 
load  device  performance. 


5.2  Recommendations 

Recommendations  based  on  this  study  are  given  for  the  ED-WOW.  Some 
improvements  of  the  design  of  the  ED-WOW  can  be  made  to  achieve  better  performance. 


5.2.1  Design  Improvements 

Presented  below  are  some  suggestions  to  achieve  better  performance  or  optimize 
the  design  of  the  ED-WOW,  which  may  be  topic  of  future  research  using  the  ED-WOW: 

•  The  rotational  speed  of  the  ED-WOW  needs  to  be  increased,  reducing  the  torque 
loading  on  the  device.  This  can  be  achieved  with  larger  intake  and  exhaust  ports, 
and  either  a  mechanically  driven  valve  train  or  faster-acting  solenoid  valves. 

•  Improvement  of  the  piston  ring  reducing  the  leakage  of  CO2  through  it  should  be 


considered. 
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•  A  hermetically  sealed  device  should  be  constructed,  eliminating  the  need  for  a 
large  enclosure  and  reducing  the  device  weight. 

•  Heat  transfer  needs  to  be  considered  in  the  simulation  model.  The  current 
simulation  model  does  not  have  the  capacity  to  calculate  the  heat  transfer  between 
the  C02  inside  the  device  and  the  atmosphere  through  the  body  of  the  ED-WOW. 
As  the  C02  expands  from  high  pressure  and  temperature  inside  the  device,  the 
temperature  of  C02  decreases  and  it  becomes  much  lower  than  the  atmospheric 
temperature.  Therefore,  heat  will  be  transferred  from  the  atmosphere  to  the  C02 


inside  the  device. 
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APPENDIX  A:  Reynolds  Equation 


The  simplified  equations  appropriate  for  lubrication  theory  are  known  as 
Reynolds  equation.  Figure  A.l  shows  flow  of  a  lubricant  under  a  piston  ring.  For 
simplicity  let  the  piston  ring  have  only  a  vertical  velocity,  V(t),  which  may  be  a  function 
of  time,  and  the  lower  wall  (cylinder  wall)  a  steady  horizontal  velocity,  U.  Assumptions 
are  as  follows: 

0 

•  The  flow  is  axisymmetric.  That  means  —  =  0 . 

aO 

•  Diameter  of  the  cylinder  is  very  large  compared  to  the  thickness  of  the  piston  ring 
so  that  x-y  coordinate  can  be  applied. 

•  The  origin  of  the  coordinate  is  fixed  as  shown  in  the  figure. 

•  The  width  of  the  piston  ring,  b,  is  very  large  compared  to  the  film  thickness  of  the 

h 

lubricant,  h,  so  that  —  0 . 

b 

•  Body  force  of  the  lubricant  is  negligible. 

•  The  lubricant  behaves  as  a  Newtonian  fluid. 

•  Viscosity  is  uniform  over  the  ring  thickness. 


The  following  scales  are  appropriate: 


x~b 

u  ~  U  =  Ro) 
y~K 


(A-l) 


R  and  (0  stand  for  crank  radius  and  revolutionary  speed  in  the  equation  above, 


respectively.  Continuity  equation  is 


du  dv  „ 
—  +  —  =  0. 
dx  9 y 


(A-2) 


The  terms  have  the  following  sizes  (omitting  the  order  symbol): 


W«0. 

b  K 


(A-3) 


The  size  of  the  y-component  of  the  velocity  becomes 


The  x-component  of  the  momentum  equation  is 


(A-4) 


du  du  dp  d2u  d2u 


(A-5) 


The  sizes  of  the  terms  (omitting  the  order  symbol)  are 

U  h  U  U  dp  U  U 
pU  —  +  p +  M  VT  +  M  tt 
b  b  R  ^  b2  K 


(A-6) 


Reorganizing  above  equation  gives 


U 2  U2 

P-r  +  P~r  = 

b  b 


dp  U  U 

- —  +  U— T  +  p— 7 

dx  ^b2  Hh2 


(A-l) 


fl 

By  multiplying  above  equation  by  h„  and  considering  the  limit  -»  0 ,  the  following 


equation  is  obtained: 
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_  dp  d2u 


(A-8) 


The  y-component  of  the  momentum  equation  is 


3v  dv  dp  d2v  d2v 

/*3?+^=_^+/'a?+v 


(A-9) 


The  sizes  of  the  terms  (omitting  the  order  symbol)  for  above  equation  are 


UK  UK  UK  UK 

rr  h  KU  h  dp  h  h 

PU-J-  +  P  -,  -T-  =  ~T+v-$r+M-rr 
b  b  K  dy  b  K 


(A- 10) 


Reorganizing  above  equation  gives 


U2  K 


U2  K 


p - - +  p - -  = 

b  b  b  b 


dp  UK  u  K 

dy  ^  b 2  b  *  hi  b 


(A- 11) 


h  h 

Each  term  of  Equation  (A-ll)  contains  — .  By  considering  the  limit  —  — »  0 ,  Equation 

b  b 

(A-l  1)  can  be  reduced  as  follow: 


&  =  0. 


(A- 12) 


Meaning  of  Equation  (A-12)  is  that  there  is  no  pressure  gradient  in  y-direction. 
Boundary  conditions  for  the  flow  are 

u(y  =  0)  =  U 
u(y  =  h(t))  =  0 


(A- 13) 


Partial  integration  of  Equation  (A-8)  over  y  produces 


"=ily!+Ciy+c- 


(A- 14) 


Applying  the  boundary  conditions  gives 
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Integrating  the  continuity  equation  (Equation  (A-2))  over  y  gives 

hr  du  hrdv 


Rearranging  Equation  (A-16)  produces 

J-|ii  ^  -  ;|H  rfy  =  -[v(A  )  -  v(0)]  =  -V  (,  ) . 


(A- 15) 


(A-16) 


(A- 17) 


From  Leibnitz’s  theorem  [Panton,  1996] 


9  ‘f  hrdu  dh  \  dO  /  \ 
-judy=l-dy  +  -u(h)-~u(0). 


0  0 


(A-18) 


The  2nd  and  3rd  terms  of  Equation  (A-18)  are  zero  because  u(h)  =  Oand— • -  -  0  . 


Substituting  Equation  (A-17)  into  Equation  (A-18)  produces 


—  \udy=  J-r-dy  =  -V(t) . 
dx  J  J  dx 


(A-19) 


Integration  of  Equation  (A- 15)  gives 


uyJ±^(^->Lh)+u(y-yL)]=-LaJL 

J  J  [ 2judx{3  2  )  {  2h)\o  2 ju  dx 


*1 U* 
6  2 


h 3  dp  i  Uh 
12  fj.  dx  2 


By  substituting  Equation  (A-20),  Equation  (A-19)  becomes 


d  ' 

h 3  dp  |  Uh 

dx  _ 

\2p  dx  2 

(A-20) 


(A-21) 
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Reorganizing  Equation  (A-21)  produces  the  Reynolds  equation  for  lubrication  theory  as 
follows: 


±9r*>^  =  «,*+Iar«. 

fl  ox  \  ox )  ox 


(A-22) 
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Figure  A.  1:  Flow  under  a  piston  ring. 
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APPENDIX  B:  Analytical  Integrations  of  the  Terms  in  the  Lubrication  Theory 


Results  for  the  analytical  integrations  of  the  terms  in  Equations  (2-56)  and  (2-57) 
are  presented  here.  The  profile  of  the  face  of  the  piston  ring  is  expressed  as: 


K(x)  = 


4  C,2 


x  . 


(B-l) 


The  mathematical  expression  for  the  film  thickness  is 

*(*»*)  =  + 


4Cfi  . 


is  a  constant  and  let  it  be  k  here  for  manipulation  convenience. 


(B-2) 


I,  =  l~dx  =  f7 - ! — T 

>hl  !k*^2) 


-1 


tan 


ldX  2'’.k+fa2J+  2 


I2~  f —dx=  - — — ^dx= — 7 - - - 77 

>hi  ikk+te2) 


3  tan 


-i 


(B-3) 


(B-4) 


j  _  r_l_  ,  _  r  1  ,  _ x _  3x  _ [  V 

3"J'-3  _J k+kx ~  4 hmk+kx2f  +  Shl(h,  +fa2)+  8 4lh% 


(B-5) 
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APPENDIX  C:  Program  to  Control  Solenoid  Valves 


This  program  has  been  developed  to  control  the  solenoid  valves.  There  are  four 
groups  of  the  valves  which  should  be  controlled  separately.  Those  are  1st  Intake  Valve, 
1st  Exhaust  Valve,  2nd  Intake  Valve,  and  2nd  Exhaust  Valve,  respectively.  The  revolution 
speed  of  the  ED-WOW  is  assumed  to  be  120  RPM. 

The  encoder  generates  2048  pulses  per  revolution.  The  intake  valves  open  10 
pulses  after  the  top  dead  center  (TDC)  of  the  piston  and  for  60  pulses.  The  opening  time 
is  very  short  because  the  time  delay  for  the  input  signal  of  closing  is  about  0.053  s.  The 
exhaust  valves  are  closed  10  pulses  after  TDC  and  open  10  pulses  before  the  bottom  dead 
center  (BDC)  of  the  piston.  This  means  that  the  exhaust  valves  are  closed  for  1004 
pulses. 

The  time  delay  is  considered  only  for  the  closing  of  the  intake  valves  in  the 
program.  The  time  delays  for  the  opening  of  all  the  valves  and  the  closing  of  the  exhaust 
valves  are  not  considered. 

The  program  below  shows  the  source  program  written  in  ‘Visual  C++’  to  control 


the  solenoid  valves  based  on  the  discussion  above. 
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/********************************************************************* 


*  Program: 

*  Control.c 

* 


*  Description: 

*  Generation  of  a  retriggerable  delayed  pulse  after  each  transition 

*  on  the  gate  input. 

* 

*  All  valves  (Intake  and  Exhaust  valves)  are  considered  at  the  same  time. 

* 


*  The  intake  valves  are  open  10  pulses  after  TDC. 

* 

*  The  exhaust  valves  are  open  10  pulses  before  TDC  and  close 

*  10  pulses  after  TDC. 

* 

*  Encoder  generates  2048  pulses  per  revolution.  1024  pulses  indicates 

*  1  stroke  of  piston  (down  or  up). 

* 

*  ***  This  program  was  created  by  J.  S.  Baek  on  November  13,  2001.  *** 

* 


* 

*********************************************************************/ 


/* 

*  Includes: 
*/ 


#include  "nidaqex.h" 


/* 

*  Main: 

*/ 

void  main(void) 

{ 

/* 

*  Local  Variable  Declarations: 

*/ 

/*  1st  Intake  Valves  */ 

u32  ulLOWcountO  =  10;  //  Valves  are  closed  for  10  pulses  after  TDC. 

u32  ulHIGHcountO  =  60;  //  Valves  are  open  for  the  next  60  pulses  and  then  closed. 

/*  1  st  Exhaust  Valve  */ 

u32  ulLOWcountl  =  10;  //  Valves  are  open  for  10  pulses  after  TDC. 

u32  ulHIGHcountl  =  1004;  //  Valves  are  closes  for  the  next  1004  pulses  and  then  open. 

/*  2nd  Intake  Valves  */ 

u32  ulLOWcount2  =  1034;  //  Valves  are  closed  for  1034  pulses  after  TDC. 

u32  ulHIGHcount2  =  60;  //  Valves  are  open  for  the  next  60  pulses  and  then  closed. 

/*  2nd  Exhaust  Valve  */ 

u32  ulLOWcount3  =  1034;  //  Valves  are  open  for  1034  pulses  after  TDC. 

u32  ulHIGHcount3  =  1004;  //  Valves  are  closes  for  the  next  1004  pulses  and  then  open. 

il6  iStatus  =  0; 
il6  iRetVal  =  0; 
il6  iDevice  =  1; 
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u32  ulGpctrNumO  =  ND_COUNTER_0; 
u32  ulGpctrOutputO  =  ND_GPCTRO_OUTPUT; 
u32  ulGpctrNuml  =  ND_COUNTER_l ; 
u32  ulGpctrOutputl  =  ND_GPCTRl_OUTPUT; 
u32  ulGpctrNum2  =  ND_COUNTER_2; 
u32  ulGpctrOutput2  =  ND_GPCTR2_OUTPUT; 
u32  uIGpctrNum3  =  ND_COUNTER_3 ; 
u32  ulGpctrOutput3  =  ND_GPCTR3_OUTPUT; 

il6  ilgnoreWaming  =  0; 
il6  iYieldON  =  1; 


/*  . .  1  st  Intake  valves - */ 

iStatus  =  GPCTR_Control(iDevice,  ulGpctrNumO,  ND_RESET); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Control/RESET",  ilgnoreWaming); 

iStatus  =  GPCTR_Set_Application(iDevice,  ulGpctrNumO,  ND_RETRIG_PULSE_GNR); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Set_Application",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNumO,  ND_SOURCE,  ND_PFI_39);  //  counter  input  source 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/SOURCE",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNumO,  ND_COUNT_2,  ulHlGHcountO); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNTr',  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNumO,  ND_COUNT_l,  ulLOWcountO); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNT2",  ilgnoreWaming); 

/*  Each  time  a  pulse  arrives  in  the  gate,  a  new  value  will  be 
latched  into  the  counter  and  sent  to  the  data  buffer.  */ 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNumO,  ND_GATE,  ND_PF1_38); 
iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/GATE",  ilgnoreWaming); 

/*  To  output  a  counter  pulse,  you  must  call  Select_Signal.  */ 

iStatus  =  Select_Signal(iDevice,  ulGpctrOutputO,  ulGpctrOutputO,  ND_LOW_TO_HIGH); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "Select_Signal/GpctrOutput",  ilgnoreWaming); 

//  iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNumO,  ND_OUTPUT_POLARITY,  ND_NEGATIVE); 

//  iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/ND_OUTPUT_POLARITY", 
ilgnoreWaming); 


//  -  1st  Exhaust  valve - 

iStatus  =  GPCTR_Control(iDevice,  ulGpctrNuml,  ND_RESET); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Control/RESET",  ilgnoreWaming); 
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iStatus  =  GPCTR_Set_Application(iDevice,  ulGpctrNuml,  ND_RETRIG_PULSE_GNR) ; 
iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Set_Application",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNuml,  ND_SOURCE,  ND_PF1_35);  //  counter  input  source 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/SOURCE",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNuml ,  ND_COUNT_l ,  ulLOWcountl); 

iRetVal  =  NlDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNTl ",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNuml,  ND_COUNT_2,  ulHIGHcountl); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNT2",  ilgnoreWaming); 

//  Each  time  a  pulse  arrives  in  the  gate,  a  new  value  will  be  latched  into  the  counter  and  sent  to  the  data  buffer. 
iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNuml,  ND_GATE,  ND_PFI_34); 
iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/GATE",  ilgnoreWaming); 

//  To  output  a  counter  pulse,  you  must  call  Select_Signal. 

iStatus  =  Select_Signal(iDevice,  ulGpctrOutputl,  ulGpctrOutputl,  NDJLOW_TO_HIGH); 

iRetVal  =  NIDAQErTorHandler(iStatus,  "Select_Signal/GpctrOutput",  ilgnoreWaming); 

iStatus  =  GPCTR_ChangeJParameter(iDevice,  ulGpctrNuml,  ND_OUTPUT_POLARITY,  ND_NEGATIVE); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/ND_OUTPUT_POLARITY", 
ilgnoreWaming); 

//  . .  2nd  Intake  valves . . 

iStatus  =  GPCTR_Control(iDevice,  ulGpctrNum2,  ND_RESET); 

iRetVal  =  NIDAQEiTorHandler(iStatus,  "GPCTR_Control/RESET" ,  ilgnoreWaming); 

iStatus  =  GPCTR_Set_Application(iDevice,  ulGpctrNum2,  ND_RETRIG_PULSE_GNR) ; 

iRetVal  =  NID AQErrorHandler(i S tatus ,  "GPCTR_Set_Application",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum2,  ND_SOURCE,  ND_PFI_31);  //  counter  input  source 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/SOURCE",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum2,  ND_COUNT_l,  ulLOWcount2); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNTl",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum2,  ND_COUNT_2,  ulHIGHcount2); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNT2",  ilgnoreWaming); 

//  Each  time  a  pulse  arrives  in  the  gate,  a  new  value  will  be  latched  into  the  counter  and  sent  to  the  data  buffer. 
iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum2,  ND_GATE,  ND_PFI_30); 
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iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/GATE",  ilgnoreWaming); 

//  To  output  a  counter  pulse,  you  must  call  Select_Signal. 

iStatus  =  Select_Signal(iDevice,  ulGpctrOutput2,  ulGpctrOutput2,  ND_LOW_TO_HlGH); 

iRetVal  =  NIDAQErrorHandler(iStatus,  ''Select_Signal/GpctrOutput",  ilgnoreWaming); 

//  . —  2nd  Exhaust  valve - - 

iStatus  =  GPCTR_Control(iDevice,  ulGpctrNum3,  ND_RESET); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Control/RESET",  ilgnoreWaming); 

iStatus  =  GPCTR_Set_Application(iDevice,  ulGpctrNum3,  ND_RETR1G_PULSE_GNR); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Set_Application",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum3,  ND_SOURCE,  ND_PFI_27);  //  counter  input  source 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/SOURCE",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum3,  ND_COUNT_2,  ulHIGHcount3); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNTl ",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum3,  ND_COUNT_l ,  ulLOWcount3); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/COUNT2",  ilgnoreWaming); 

//  Each  time  a  pulse  arrives  in  the  gate,  a  new  value  will  be  latched  into  the  counter  and  sent  to  the  data  buffer. 
iStatus  =  GPCTR_Change_Parametcr(iDevice,  ulGpctrNum3,  ND_GATE,  ND_PFI_26); 
iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/GATE”,  ilgnoreWaming); 

//  To  output  a  counter  pulse,  you  must  call  Select_Signal. 

iStatus  =  Select _Signal(iDevice,  ulGpctrOutput3,  ulGpctrOutput3,  ND_LOW_TO_HIGH); 

iRetVal  =  NIDAQErrorHandler(iStatus,  "Select_Signal/GpctrOutput",  ilgnoreWaming); 

iStatus  =  GPCTR_Change_Parameter(iDevice,  ulGpctrNum3,  ND_OUTPUT_POLARITY,  ND_NEG ATI VE) ; 

iRetVal  =  NIDAQErrorHandler(iStatus,  "GPCTR_Change_Parameter/ND_OUTPUT_POLARITY", 
ilgnoreWaming); 

//  - - 

//  NOTE:  On  USB  devices,  this  call  takes  longer  than  normal,  so 
//  you  may  not  be  able  to  output  a  pulse  as  short  as  1ms. 

iStatus  =  GPCTR_Control(iDevice,  ulGpctrNumO,  ND_PROGRAM); 
iStatus  =  GPCTR_Control(iDevice,  ulGpctrNuml,  ND_PROGRAM); 
iStatus  =  GPCTR_Control(iDevice,  ulGpctrNum2,  ND_PROGRAM); 
iStatus  =  GPCTR_Control(iDevice,  ulGpctrNum3,  ND_PROGRAM); 

iRetVal  =  NlDAQErrorHandler(iStatus,  "GPCTR_Control/PROGRAM",  ilgnoreWaming); 


//  HINT:  If  you  don't  see  pulses  at  the  OUTPUT  of  your  counter,  check  your  connections. 
printf("  Hit  any  key  to  stop.\n"); 
iRetVal  =  NIDAQWaitForKey(O.O); 


//Reset  GPCTR. 

iStatus  =  GPCTR_Control(iDevice,  ulGpctrNumO,  ND_RESET); 
iStatus  =  GPCTR_Control(iDevice,  ulGpctrNuml,  ND_RESET); 
iStatus  =  GPCTR_Control(iDevice,  ulGpctrNum2,  ND_RESET); 
iStatus  =  GPCTR_Control(iDevice,  ulGpctrNum3,  ND_RESET); 

//  Disconnect  GPCTRO.OUTPUT. 

//  Note  that  this  will  cause  the  output  to  be  high  impedance  which 
//  will  most  likely  bring  the  logic  level  HIGH  if  there  is  a  pull-up 
//  resistor  on  this  pin.  (Check  your  hardware  user  manual.)  If  you  do 
//  not  want  this  behavior,  comment  out  the  next  line. 

iStatus  =  Select_Signal(iDevice,  ND_GPCTRO_OUTPUT,  ND_NONE,  ND_DONT_CARE) 
iStatus  =  Select_Signal(iDevice,  ND.GPCTRl  .OUTPUT,  ND.NONE,  ND.DONT.CARE) 
iStatus  =  Select_Signal(iDevice,  ND.GPCTR2.0UTPUT,  ND.NONE,  ND_DONT_CARE) 
iStatus  =  Select_Signal(iDevice,  ND_GPCTR3_OUTPUT,  ND.NONE,  ND_DONT_CARE) 

printf("  Single  pulse  has  been  generated  at  GPCTRO_OUTPUT!\n”); 


/*  End  of  program  */ 


